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Abstract 

We carry out a systematic classification and computation of next-to-leading order kinematic 
lO . power corrections to the fully differential cross section in the parton shower. To do this we devise 

L». I a map between ingredients in a parton shower and operators in a traditional effective held theory 

I • I framework using a chain of soft-collinear effective theories. Our approach overcomes several difh- 

f^ ' culties including avoiding double counting and distinguishing approximations that are coordinate 

^P ■ choices from true power corrections. Branching corrections can be classified as hard-scattering, 

that occur near the top of the shower, and jet-structure, that can occur at any point inside it. 

Hard-scattering corrections include matrix elements with additional hard partons, as well as power 
rS ' suppressed contributions to the branching for the leading jet. Jet-structure corrections require 

C^ ' simultaneous consideration of potential 1^-2 and 1^-3 branchings. The interference structure 

induced by collinear terms with subleading powers remains localized in the shower. 
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I. INTRODUCTION 

For scattering problems involving strongly-interacting particles, we are often interested 
in final states with large multiplicities, sometimes including thousands of hadrons. To get to 
this level, we cannot rely solely on full fixed-order calculations. Tree-level event generators 
pj-l^ only go up to 8-10 external particles as Monte Carlo for higher multiplicity phase 
space is increasingly intractable. At one-loop, the frontier is 2 — )■ 4 processes, which have 
been done at the level of differential cross sections for W+ 3 jets p, l6[ and tibb |7|]. At 
two-loops, there are 2 — )■ 1 exclusive calculations for weak boson production by hadrons 



followed by decay (W and Z [8| and W [9i] to leptons, and H decaying to photons |lOl. Illl|.) 



Additionally, e"'"e~ — )■ 3 jets to NNLO is known |12Hl5l|. In any case, a strict fixed order 
counting is not suitable for exclusive observables with large multiplicities, nor for many 
inclusive observables where certain regions of phase space receive kinematic enhancement 
by large logarithms. If Q is a hard scale in the process, then a subset of the amplitude 
gets enhanced so that its coefficient is (asln^((5/p))™, where p <^ Q refers to a small scale 
that is induced by the choice of observable or cuts. Since we can resum these large logs by 
systematically treating real radiation, we can give a leading log (LL) description of these 
observables without performing multiloop computations. The soft and collinear limits that 
yield these large logs also allow us to simplify the amplitude. Therefore, capturing the 
dominant contributions to these observables and simulating processes with a large number 
of particles becomes feasible. This is a main goal of parton Shower Monte Carlo (SMC). 

A final state SMC is based on the "strongly-ordered limit," which describes the leading 
log contribution (accounting for soft emission by angular ordering or other approximations). 
In this kinematic configuration, each radiated particle comes off much more collinear to its 
parent than the previous one, a situation that can be formulated in terms of perpendicular 
momenta or virtualities, i.e. 

qo± > gi± > g2± > • • • , or go^ > g^ > ^2^ > . . . . (1) 

Furthermore, and important for practical computation, in this limit each collinear emission 
is independent of the previous one. Thus, if we have calculated the differential cross section 
for z-parton emission, dai, then we can obtain the {i + l)-parton case as 

p(0) 

dcTi+i oc ^^ dai, (2) 

where P'-^^ is the leading order (LO) "splitting function" that captures the probability for 
the i*'^ emitted parton, of type j, to split into two others, kl, and qf is its virtuality. We 
can therefore formulate the process in terms of a probabilistic Markov chain of i 1 — )■ 2 
particle splittings. The probabilities are determined by the functions P\_^kii which are the 
LO Altarelli-Parisi kernels. As an example, for q — )■ qg, after averaging and summing over 
spins, 

na,-^c.i±^. (3) 

where z is the longitudinal momentum fraction of the daughter with respect to the parent. 
This classical, probabilistic process gives rise to the SMC algorithms used by event generators 



such as P5d;hia [ly, [l7(] and Herwig 18|, ll9| to model radiation. For a virtuality-ordered 
shower, such as the original version of Pythia, given some initial offshellness, q^, and an 
initial momentum fraction, xq, SMCs generate the virtuality and the momentum fraction 
of the daughter particle after the spitting. The former is determined by a Sudakov factor, 
A(g^, gg), which gives the probability of a parton to evolve from ql to q^ without branching. 



H(f^<il) = exp 



dq 
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The traditional LL parton shower makes the multiplicity problem tractable, but it has 
shortcomings related to the leading log approximation. Even though Eq. ([2]) is only correct 
in the collinear limit, the shower is used everywhere in order to generate events that cover the 
full phase space. In addition, since each collinear emission is independent from the previous 
one in the shower, the LL approximation does not include their spin or color correlations, 
nor any of their interference. The situation is different for soft gluons where the inclusion 
of color effects allows one to work in the simplifying limit of angular ordering. 

The hierarchy of scales in the parton shower makes it amenable to an effective field theory 
treatment. Since the shower regime occurs for particles in the soft and collinear regions, we 
can describe it with Soft-Collinear Effective Theory (SCET) |20l-l23l|. Like any EFT, SCET 
comes with an expansion that allows, in principle, for systematic improvement. The first 
work on parton showers using SCET came in 2J, |25|, which we review in Sec. IIIBl where 
the authors showed how the splitting functions and corresponding Sudakov factors, along 
with the factorization of emissions emerge naturally. Furthermore, they could include virtual 
corrections by matching to QCD at higher order in a^. Unfortunately, in reproducing the 
LL shower in SCET, they introduced many conventions whose extension to higher orders in 
the kinematic expansion is unclear. We therefore develop a modified approach to alleviate 
these difficulties. 

Before discussing our setup, we give an overview of advances in the parton shower litera- 
ture beyond the basic LL picture. The structure of these advances depends on what aspect 
of the shower one aims to improve. Possible motivations include accuracy at higher orders 
in as-, higher order in logs, and higher order in powers of the kinematic expansions. We first 
introduce some terminology for higher order log resummation. If the resummation of large 
logs 



, L, is at the cross section level 



da 



j^i^sir 



^^+[j]a,L(a,L2)^' 



+ 



NLL 



(5) 



then we will refer to it as LL, NLL, etc, as indicated. If the cross section transformed to an 
appropriate set of variables has a resummation of logs in the exponent 






+ ... 



NLLc 



(6) 



then we will attach a subscript "exp" to the orders to indicate this. 

A major concern with parton showers is how one handles the merging with matrix element 
(ME) calculations that describe the initial underlying hard process. One can consider a 



simple setup where one declares that a scale, /zq, divides collinear from hard radiation. Here, 
emissions above /io are described through tree-level ME calculations, and those beneath by 
running SMC. Each regime would get a reasonable treatment, but naively interfacing the 
two leaves leading-log sensitivity to /io- This is because the LO (in as) result contains no 
Sudakov log resummation. Methods for carrying out matrix element and parton shower 



merging including this information have been considered in Refs. [26|-|28[ and are referred 
to as CKKW-L and MLM. In CKKW-L, one distributes the particles in an event according 
to the probabilities given by the exact tree- level matrix element, with /Iq as a lower cutoff 
related to the perpendicular momentum between any two particles. One then clusters the 
event using the kx algorithm 29| to determine the splitting scales, qfrp. With these in hand. 



one reweights the event by multiplication by appropriate Sudakov factors, as well as factors 
of as{(liT)/c(s{Q), where Q is some hard scale. We can then run a parton shower algorithm on 
these squared amplitudes, vetoing any splitting qiT harder than /Iq to avoid double counting. 
It was demonstrated that the n-jet rate depends on /ig only beyond NLL order, with the 
first missing term being a"^ In^ {Q / ^q) . CKKW-L has been built into Sherpa j30i] . 

Another important effect concerns soft gluons, which are also kinematically enhanced. 
Collinear emissions reinforce the picture of partonic radiation as an isolated jet since they 
get distributed within some narrow cone about the original hard parton. A priori, soft 
gluons have no preferred direction and can communicate between elements of the shower. 
Fortunately, wide-angle radiation only observes the net color charge contained in the cone 
of emission. Therefore, the pattern of soft radiation far from the collinear jet is not sensitive 
to splittings that have taken place within it. This coherent branching and angular ordering 
can be accommodated by methods such as evolving the shower by decreasing angle mono- 



tonically, as is done in Herwig 3l|, or by enforcing it with a veto in a virtuality-ordered 



shower (the rightmost expression in Eq. ([T])), which is an option in Pythia 16] . Accounting 



for coherence properties leads to LL resummation for the soft emissions [32|-|36| . Additional 
considerations treated in shower programs include putting a^ at the kx scale of each split- 
ting, and encoding momentum conservation at each vertex, which give the parton shower 
information beyond an analytic LO/LL calculation. These along with the overall choice in 
evolution variable (mass, k±, angle, etc.) are treated in different fashions by different SMC 
codes. 

There are of course further corrections to include to go to NLO in as, denoted NLO(as), 
NLL in kinematic logs, and/or NLO in power corrections to the strong ordering, denoted 
NLO (A). The most effort to date has gone to working out the NLO(as)/LL contribution to 
incorporate one-loop corrected amplitudes at the top of the shower. Adding as corrections 
involves the numerical challenge of combining real and virtual results which separately have 
IR divergences. The basic resolution is to extract the pole-portion of the real emission of 
i-partons and include it along with the virtual contributions to the i — 1 case. Unfortu- 
nately, this does not sum leading logs. One cannot blindly extend the CKKW procedure to 
NLO(as)/LL, as it leads to double-counting problems; the Sudakov factors in the reweight- 
ing contain a portion of the one-loop contributions. Separately adding the full one-loop 
result would clearly overcount. 

There are two main solutions to the NLO(q;s)/LL merging problem in the context of 
standard 1 — )■ 2 splittings. MC@NLO [37i] works by means of subtraction, finding the places 



where the Sudakovs will contribute at NLO(as), and removing the splitting function contri- 
bution. This approach is conceptually clear. Since the full amplitude and splitting function 
portions are calculated separately before subtraction, the latter for each SMC program, this 
is time-consuming. Furthermore, since the subtractions occur for the amplitude squared, one 
cannot guarantee positivity of the result and must deal with negatively weighted events. To 
avoid the computational difficulties of process-by-process subtraction and negative weights, 
an alternative is the POWHEG algorithm 38|. It keeps the IR-safe NLO(as) cross section 



manifest, and defines a Sudakov factor based on a modified splitting function to handle 
LLexp and a subset of NLLexp resummation for the hardest emission. In this way, it makes 
use of quantities already obtained in the fixed order NLO(as) calculation, requiring fewer 
additional steps for its implementation for each known process. The conservation of proba- 
bility obeyed by the splittings and related Sudakov factors avoid double countings and give 
back (Tnlo upon integration. 

A separate set of approaches goes beyond the 1 — )■ 2 formalism to consider the radiation's 
effects on one or more "spectators." The consideration of an additional parton in the pre- 
emission configuration has led to work known as dipole subtraction and dipole antennas. 



The former was initially developed in 39|, |40|. It ex plic itly subtracts the IR divergence from 
real emission via a simplified "dipole" term. Refs. 4ll |42|] have proposed algorithms based 
on these techniques. There has also been development on the theoretical side of subtractions 



by Nagy and Soper 43l446l| , with the aim of including spin and color effects, while improving 



the efficiency of implementation 47|. The original use of antennas came in the ARIADNE 
program, which treats the 2 — >■ 3 splitting as its basic unit 48l-l51| and allows for exact 
momentum conservation. There have since been more systematic attempts to extract the 2 



—7- 3 "antenna" functions from QCD and implement them in a shower, e.g. VINCIA [52 
Ref. 53| even derives spin-dependent antenna functions, though its SMC implementation is 
yet to appear. 

A different approach is the GenEvA framework 5J, |55| which allows the issues of phase 
space double counting and combining matrix elements and log resummation to be treated 
independently. This is done using effective theory ideas for how to separate scales. In this 
setup, one manifestly avoids negative weights and double counting by using multiplicative 
merging. For example, GenEvA yields a calculation that is equivalent to POWHEG for the 
NLO(as)/LL matching and at the same time a CKKW-L type LO(as)/LL matching for ad- 
ditional emissions. In a similar fashion, the power suppressed matrix element computations 
and subleading no-branching probabilities derived here could be implemented in GenEvA, 
and work in this direction is commencing. 

Another approach to go beyond LL is to incorporate the contribution of the 0{a1) cor- 



>(i) 



rections to the Altarelli-Parisi splitting kernels, Pqq . This was done to resum soft logs to 



NLL for semi-inclusive variables in DIS and Drell-Yan [56j . In order to conserve probability, 
these corrections must be correctly accounted for in both the probability for real emission in 
Eq. ([3]), as well as no-branching branching probabilities. This is related to why POWHEG 
only implements them for the hardest splitting, where they have information from the full 
fixed-order computation. The KRKMC group incorporates the subleading real emission 



contributions into fully exclusive partonic configurations in SMC 57H59|. Some of the sub- 



leading contributions take the form of 1 — > 3 splittings, requiring a modification of the usual 



1 — )■ 2 algorithm. Similar to CKKW, the KRKMC groups corrections take the form of a 
multiplicative reweighting. For a particular configuration of partons in phase space, they 
reweight by a factor that includes the insertion of 1 — )■ 3 "defects" and loop-corrected 1 — )■ 2 
splittings that account for the effects of Pqq ■ If p is the fully differential cross section, they 
define a corrected weight for n partons, Wn as: 

_ pLojkl, ■■■,kn) + ^"£^1 PN^Lo(fcl, • • • , fcn) , s 

pi^o{ki, . . . , kn) 

where r determines the number of defect insertions in any configuration. Since this reweight- 
ing involves splitting probabilities and not subleading no-branching probabilities, it does not 
clearly improve the level of log resummation. 

In this work we set up an EFT framework to classify and study perturbative as correc- 
tions, higher order log resummation and/or kinematic power corrections to parton showers. 
While the ultimate goal is to facilitate the implementation of a NLL/NLO(as) parton shower 
algorithm accounting for the leading deviations from strong ordering, our task here is much 
more modest.^ We focus primarily on kinematic power corrections in the fully differential 
cross section for an arbitrary number of final state emissions. That is, our main goal is to 
compute 

^^^3 1^ + l^^S 1^ ■ i^) 



Here NLO(A) is the next-to- leading order power correction in the cross section, which in 



we 



volves terms that are NLO(A) and NNLO(A) in the amplitude. Similarly to [2J, |25 
use an operator approach based on SCET. A main issue to resolve is taking into account 
different possibilities for the kinematic configurations of subsequent emissions, to go beyond 
the strong ordering described in Eq ([T]). The hierarchy between regions is expressed by the 
power counting parameter A <^ 1. We overcome this issue by setting up a tower of related 
soft coUinear effective theories, called SCETj, which also helps us deal with several technical 
obstacles. We formulate the shower description as a standard matching procedure between 
operators in different SCETj. Power corrections are encoded by performing matching com- 
putations at subleading order in the kinematic expansion. These corrections modify the 
processes that initiate the shower, modify certain early branching probabilities, and open 
up the 1 — )■ 3 splitting channel. Virtual perturbative a^ corrections are included by perform- 
ing matching calculations beyond tree level between SCETj theories. Finally, corrections 
to the Sudakov no-branching probabilities are encoded through anomalous dimensions of 
leading and subleading operators at the appropriate order within different SCETj 's. When 
we refer to a parton shower in the context of our calculations, we mean an explicit amplitude 
formula that would agree numerically with a corresponding shower algorithm. We will carry 
out the necessary computations for the power corrected matching equations, and a subset 
of the required calculations for anomalous dimensions occurring for operators beyond the 
LL shower. This analysis includes the leading corrections to the shower from interference 



^ In particular we note that soft NLL resummation may only be feasible at leading orders in 1/Nc 60l.l61|. 



and from spin correlations. As much as possible, we attempt to give pointers for additional 
computations that are needed in places where our analysis is incomplete. For example, to 
simplify things we have not treated color correlations since doing so increases the basis of 
operators and the number of computations, but does not change the conceptual setup. 

The outline of our paper is as follows. We present a brief overview of SCET in Section Hi A[ 
We review the Bauer-Schwartz SCET shower method in Section lll Bl and discuss the technical 
obstructions to extending it to include power corrections. In Section IIIC[ we present our 
SCETj framework to resolve these issues. In Section llllt we analyze the LL shower in the 
SCETj framework, and show that the transition between SCETs, SCETj — )■ SCETj+i, can 
be encoded by operator replacement rules on single parton coUinear fields. Soft emissions in 
SCETj are discussed, and we summarize the correspondence between SCETj objects and LL 
shower ingredients. In Section llVt we use the SCETj formulation to classify and compute 
various corrections to the shower to (9(A^) in the cross section. Two main categories of 
branching corrections emerge, which we refer to as "hard-scattering" and "jet-structure." 
We also discuss ingredients needed for renormalization group evolution corresponding to no- 
branching probabilities, derive all the LL anomalous dimensions for our subleading operators. 
Additionally, we mention the issues involved in obtaining NLLgxp resummation from our 
results. A summary of corrections in the SCETj framework is presented as a table in 
section lIVEt including the type of corresponding ingredients needed in a subleading shower. 
We present in Eqs. fll04p -f lT06|) a parton shower reweighting factor that should allow one 
to implement our corrections. We also discuss the correspondence of these corrections with 
those currently included in other Monte Carlos. Conclusions are given in Section |Vl At the 
present time, we do not have an algorithmic implementation of our power suppressed shower 
results, but work in this direction is in progress. 

Many details are relegated to the Appendices. Further details about SCET can be found 
in Appendix |Al We describe finite reparametrization transformations in Appendix |Bl which 
is an important symmetry that we use in our matching computations to disentangle kine- 
matic coordinate conventions from kinematic power corrections. Details on the matching 
of QCD -^ SCETi, SCETi -^ SCET2, and SCET2 ^ SCET3 can be found in Appendices 
ICl [Dl and El respectively. A complete list of the operators needed to compute Eq. (jH]) in 
SCET AT is given in App. [El Appendix [Fl contains a cross-check on our results, where we 
integrate a subset of our power suppressed terms to rederive the abelian terms in PqUqg, 



namely the 0{as) correction to the q ^ qg splitting function [62 . 

Those readers looking to find a quick summary of our results should look in Sees. IIIICI 
andHVB 

II. OBTAINING THE PARTON SHOWER WITH SCET 
A. SCET Basics 

Soft-Collinear Effective Theory is an effective field theory of QCD that describes the 



interactions of collinear and soft particles 20l-l23|. We present here the basic ideas needed 



for our analysis of the parton shower, including how collinear sectors are organized into 
equivalence classes by the power counting parameters. Further SCET concepts are reviewed 

8 



in Appendix \M 

The momentum, p, of any particle can be decomposed along two light-cone vectors, n 
and n, with n^ = 0, n^ = and n ■ n = 2, as 

p^ = p—+p^^ + n-p—, (9) 

where p = n-p and the particle's invariant mass is p^ = n ■ p p + p\. We use a Minkowskian 
notation for p]_ = —p]_, where p± is Euclidean. SCET's degrees of freedom include n^- 
collinear fields for a set of distinct directions {ui}, and soft fields.^ A particle is coUinear to 
a direction n if its momentum scales as: 

{n-p,p,p^)r^{X^l,\)p, (10) 

where p ~ Q is some hard scale in the process, and A ^ 1 is the SCET power counting 
parameter. A particle is soft if 

(n■p,p,J9x)~(A^A^A2)g. (11) 

Collinear and soft fields have virtuality ~ Q^X^ and Q^A^, respectively. We obtain SCET 
from QCD by expanding in powers of A, integrating out hard modes, and dividing the 
remaining ones into collinear and soft fields. Our collinear and soft degrees of freedom also 
contain all the IR regions that can be obtained by a rescaling of A — )■ A', for i > 1. The 
leading order SCET Lagrangian is 

riO) _ riO) , Y^ r{0) /.px 

-'-SCET — '-s ^ 2-^ >~'n 5 \^^) 

ne{ni} 

where Dn is defined in Eq. (lA9p and has only interactions among particles collinear to the 
same n. Cs is the Lagrangian for soft interactions discussed further in App. |X1 Particles 
collinear to different directions can interact either by the exchange of soft modes, or from 
their coupling to other sectors in external operators. Two collinear sectors in SCET, n\ and 



n2, are distinct if 63 



r^l■n2>A^ (13) 

so any particle is collinear to at most one direction within a given SCET. The collinear 
sectors \n^ in SCET are really sets of equivalence classes of null vectors, {[^j]}, where 
the equivalence class is [rij] = \n ^ ["-j]|^ ' "^i ^ -^^l- A class [n^] consists of all light- 
like vectors connected to n^ by a type-I reparametrization invariance (RPI) transformation, 
n^, — 7- r2^ + A^._|_, where the scaling of the transformation parameter is A^._|_ ~ A (see App. IB] 
for a detailed discussion of RPI). Physically, the class [rij] corresponds to light-like vectors 
for particles whose momenta is in a cone centered on n^ with an opening angle ~ A (c/. 
Fig. [15]). 



^ Our primary interest here is the perturbative structure of jets, so we use SCETj theories with colhnear 
and ultrasoft modes. For simpUcity we wiU always use the phrase soft in place of ultrasoft. 



Thus, the defining concepts of a SCET-theory are its hard-scale Q, its colhnear sectors 
{[rij]}, and its power counting parameter A which governs the importance of operators and 
the size of the colhnear sectors in phase space. 

Most of our discussion will involve interactions with colhnear fields, and we use the nota- 
tion Xn for quarks and B':^^ for gluons (definitions of these fields can be found in Eq. (lAlSp . 
and they incorporate colhnear Wilson lines built out of n ■ An fields). We can match QCD 
onto a series of SCET operators organized by powers of A. The key building blocks are: 
Xn, -B^x, and V^^ (a type of derivative operator that yields the perpendicular momentum 
of an n-collinear field), each of which scale as A in the kinematic power counting. A general 
notation for the z-parton operators we will consider is: 



i/2 

^b=j /2+1 



'-a=l -I ^h=j/2+l -' '-c=l 

(14) 
where the number of partons is the sum of quarks and gluons, j + k = i, and the total 
number of _L derivatives is ^ = J2m=i ^rn- In the operator argument, we list the index labels, 
Ud, of the parton fields on the RHS. The superscripts in the argument on the LHS denote the 
number of derivatives acting on the field with the corresponding direction. There may be a 
degeneracy among the index labels, rid, and so the operator has at most i distinct colhnear 
directions. The scaling of these operators is O'^^'^'^^ ~ A-'"'"'^"'"^. They are tensors in the space 
of spinors and Lorentz vectors, and the indices get contracted with structures contained in 
the Wilson coefficient C for the operator. If CO is a Lorentz scalar, then j is even. Since the 
colhnear fields carry a label referring to a specific light-cone vector, these operators describe 
particles in a specific region of phase space. SCET therefore distinguishes situations with 
the same particle content, but different kinematics, in a straightforward way. 

For example, one can take an amplitude for three external particles: a quark, gluon, and 
antiquark. We can consider two different configurations, lo'nofi'no^n) and Iqmgn'Qn)- In the 
first, shown in Fig. W^I), the quark and the gluon are no-collinear, and the antiquark is 
colhnear to a different direction, n. Here the amplitude is described by operators with two 
distinct directions, say 

O(2A0)(^^^^) ^ ^^^Yxn ~ A^ , O^''''''\no,no,n) = Xno9B'n,^T'xn ~ A^ , (15) 

where the form of the Dirac structures F and F' are not central to our discussion here. (9'^^'°'°) 
can emit n ■ An^ gluons from the Wilson line in x„q, but requires a Lagrangian insertion to 
emit an A:^^ gluon. Schematically, the amplitude for a transverse gluon has contributions: 

A^ = dx (0|r {C^noi^) Xno^XniO)}\qnognoqn) + {0\XnogKo±'^' XniO)\qnQgnoqn) ■ (16) 

In Fig. ^11), each of the particles is colhnear to a distinct direction, so no cone of size ~ A 
fits two of the momenta. In this case, the amplitude can only come from an operator with 
three distinct labels, such as XniB^, jT"xn'- 

A" = mn,9B';^,^J"xn\qn,9n',qn) ■ (17) 
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FIG. 1: Different kinematic configurations of a final state witli a quark, antiquark, and gluon are 
described by different SCET operators. In (7), the quark and tlie gluon are collinear to the direction 
no, represented by their sharing a common cone. In (77), the vectors q'l and A;^ are too far apart 
to be collinear. The Feynman diagrams show that collinear particles can come from Lagrangian 
insertions, whereas non-collinear ones arise exclusively from higher-multiplicity operators. The 
Feynman diagram in (7) only depicts the first term on the RHS of Eq. (J16p . 



The ability of SCET to cleanly separate contributions such as those in Eqs. flT6l) and ([Tj 
will be useful for formulating a complete set of power suppressed corrections to the parton 
shower. 



B. Bauer-Schwartz Method 



The original application of SCET to study and improve the parton shower was carried 
out in 23, l25| by Bauer & Schwartz. The main reasons why SCET is useful for this are: 



• The SCET fields, soft and collinear quarks and gluons, have support in the infrared 
exactly where the parton shower amplitudes have their dominant contributions in 
phase space. 

• Since SCET is improvable order-by-order in the kinematic expansion parameter, A, 
one has the potential to systematically correct the shower. 

We will give a short overview of the Bauer-Schwartz approach, and then discuss the com- 
plications that arise when trying to extend the analysis to NLO in the A expansion, namely 
NLO(A). In this section we will use notation that is not found elsewhere in the paper to 
retain consistency with 2J, |25 . 
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The procedure of [2J, |25| starts by constructing i-parton operators, Oi, through matching 
SCET to QCD at a hard scale. For example, their O2 will equal O^"^'*^'^' (rii, ^2) in the 
notation of Eq. (HM . and C3 will be 0^'^'^'^\ni, n2, ns). As we run Ci(/i) down, the leading 
log renormalization group evolution (LL RGE) does not mix operators and the exponential 
evolution kernel encodes the no-branching probability. The evolution continues until another 
parton becomes apparent at a scale fi = pt- 

If we have an z-parton operator, Oi = 0^^''^~^''^\ni, . . . ,ni) with all n's distinct, then it 
has the RG solution Oi{fi) = U'-^^'-^^'^\Q, fi)Oi{Q) with 



U^'^-^'^^XQ, fi) 



exp 






(18) 



where ^(•''^"J'O) is the operator's anomalous dimension. The leading-log resummation effects 
of the Sudakov factor in the PS enter through one-loop operator running in SCET, as 
dictated by the cusp anomalous dimension. The one-loop cusp portion is especially easy to 
calculate in SCET as it depends solely on the number of coUinear fields, even though the 
calculations have loops involving soft ones as well 2M, I25| . 



(n,,7ig,0) 

7ll 



(/^) 



TT 






1 -" 



(19) 



This form of the kernel gives a product of Sudakov factors which are the no-branching 
probabilities for each parton in the operator: 



U, 



LL 



iQ,fx) = A|(Q, /i)A/ (Q, ^). 



(20) 



Here, as in 27|, one accounts for leading-log effects for any particle multiplicity by simply 
multiplying matrix elements by appropriate Sudakov factors. 

As we run Oi{fi) down, another parton becomes apparent at a scale /i = pt- To account 
for this, Bauer-Schwartz devised a "threshold matching" of Oi to a new, higher multiplicity 
operator, Cj-+i, where the subscript still denotes the number of partons in the operator and 
the superscript tracks the parent operator. The general threshold matching equation is 



H'\oi 



J) 



A«=PT+<: 



[ciUo.. 



n+1 



fi=prp-e 



(21) 



After further running and threshold matching, we eventually have On for various n > i. 
The n — i particles emitted at increasingly lower scales by this process correspond to the 
parton showering of the original fields created at the hard scale by Oi. Additionally, they also 
showed that an appropriate list of SCET operators (0j's and O^ s) can interpolate between 
fixed-order QCD and parton shower (PS) calculations of IR-safe observables. Furthermore, 
they derived the 0{as) effects from matching QCD to SCET at one-loop for O3. 

That subsequent emissions reproduce the usual parton shower splitting function emerges 
easily from SCET. Consider an operator Oi = Xno^; where Q is arbitrary and we have 
made explicit a single coUinear quark field, Xno ■ If we emit a coUinear gluon from this quark, 
q{Qo) ~^ lil'Ddiki), the amplitude for the process is 






(22) 
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where Ung is the colhnear quark spinor, and p" is the combination of the SCET single gluon 
emission Feynman rule plus the Xno Wilson line emission (the quark Cn can be found in 
Eq.lXQJ), 



„ (^l)«0±7no± 7no±(^o)no± n° 

qi qi qo 






Jno- 



Mo- 



qi 



(23) 



Note that p" in SCET comes entirely from Xm without reference to anything residing in Q. 
The subscript (tiq-L) refers to components perpendicular to n^ and n'^, which we denote by 



_L for the remainder of this computation. The amplitude in Eq. fl2^ is gauge invariant and 
kfpa = 0. Squaring A^^'^^ and summing over spins we have Y.spin'^no{qi)uno{qi) = qiiof^, 
and the gluon polarization sum denoted J2spin^a^*p ~ ^ai3- Since p"" commutes with ]?4q, 
we get an answer proportional to p^'p'^^ c/q,^, where without loss of generality we can use a 
light-cone gauge, dap = —Qap + (^a^i^ + ^ia^/3)/^i- Crucially, this is a Dirac scalar: 



p'^p^'^dap = IpI 



2gi 
kiqo 



qj± , 2gox-gi± g^^\ 



qoqi 



qo 



(24) 



where we have used the on-shell conditions qf = and kf = 0. 

In a frame where g^ = we have qi± = —ki± and qo/qQ = 1/(^0 ■ qo). Here no ■ go = 
uq ■ ki + riQ ■ qi = —kfj_/[qo z{l — z)], where z = qi/qo- Thus we have the simpler expression 



P^ 



< + 



i^^)no±lno± 



qi 



(25) 



which we have written in light-cone gauge without the Wilson line contribution (oc n'^qQ), 
and 



7? / "'' On 2:^(1 

Putting these properties together in the full amplitude squared we get 



(26) 



\A^^''\' 



9^Cf qi 
(no • go)^ 2 

/C'i.2^^ii4^Tr 



Tr [i^^p'^nn^ p^^] d 



■a/3 



^^Cf gi 
(«o ■ ?o)^ 



IppTr 



i'o 



nn^ 



\k 



iJ-l 



2 



go^fifi^ 



(27) 



Thus, all information about the emission factors out to the front and is independent of the 
rest of the process encoded by Q. Since the power expansion is built into SCET, there 
was no need to expand terms in the amplitude to obtain this result (unlike the analogous 
computation in full QCD). In order to recover Eq. ([3]), we still need to include the z- 
dependence from phase space, since P-^ (z) operates at the level of the cross section. Using 
d^k/{2Ek) = dkd'^kj_/{2k), for gi and ki we have 



dqid'^qi±dkid'^ki± dqod'^qo^ dzd'^ki± 



2gi 



2ki 



-)■ 



2go 2z{l 



(28) 
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where the arrow means that we insert d^qoS^^\qQ — qi — ki) and integrate d?qi along with 
dijiQ-qo). Thus, we recover the expected 1/(1 —2) dependence from the measure. Combining 
pieces and performing the trivial azimuthal integral (i0fci, we get the expected expression: 

dk'^ 
dax+gg = dz -l^P^\{z) dax+g, (29) 

where Pq4qg{z) is the quark splitting function in Eq. ([3]). Here dax+q is the cross section 
for the rest of the process with emission of a momentum go quark, and the corresponding 
amplitude squared is Tr [^^^nfi"'"] . Whether fl represents a simple hard current or an entire 
chain of coUinear splittings, we see that the q ^ qg emission factors out with the expected 
soft-coUinear double pole, as in Eq. ([2]). 

In order to obtain their results, Bauer-Schwartz introduced choices and approximations at 
several points which obscure the path toward systematically computing NLO(A) corrections. 



Indeed, they concluded that obtaining these corrections may be prohibitively difficult [25 
Some of the issues one encounters trying to work at higher orders are: 

1. At NLO(A), it becomes crucial to distinguish which simplifications correspond to ap- 
proximations with power corrections, and which involve a choice of coordinates where 
a symmetry makes the final answer coordinate independent. For example, a collinear 
state typically has nonzero momentum components perpendicular to the index n of the 



field that annihilates it. Refs. |2J,|25[, however, dictated that collinear SCET fields in 
their operators only create particles whose momenta perfectly align with their index 
direction, n: 

Xn\q) = Sn,n,, whcrc n^ = q^/Eq, (30) 

leaving it ambiguous what amount of symmetry protects this choice. Eq. ( !30l) enforces 
certain kinematical restrictions on final state particles, and requires that fermion fields 
be rotated to an appropriate n^ via ^„ — )■ {TijiT/i/A)^nq- 

2. At LO, it was possible to avoid a potential double counting between collinear and soft 

fields by dropping soft emission and Wilson line emission, and taking only collinear 

emissions with transverse polarization. The threshold matching procedure is designed 

to avoid double counting of collinear operators, such as a Lagrangian emission from 

(2) 
O2 and direct emission from 0\ , since only one of these is allowed to operate at a 



time. However, the threshold matching in Eq. fl2T|) makes the technical procedure for 
incorporating power corrections unclear. 

3. Threshold matching contains another impediment to systematic improvement. 

Through this procedure, the initial operator O2 has nonzero projection onto Fock 

states of any multiplicity, but the number of particles created by an operator is a 

scale-dependent question. The matching scales are determined by the strong ordering 

kinematics, pi^ ^ . . . ^ VmL- At the scale of an emission, say pi^, one threshold 

(2) 
matches to the operator 0\ , which only adds one parton at a time. However, going 

to higher orders in the shower necessitates more general configurations. 
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In carrying out their method, Bauer-Schwartz carefully enumerated the above approxima- 
tions. They affect the ability to include corrections in A, but do not impact the terms 
necessary for a LL shower. 



Building on the work of Refs. [2J, |25[ , the main goal of the framework we develop in the 
next section is to overcome this list of issues so that we can determine power corrections to 
the shower using SCET. 

C. Using SCET, 

The main feature of the parton shower is the ability to capture the dominant physics 
of particles emitted in kinematically hierarchical regions of phase space. Our goal is to 
formulate the SCET interface with the shower using a standard sequence of matching and 
running steps in different versions of SCET, 

QCD -^ SCETi -^ SCET2 -^ > SCET^ . (31) 

We refer to this as the SCETj procedure. The key distinction between a SCET at one 
stage and the next is the definition of the corresponding resolution parameters 1 ^ Ai ^ 
A2 ^ ■ ■ ■ ^ Xn, where Aj sets the power counting for SCETj. As we move down the chain, 
the corresponding SCET resolves smaller ~ [QXj)"^ invariant masses and relative squared 
perpendicular momenta, and has a different meaning for its collinear sectors {[^i]}scET • To 
keep track of this, we will attach a subscript to the operators to denote the SCETj in which 
its fields live, 

OP'^)(ni,...,n,+,). (32) 

Effectively with Eq. (13T|) . we partition the momenta of partons in the shower history into 
classes, 

Vto D ni D . . . D VtN , (33) 

where Qj contains the momenta of all propagators having p"^ ~ [QXjY or smaller, or an 
equivalent condition on relative perpendicular momenta. The allowed momenta in fij cor- 
respond to the collinear modes of SCETj. The sequence of SCETj's is truncated when we 
resolve a scale of order the parton shower cutoff, QX^ = p™* ~ 1 GeV, that is in SCET at. 

Note that we do not associate a large hierarchy to the hard scales pi between SCETj 
and SCETj+i. That is to say we do not associate the energy loss due to splitting with 
a power of Xj. Instead if Q is the scale of the primary hard interaction then we consider 
pi ~ V'' Q i^ SCETj, where i] ^ Aj and for numerical estimates we can take r/ ~ |. (For 
each branching the geometric mean of the two daughters' p fraction averages to 0.4 which is 
roughly one half.) Parametrically, the decrease in the parton energy is not as rapid as that 
for the perpendicular momenta encoded in the power counting parameter Aj. In principle, 
we can account for r/ as a separate factor. In practice, we will be most interested in tracking 
powers of Aj and will only include t] factors in places where the corresponding powers of two 
have a numerical impact on the implementation, or if we wish to disentangle the changes in 
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offshellness due to strong-ordering effects and those coming from the more modest decrease 
in pi. 

The strongly ordered configuration of partons in Eq. ([T]) corresponds to removing a single 
g| in flj as we pass from flj -^ fij+i- However, with Eq. f l33|) . nothing stops us from having 
multiple emissions at a single scale. If two mother particles, with g| and q'|+i, are associated 
to the same Qk, then when we integrate out that scale in SCET^+i this configuration just 
contributes to an operator with a different parton multiplicity from the strongly ordered 
one. Thus, with Eq. ( 13T|) there is no obstacle to considering corrections from an arbitrary 
assignment of g^'s to QkS. This resolves issue 3. of Sec. IIIBI since we can treat emissions 
where the shower tree has momenta with the same parametric scaling in A. 

To carry out calculations in the SCETj framework, it is convenient and sufficient to take 
a specific definition of the power counting parameters, Aj = (A)*. We want the hierarchy 
between neighboring splittings to stay the same throughout the shower so as not to privilege 
any portion of it. We will see in Sec. lIVDI that this democratic setup allows us to interpret 
part of our 0{X) corrections to z-parton amplitudes as universal corrections to the splitting 
probability, given at LO by Eq. ([3]). As we go to lower scales, our definition of coUinearity 
also changes, and by analogy to Eq. (fTOj) . fields coUinear to n within Qi have: 

(n ■ qi, qi, q^^) ~ (A^*, 1, A^ g^, (34) 

and virtuality ~ (gj)^A^*. In SCETj, £„ again only couples collinear fields in the same 
direction n. Since different SCETj 's have different definitions of coUinearity, our description 
of identical physical processes changes when we switch to a theory with a lower scale. For 
convenience, we will use the same auxiliary vector fi^ for any ?7,j-collinear field in any SCETj. 
If 77, is a valid auxiliary vector for n-collinear fields in SCETi, then it is readily apparent 
that it will be a valid choice for all subsequent collinear fields in SCETj's that descend from 
an ra-coUinear mother in SCETi. Our default choice is stronger: given a set of light-like 
vectors in {rij} in SCETi we take a light-like n that is parametrically close or aligned with 
the antiquark direction. We then adjust the magnitude of n° and of n^ so that 't-j = and 
n ■ rij = 2 (for a related discussion based on RPI see Appendix ICl). 

We depict the different descriptions of the same physical configuration in Fig. [2], where 
the left panel is in SCETj and the right panel is in SCETj+i. In SCETj, the quark (gl) and 
gluon (/ci) are no-collinear. This means that at LO they are emitted from a qqg vertex in 
the LO SCETj Lagrangian (or a Wilson line interaction). Schematically, the amplitude for 
a ±-polarized gluon looks like'^ 

^,,, = C(2'°'°) ldx{0\T{CscETA^)O^'''''^}\qqg), (35) 

namely like the first term in Eq. fITB]) . The right-hand panel of Fig. [2] denotes the same 
configuration as seen by SCETj+i. The scale of this theory is lower and the definition of 
coUinearity stricter, so the quark and gluon are not collinear here. Therefore, the amplitude 
now comes from a three-parton operator, 

^qqg ^ C'(2'l'0)(0|O(2 = 1.0)|gg^), (36) 



From here on, we will drop the superscript (0) and the subscript n from the collinear Lagrangian. 
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FIG. 2: The same three-parton process as seen in two different SCETs, SCETj and SCETj+i. 
Above: Kinematic configuration of the quarks and gluon. The solid cones represent the regions 
considered cohinear to the vectors drawn. Below: Feynman diagrams for the corresponding ampli- 
tude. Note that in SCETj+i we have removed a degree of freedom that propagates in SCETj. The 
amplitude thus comes from a higher dimension operator 0\J^, rather than from a time-ordered 
product of >CscETj with 0\ , as it did in SCETj. 



as in Eq. (^. We match SCET, -^ SCETj+i to calculate C^^.i-O. 

Given the above conventions and with the notation in Fig. [2] at hand, it is worth stating 
some simple kinematic relations that we will use later on. Take an no-collinear mother 
particle of momentum q^ = q^in^ (2) + go('T'o/2). Let go decay to two onshell massless 
daughters, ki and gi, with momentum fractions x and (1 — x), back-to-back ±- momenta k±, 
and light-like directions rii and n\, then 



K 



n'' 



Un 



n 



/fj. 



kt\ + h'^ + k^. = h'^ 



?r 



n'' 



rin 



W, 






(37) 



Note that our convention of using the same n^ auxiliary vector ensures that in these decom- 
positions the momentum multiplying Uq is the same as the momentum multiplying n^ 
The collinearity of ki and gi can be determined by the size of k\, Qq, or rii ■ n[, and the 



(')m 



relation between these three choices is 



2fc2 



rii ■ n^ 



2%^ 



(go)2x2( 



X] 



(go)2x( 



X] 



(38) 



Since we take kj_/{qoY ~ Af in SCETj, we have ql/iqoY ~ ^"i/v"^ ^^^ '^i '^'i ~ ^1/v'^- Thus, 
all three choices are equivalent for counting powers of Aj, but differ with respect to how 
powers of the energy loss parameter r/ ~ 1/2 appears. 

After this introduction to SCETj, we now list some technical advantages of this framework 
for our analysis: 
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1. Collinear fields in SCET with different n-labels, as well as soft fields, do not over- 
lap in Hilbert space. This allows us to separate an z-jet process with i distinguished 
partons, from an (z — l)-jet process with i partons, where two are collinear and un- 
resolved. Lower-scale SCETj's distinguish configurations more finely based on their 
stricter definition of collinearity. This resolves issue 2, avoiding the double-counting of 
similar configurations, from Sec. IIIBI This SCET property also illuminates simplified 
structures in the power corrections, such as the form of the amplitude interference 
(cf. section IVD]). 



2. Soft modes communicate between collinear sectors and threaten the factorization of 
different jets. Fortunately, SCET constrains the interactions they have with collinear 
fields. In fact, one can decouple them using soft Wilson lines in the LO SCET La- 
grangian. At LO, using the SCETj soft Wilson lines, we maintain factorization, obtain 
angular ordering, and rederive the coherent branching of soft emissions (cf. section 
IIIIBI) . Soft interactions which are power suppressed can also be systematically stud- 



ied in SCET with Lagrangians available in the literature 64J-|66| . which we give in 
Eq. ([lOHD. 

3. In SCETj, we have a symmetry group RPIj which corresponds to coordinate choices. 
In SCETj+i, only a subset of this, RPIj+i C RPIj, remains a symmetry of the new 
theory. The kinematics in the coset portion RPIj/RPIj+i within SCETj give a set of 
higher- dimension operators in SCETj+i, and describe configurations which would not 
otherwise be contained in the SCETj+i Lagrangian (cf. section UTTl and Appendix iBj) . 
This resolves issue 1. from Sec. IIIBI making the difference clear between approximations 
and conventions chosen for simplicity. 

4. In matching between SCETj and SCETj+i, suppressed operators in the lower-scale 
theory are needed to reproduce the physics of the higher one. It can be proven that all 
higher order purely collinear operators can be built from quark fields (Xn), perpendic- 
ular gluon fields {B±n), and the perpendicular momentum operators {V±n) |67[. Thus 
the symmetries and equations of motion of SCET greatly simplify the operator basis 
one needs to consider at each order in A (cf. section HVl and Appendices ICl iDl andJEl). 

The final SCET^v corresponds to the scale where the shower stops, i.e. where Qrj'^X'^ ~ 
p^*. In SCETtv, we only need the coefficients of the operators where all collinear partons 
have distinct ?7,-labels, and which have no "Pnx's, C}^' ' O^' ' . Once we reach the physical 
resolution scale, it is only meaningful to have one collinear parton in each distinguished block 
of phase space. Using RPIat, we can set n^ = p'j/p^- This is as in Eq. fl30|) . but we only do 
this when we run up against the physical limit that requires just one parton per equivalence 
class. At intermediate stages, we allow different fields to share n-labels, which also results 
in operators containing Vn±- The coefficients C}^' ' encode the history of the shower. They 
can be written entirely in terms of: dot products n, ■ nj/, equivalent to products of final 
parton momenta, which carry the scaling in A; hard momenta pi, the renormalization scale 
/i, and collinear cutoff parameters encoded in B-functions. 

As far as the shower is concerned, A is merely a bookkeeping device which determines what 
pieces are needed beyond LO. One could try defining Ai = ki±/Q, A2 = k2±/Q, etc., but this 
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FIG. 3: Operators that reproduce strongly-ordered gluons are constructed through a series of 
matching computations with emissions in different SCETj. The horizontal dashed arrows refer to 
the radiation of a gluon from a time-ordered product of the SCETj Lagrangian with the operator 
creating fields at the point marked by (8). The diagonal solid arrows denote the matching onto a 
higher multiplicity operator in SCETj+i. 



is not ideal since there is a chance for events where ki± ~ Q or ki± ~ p§^*. The organization 



in Eq. (1551) instead exploits the fact that on average showers are strongly-ordered. Our 
expansion in A will then on average give a description of the most likely deviations from 
strong-ordering. Our goal in using the SCETj framework is to extract an amplitude suitable 
for reweighing the parton shower to this level of accuracy.^ From the SCET side, we pass to 
the shower weights built from SCET^v squared amplitudes (c/. Eqs. 1104111071) . They contain 
the information needed to describe a strongly-ordered shower and its leading kinematic 
corrections. 

Before proceeding to our computations, it is worth commenting explicitly on which shower 
ingredients we do not compute. We only treat the case of a showering quark q ^ qg and in 
general take the abelian limit of QCD {Ca = 0). We have left out gluon splittings, g ^ qq 
and g ^ gg, from this analysis, though we expect that the extension to these cases should 
be straightforward. We have also not determined the effect of NLO(A) power corrections 
from subleading soft interactions, although we briefly examine the factorized structure of 
LO softs in section flIIIB|) . These items are all left to future investigations. 



^ As a well-defined EFT, one certainly could also do standard factorized cross section computations in any 
SCET, if one wanted. 
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FIG. 4: The opening angle of the hght grey (blue) cone is ~ A^*, and the opening angle of the dark 
grey (red) one is ~ A^*-*^^'. The particle with momentum p is collinear to both n and n' in SCETj, 
but only to n' in SCETj+i. RPIj allows us to move the field label, n, to any location inside the 
appropriate cone for SCETj while keeping the theory invariant. 

III. PARTON SHOWER IN SCET VIA OPERATOR REPLACEMENT 

In the previous section, we presented our approach of using a series of EFTs, the SCETj, 
to handle processes with a hierarchy of many scales. We will now use this technique to 
calculate the leading contribution to a series of collinear emissions, as occurs in the parton 
shower. Our ultimate goal is to incorporate corrections, but as a starting point we want to 
easily reproduce the strongly-ordered configuration of Eq. ([1]). We can do this if we declare 
that in a shower, the z*'^ particle decomposes as: 

(n ■ qi, qi, q^) ~ (A^*, 1, A^ g^, (39) 

and therefore has virtuality qf ~ {qiY^^^ {cf. Fig. |8]). This is exactly the same condition as 
Eq. dMD, which we used to define the EFT, SCETj. 

To calculate the operators that describe i emissions in the strongly-ordered limit, we will 
perform a series of matchings SCETj — )■ SCETj+i. We will find that the most efficient way 
to describe the process at LO in A is to be in SCETj+i for i-parton radiation. Thus, we 
emit and match i-times in series, as shown by Fig. [3l At LO, we will show that one can 
implement this using an operator replacement rule. In the case oi q ^ qg emission, it takes 
the form: 

Xni^cgl3Z,±Xn2, (40) 

where Xn and B"j_ are the SCET fields associated with collinear quarks and gluons, respec- 
tively, and c is the Wilson coefficient whose spin and color indices are suppressed. Though 
we do not compute them, there are similar B"^j_ — )■ c'xn2Xn3 + '^" ^n2±^n3± rules as well. 
In SCET, each collinear field carries the label n, which gives its direction of collinearity. 
Note that the quark field on the LHS of ( HOl) has a different one from those on the RHS. 
This relates to the stricter definition of collinearity in SCETj+i shown in Fig. HJ In order to 
perform the matching, we will make use of the reparametrization invariance (RPI) discussed 
in point 3. of Sec. Ill CI to change fields' n-labels. 
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A. Leading Sho^ver Revisited 

We first want to reproduce the strongly-ordered contribution to z-gluon radiation from 
tlie quark in an initial 7* — )■ qq pair production. Our iterative matching procedure for 
multiple EFTs takes a particularly simple form at LO in A. For our standard example, we 
take the process e"'"e~ — )■ jets. Starting in QCD, we couple the quarks to another sector 
via the operator, Jqq-q = qV^q. This allows us to avoid complications that come from the 
initial state such as backward evolution. In SCETi (which is equivalent to the usual SCET), 
matching to QCD at tree-level converts the quark coupling to the following operator at LO: 
Xno^^Xn, which produces q and q in different collinear directions. Details on the matching 
of QCD to SCETi are given in App. O Using the notation in Eq. (in|) . we write the SCETi 
operator in the following way: 

where 

(of'°'°'(n„,fi))_^ = (x„.),(x-i,)j, (42) 






and i and j are spinor indices. The subscripts 1 in Eq. fH2|) indicate that the fields are 
defined in SCETi. Our focus is on gluon emissions from the quark, and we always take the 
antiquark in the same direction, n, therefore we drop it from the list of n-labels. Also, we 
will use the following shorthand notation for the most common operator, 

of '^'°Vi, n[, . . . , nl n) ^ Of\n,, n[, . . . , n',) , (43) 

where the subscript marks these as being in SCETj. In the rest of the paper, we will often 
drop the spinor indices. Using the above convention, we write the operator in Eq. ( 14T1) as: 

Xno^^Xn = C^^^lo0^f\no) . (44) 

The LO derivations are independent of the exact structure of F'^. In fact, even the antiquark 
is a spectator, and we could just as easily use O^''^ = Xno^; where Q is arbitrary. However, 
as we will discuss in Sec. HVl matching QCD to SCETi at higher orders requires us to specify 

n. 

To calculate operators in SCET2, we start with single gluon radiation. In this case, shown 
in Fig. [5l the emission amplitude is:^ 

A^S = C[%o{0\ldxT{C^''''^^ix)Of>\no)}\qnognoqn) (45) 

= gunM U + ^^^^"-^^"°^ ) %r^v,ip,), (46) 



All the amplitudes we write in this work refer only to the hadronic part of e^e — ?■ jets, thus A'^q is the 
amplitude of 7* ^ qqg. 
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FIG. 5: Momentum labels for single {A) and double (B) gluon emission. 



where we have labeled the collinear directions of the particles in the state IgnoOmQu) for later 
convenience. The SCETi Lagrangian is given in Eq. flA9|) . Here we study the process in the 
center of mass frame with p^ = {Q, 0, 0, 0) and the quark (go) and antiquark (pg) along the 
directions uq = (1, 0, 0, 1) and n = (1, 0, 0, —1), respectively: 
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We decompose the emitted quark (gi) and gluon {ki) along the directions (no,n) 
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2 u ' \ J-/no-L ' 2 

The variables are illustrated in Fig. |5l By momentum conservation we have (fci)no± = 
~('?i)no±) Q = ^0 = ^1 + ^1 and riQ ■ pq = Q — tiq ■ qi — riQ ■ ki. We take all the external 
particles on-shell, thus uq ■ qi = —{qi)nQ±/Qi and similarly for uq ■ ki. As we discussed in 
Section IIIBt 24j . |25| showed that single gluon emission in SCET reproduces the splitting 
function, Eq. ([3]), and factorization behavior, Eq. (127|) . of the standard parton shower. This 
simple behavior for a single radiation will reproduce the shower for an arbitrary number of 
gluons. 

We now want to match the single emission to SCET2 (c/. Fig. |6]). There is a slight 
technical complication due to the different definitions of collinearity in the two theories, as 
illustrated by Fig. |H In SCETj, a collinear field with label n can annihilate a state containing 
a particle whose momentum vector lies anywhere in a cone with angle ~ A* about n. When 
we change to a lower-scale theory in a matching equation, we have to take care that the 
operators' n-labels are appropriate for the desired amplitude. Using the terminology of 
Fig. m while any label vector in the light grey (blue) cone is sufficient for a particle with 
momentum p in SCET,, for SCETj+i we need one in the dark grey (red) cone. This is where 
RPIj comes in, as mentioned in Sec. IlICi We use it in SCET, to transform all quantities 
in the amplitude (spinors and vectors) that depend on the label vectors, such that the label 
after rotation lies within a collinear cone with angle ~ A*^^ about the particle momentum. 



22 



@ 



no 




ni 



n\ 



FIG. 6: (Left panel) Single gluon emission in SCETi conies from the time-ordered product of 
the Lagrangian with a quark-creating operator, A = (0|T{£gcETiC'i }\QgX). (Right panel) For 
parent quarks with virtuality ^ Q^A , the gluon comes from the central vertex in SCET2 via a 
higher-dimensional operator, A = (OIO2 iQd^)- 



The simplest convention is to choose the n-Iabel to align perfectly with the particle. If 
desired, we could make any choice consistent with RPIj+i transformations. For the process 
under consideration, we define labels, rii, n[ such that. 






(49) 



In SCETi, we are free to use uq or tt-i to describe the gi quark and ki gluon because of the 
RPIi symmetry. Since ni is a valid index for the quark field in SCET2, we do the matching 
computation using the same spinor, M„^(gi), in both theories. In App. [Bl we derive the RPI 
transformations we use here and other rotation formulas. For now, we quote the results we 
need: 
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As required, the two different n^- vectors' directions lie within cones of size A about hq. It is 
simple to check that in the new basis, {qi)ni± = <?! ~ (^1 ■p)n/2 — qi ni/2 = and similarly 
for {ki)n'±. Acting on Eq. ( H6|) . we get: 



A qqg "0 — 
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OpM 



r^t;. 



(51) 



where go = ^i + ^i- Having changed bases, we can easily write the SCET2 operator that 
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reproduces Eq. f l5T]) . CgLQOj (^i)"-i)) where:^ 

<^2,LoK'^i)=f^LL K;Q,/ii) — Uo + = I A^ e^jni-nj. (52) 

We note that we have also given the Wilson coefficient the n-labels of the operator it mul- 
tiplies. In cases where it is clear, we will only explictly label one of C or O. In addition to 
the expected tree- level amplitude term in brackets, we also give the RG kernel, f/^L ' ' ^^'^ 
an angular phase-space cutoff, Qs2 ■ We discuss each of them in turn. 

The former comes from running the SCETi operator OJ from Q to the scale /^i ~ XQ. 
When t/LL refers to an operator where all collinear directions are distinct, we will drop n's 
from the notation. From Eq. ( l20l) . we have 

[/if'°)(no;Q,^i)=A,(g,/ii), (53) 

where LL refers to the fact that we take the one-loop cusp anomalous dimension, which 



resums the leading logs of this running. As mentioned in Sec. IIIBt 24l . 1251 ] showed this 
resummation to be equivalent to that of no-branching Sudakov factors of CKKW-L. We 
discuss the running of our operators in more detail in Sec. IIVCI 

The phase-space cutoff G^j [ni ■ n[] encodes that rii ■ n[ < A^/r/"^ (the power of r]'"^ was 
discussed in Sec. IIICp . The SCET2 operator, O2 {ni,n[), can only distinguish that the 



quark and gluon are not collinear in SCET2, but does not know that they were collinear 
in SCETi. Thus, we put a cutoff on how far apart they are using rii ■ n[ to ensure that 
this SCET2 operator cannot create them in a region of phase-space where they would have 
been non-coUinear, even in SCETi. As an example, we could choose O to be the usual 
step-function 



Qs^lm-rij] 



1 Ui ■ Uj < 4 

rii ■ rij > 6k 
esk = l- Qs, ■ (54) 

In practice we will use a smooth version of the above step. For later convenience, we 
defined the complement, B^^.. In working with SCETj operators, we relate Sk to A. In 
general, the Wilson coefficient in SCETj has to encode whether rii ■ "^j < A^*-*"^^''?'^ or 
Hi ■ Uj > X^^"^'^^ /r]'^, in order to do it we will set Si = A^*~^/?7^. This satisfies the necessary 
criteria since A^*~^ <^ A^*~'^ ^ 1 (and recall that rj is the parameter that accounts for the 
decrease in p of a daughter relative to its mother). For Cg/o above, this means §2 = X/rj^. 
At the end of Sec. Ill CI we discussed how A gives us a way to parametrize strong-ordering 
and deviations from it. To this end, we did not need to assign it a numerical value beyond 
A ^ 1. Here for the implementation, we do have to make an explicit choice as to where 



^ See Appendix [D] for more detail on this matching. Though we have written Eq. (j52p to look as much like 
the SCETi amplitude as possible, we can rewrite it purely in terms of external momenta, as in Eq. (|D13[) . 
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FIG. 7: Plot of the smoothed ©-function, @{x)\^a, defined in Eq. ()D2ip . taking A = 0.1 and 
a = 0.016. The parameter A determines the value of x where the function switches from to 1, 
and 2a is the range in x over which the transition is made. Comparing this smoothed G-function 
to G^ in Eq. (j54p . we have parametrically A ~ 5^ and a <C (5^. This plot is for the case 5^ = 0.1. 



our G functions turn over, and for this purpose we will use fixed values such as A = 0.1 
and rj = 1/2. This means 62 = 1.6 and since the 77"^ is a common overall factor that all 
Sk>3 ^ 0.16. The smoothness of both G and our physical processes gives us great leeway in 
the choice for A, and we expect that any A ~ 0.1 will suffice (c/. Fig. [TOl) . 

Once we square and integrate our operators, we have certain practical considerations to 
take into account. For example, it is better to use a smoothed step. We give an example 
of such a function in Eq. (1D21I1 . and plot it in Fig. [7l where we choose an appropriate 
numerical value for 6k- If one only wishes to recover the LL shower, then one should use 
G = 1, as the errors induced by this do not affect the leading resummation. Furthermore, 
taking G = 1 ensures that the LL shower can cover all of phase space. Once we include 
corrections, though, then it is important to keep different types of coUinearity distinct and 
include non-trivial G's. In the presence of corrections, there will always be amplitudes with 
a G and others with a Q, which together cover all of phase space (see also Fig. [TOl) . 

Unlike standard SCET, where all the coefficients are of order A", CgLQ has an over- 
all weight of A^^. We get A^^ from the SCETi propagator, l/gg- The numerator is 
proportional to A and comes from the vertex: (tIq + (^i)no±7°'_L/^i ) • The second term 
is straightforwardly 0{X) from (^ )no±- Since n^ gets contracted with -B", _,_, it only 
contributes its perpendicular component in the n[ frame. From Eq. (150|) . we see that 
ino)n[± ~ no-n[ ~ (fci)„o±/A;i ~ A. 

C2 Lo is gauge invariant despite the presence of the G^^. function. This follows from 
writing Eq. (152|) only in terms of scalar products of n vectors, (cf. Eq. IDlSp . since collinear 



directions are invariant under collinear gauge transformations 23 



^(1) 



We note that we can obtain C2LQO2 

^(0) xo(0) 



^'^^n„n[] 



from the original two-parton operator, 



C{ IqO\ (no), in two steps: first we multiply it by the running factor 



(55) 
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where the formulas for f/^L are given in Eqs. fll8H19l) . Secondly, we apply the replacement 
rule 



(Xno)i^ (CLoK))ii(XnJi^i3°;±, (56) 



where c^q is: 



CloK) = — Uo + = ] —r ^s^lni ■ n^] ■ (57) 

% \ 9i y 4 

The relation ( l56ll is the operator statement of splitting in the parton shower. The scale 
Hi defines the endpoint of running in the UV theory. As we evolve down, more partons 
become apparent. We can see this here by the presence of two fields where there had been 
one. It makes the basic aspects of the shower manifest. The replacement rule affects the 
quark alone, and so we see that the amplitude for splitting factorizes off from the rest of the 
process. The RG kernel reflects the no-branching probability. Lastly, we can interpret the 
vertex portion of c£o as the "square root" of the splitting function. The spinor projector 
(fl^g/4) in Eq. ( 1571) rotates the spin-sum from ^-^ to ^g in accordance with Eq. ([2]). The 
remaining part of c^q after stripping off the G^j is: 



^/C^ 



which squares to a trivial Dirac structure. Furthermore, even though PaiQ/lo) ¥" Pa because 
of the RPI rotations we performed (where p is defined in Eq. |25|) . we have |pP(^o/?o)^ ~ l-^P 
with respect to the gauge polarization sum, da/3, so 






\P\' = -i^- (59) 



"11. 



Just as before, including the 2;-dependence from the measure and spin-sum, we recover the 
the standard splitting function oc {l + z'^)/{l — z). Thus, c^q weights the probabihty assigned 
to the expectation value of C*2lo^2 (^i?""-!) appropriately. 

Having computed the LO result for a single gluon, it is straightforward to proceed to an 
arbitrary number of emissions. In SCET2, we know that a two-gluon process comes from 
the T-product of the Lagrangian with 6*2/0^2 (^15 ''^i)- Similarly to before, the amplitude 
has the contribution, 

Af^' = C2o(0| JdxT{CscETMOi'\nr,n[)}\qn,gn,9n',qn) • (60) 

The vertex for gluon emission in the SCET2 Lagrangian is identical to that in SCETi. 
Thus, integrating out the parent of the Lagrangian-emitted gluon, we obtain a two-gluon 
SCET3 operator, C^^^qO^ {n2,n[,n'2), similarly to before. Also like in the matching 

SCETi -^ SCET2, we can obtain C^%0^^\n2,n\,n'2) from the SCET2, C^2%'^2\nu K) , 
by multiplying it by the running factor for O. 



(1) 
2 ) 



f/i;^) = A,(/io,^2)A,(/ii,/i2)i/2^ (61) 
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FIG. 8: (Left panel) Our kinematic convention for a strongly ordered process. Quark momenta 
are denoted by qi and gluon momenta by ki. (Right panel) Power counting of the LO coefficient 
in SCETjy. The powers of A with negative exponents refer to the propagator contribution to the 
amplitude. Those with positive exponents refer to the perpendicular momentum of the gluon with 
respect to its parent, which appears in the SCET vertex Feynman rule. 



with /^o ~ Q ^-iid applying the replacement rule: 



(62) 
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n" + 



Q2 



- ©53 [^^2 ■ n'2] 



where n2 and n'g are directions proportional to the quark and second gluon momenta, defined 
in Eq. ( JB15I) . and ^3 = X^/r]'^. One can iterate this procedure to obtain the LO result for 
(A^ — l)-gluon emission. If we use the replacement rule A^ — 1 times we go down to the 



SCET AT operator C^^lq C^ ^\nN_i,n[, 



i^N-l)i 



after which Lagrangian emissions are 



no longer distinguished as separate particles. We have: 
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The variables for A^— 1 emissions are illustrated in Fig. [HI where qk-i = [qN~i + X]j 

and 6k = \^^~^ /rj'^. From the power counting one knows that /iq = Q, and yU^ ~ Q^^^ where 
the latter scaling determines how /i^ depends on p-^ momenta, but not how it depends on 
ratios of the large Qj momenta. To sum LLgxp the approach taken by CKKW and elsewhere 
is to use /x^ = k\, namely the transverse momentum squared of the emission [29|, IGSMTI . 
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This accounts for soft interference effects and coherent branching, see Ref. [69[ for a review. 
To investigate this scale choice in the SCET, framework requires an examination of the logs 
in the one-loop matching computation for clo, and consideration of soft gluons in SCETj 
and SCETj+i. Having not carried out this computation ourselves, we rely on the previous 
literature. For our variables using Eq. (|38|) this implies 



2 i'qkkk\'^ \nk ■<! 



/^I=(P) ^V^- (64) 

\qk-i'' ^ 

(In contrast, the choice of invariant mass q\_^ would have yielded ^\ = {qkkk)\nk ■ ""-fcl/S, 
but this leads to incomplete cancellations of soft divergences, and therefore problems with 



the resummation of soft logs (681].) The directions Uk and n'^, are aligned with the exter- 
nal quark, g^, and the gluon momenta, kk- They are related to Uk-i through an RPIfc 
transformation. We can extend the argument to calculate the scaling of C^ Iq to the 
SCETtv coefficient in Eq. (!63|) . Counting the contributions from the tree-level terms, c^q, 

CiS^ - Uf ' l/A-^ = A-^(^-l)/^ cf. Fig. El 

Similarly to the discussion above Eq. f l59p . we can extract the vertex part of d^Q to define 
P"fc. We get that: 

|P-P = ^1±^, (65) 

[qk)nk-i± 

where Zk = qk/qk-i- Thus, the amplitude squared goes like the factorized product of the 
appropriate 1 — )■ 2 splitting functions. Since CJy ~ (^Af-i, ""-i, • • • , ^'n-i) i^ just built up from 
the repeated use of Eq. (1561) . we see that it requires no added information after we compute 
the first q —^ qg splitting. Thus, what we need to pass to a shower algorithm comes just from 
single real and single virtual gluon computations, as we list below in Sec. IIIICI in Table [B 
The collinear splitting needed for a LL shower is entirely handled by the replacement rule 
in Eq. ^7 

Lastly, we note that at higher orders in SCET^v, we will only ever need to compute the 
Wilson coefficient, C^ , of CJy ~ . Since each field in this theory has its own direction by 
the physical resolution constraint, we can use RPIat to make all operators with Vn± equal 
to zero. 



^ It is straightforward to see that we do not have additional contributions at LO in A. Firstly, consider the 
possibility of operators that do not take the form of a single-field replacement rule. These would depend 
on the details of the hard process that produced the quark in the first place and could threaten the 
factorization of the shower. In fact, we will get such terms when we match QCD — > SCETi, but they are 
always suppressed, as we discuss in Sec. lIVI Returning to single-field replacement, let us consider matching 
SCETi -> SCET2, as results in this case will generalize to all SCET,. Rule ^ sends Xm -^ C'-^^jj. Xna- 
At LO, we cannot get such a replacement involving multiple gluon fields, Bn ±, as this implies that we 
have integrated out multiple, hard {^ QA^) propagators. Such a contribution would not be strongly 
ordered, and is suppressed. In Sec. IIV( we will also see that we do have such contributions at higher 
orders. 
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B. Soft Emissions 



SCET describes soft degrees of freedom using soft quark and gluon fields: qs{x) and As{x). 
In this work, we focus on fully differential cross sections where we can always distinguish 
collinear and soft modes. In an integrated cross section in SCET, we have to implement soft 
emissions with some form of zero-bin subtractions 72(| to avoid double counting between soft 
and collinear radiation. (In the shower literature a proper treatment of softs is also often 



implemented by subtraction methods 391, |40|, l43l-l46 
sector couple through the covariant derivative. 



731.) The collinear sector and the soft 



iD^ = id^" + gA", 



(66) 



acting on the collinear fields. At LO in A, the collinear particles only couple to the n ■ Ag 
component of the soft gluons and the soft-collinear factorization guarantees that we can 
absorb this interaction into a Wilson line, Y{x), along the direction of the collinear particle, 



.0 

yn{x) = Pexp ig / dsn-As{x + sn) 

J —QO 



(67) 



In SCET, this is accomplished by making field redefinitions 23|, so that the new collinear 
fields no longer couple to soft gluons through their kinetic term, as we review in App. |Al 
The outcome for the composite fields considered here is that 



Xn ' ^nXr, 






(68) 



Note that here we consider nonabelian soft interactions, which is why the soft Wilson lines 
do not cancel for the i3^ field. 

In matching SCETj to SCETj+i, we will only consider external soft modes in SCETj+i 
with momenta k ~ QA^^*^^^ These are contained as a subset of the softs in SCETj. We do 
not consider particles with soft momenta k ~ Q\^^ that could not be encoded by onshell 
modes in SCETj+i. Such modes are forced to have larger momenta than the soft fields in 
SCETj_|_i, and they are not responsible for IR divergences. Any contributions from momenta 
of this type can be encoded in the Wilson coefficients of our SCETj+i operators. 

In a given SCETj, after making the field redefinition, the effect of soft gluons is encoded 
by Wilson lines Yn in the operators, with the form 



TV 






Y^ Y\ Y > E^^^°'''Y'^ V Y- 



Xf. 



(69) 



fc=i 



The angular ordering property and the coherent parton branching formalism for soft emis- 
sions with multiple hard partons emerge naturally from such operators in SCETj+i. If we 
take the Fourier transform of Yn{x) we get 



m=l pemis. 
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n- At^ ■ 
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(70) 
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where ki, k2, ... kn are the momenta of the gluon fields. The eikonal structure of (1701) leads 
to angular ordering. If a collinear particle with momentum qi in the rii direction emits a 
soft gluon of momentum kg, the amplitude acquires a term proportional to 

Fsoft = ^ = ^ + 0(A), (71) 

f^i ' ^s Qi' i^s 

where Eg is the polarization vector of the soft radiation and gf = qn^jl up to power 
corrections. \i AJ^q^^ q^^ ■ ■ ■ , qn) is the amplitude to emit n collinear particles with momenta 
?i) Q'25 ■ ■ ■ -.Qn and v4„+i the amplitude with one more emission, kg-, in the soft region, we get 
^n+i(gi, g2, ■ ■ ■ , <?«, ^s) ~ ^n(gi, g2, ■ ■ ■ , 9n) I]r=i C** Qi ' ^s/*?* ' ^' whcrc Ci is a color factor. 
For the cross section this implies 

(ia„+i = ^^"-^^^ 2^ C-^jPFij , (72) 

where ^^2^ and Eg are the element of solid angle and the energy of the emitted soft gluon, 
and Cjj is a color factor. Here 

Q'j ■ "^s Qi ' ^s 
is known as the radiation function. Without color weights, the integration of Wij over 
azimuthal angular variables would imply that soft gluons only contribute when the gluon is 
confined to the cones centered in the directions of particles i and j, and are hence angular 
ordered. 

To see how coherent branching emerges, we consider effects encoded by operators with 
exactly the same collinear field content in SCET, and SCET,;_|_i. Graphs involving soft gluons 
will agree, and there is no contribution to the matching. If we consider instead the collinear 
calculations that lead to the LO replacement rule Xno ~^ ci^oXni^n' ' then the soft gluons 
are encoded by 

SCET,: XnoYl, SCET,+i: CLoXnXyn[B^'X[ ■ (^4) 

For soft gluons at wide angles relative to no, rii, and n[, the effect of attachments to Y^^Y^^t 
are power suppressed because soft emission from these two lines cancels up to terms that 
are power suppressed by rii ■ n\ ~ y?"^ jif. The remaining attachment to yj, looks the same 
as those to Y}^^ at leading power, since no ■ n!^ ~ y?"^ jrf. Thus, wide angle soft gluons do not 
resolve the substructure revealed by matching to SCETj+i and effectively only couple to the 
overall color charge of the parent quark Xn^- Soft radiation that is close in angle to nx and 
n']^ resolves the split into quark Xni and gluon B^i , compensating for the n\ ■ n'^ suppression 
by additional collinear singularities in its propagator factors. Thus, the coherent branching 
formalism for soft gluons emerges naturally for amplitudes in our SCETj picture. 

From the SCET point of view, it would be natural to distinguish soft and collinear 
radiation in the shower and treat them independently, being careful not to double count. 
For simplicity, all available shower codes treat them in a simultaneous fashion. Accounting 
for soft coherent branching in the shower typically leads to modifications of the Sudakov 
probability factors (see for example Ref. 12^]), and affects the choice of evolution variable 
or adds additional vetoes. In the context of SCET, the implications of this were discussed 



recently in 75 . 
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C. Summary for LO Parton Sho^ver 

In Table [U we summarize results for the mapping between the LL parton shower and our 
SCETj picture at LO in A. In the first column, we put the elements needed for showering, and 
in the central column the translation to elements in the SCETj setup. The usual splitting 
function is related to our replacement rule Xno ~^ ci^oXni^n' ' ^^^^ ^^ ^^^^ ^^ related to the 
SCET2 coefficient of the operator O2 ■ The LL Sudakov comes from LL running factors 
related to the one- loop cusp anomalous dimension as in 24, l25|]. At leading order, soft 
emission in SCETj is taken into account by adding soft Wilson lines Yn into our operators. 
This leads to angular ordering and coherent branching, which must be accounted for with 
modifications to the shower to account for the soft singular regions. Finally, showers are 
constructed with different choices of evolution variables and the choice effects the structure 
of power corrections. In SCETj, we have seen that we can write all coefficients in terms of 
the large momenta (g) and dot product of n vectors {rii ■ Uj), which are natural variables in 
the SCETj picture. One can convert these variables to k^, virtuality, or angles as desired. 
At LL this translation is straightforward. 



Shower Concepts 


Quantity in SCET, 


Found In: 


Splitting function 


Replacement rule 


Eq.dMD 


LL Sudakov factor 


One-loop cusp 
anomalous dimension 


Eq. (poll 


Soft emission 


Soft amplitude 


Eq. dMD 



TABLE I: Mapping between parton shower and SCETj at LO/LL. 



IV. SCET POWER CORRECTIONS TO THE SHOWER 

As we have seen in the previous section, we reproduce the usual parton shower by match- 
ing collinear gluon emissions to increasingly lower-scale EFTs, the SCETj. Our goal is to 
catalog the leading power corrections (in A) to the differential cross section for the emission 
of an arbitrary number of collinear gluons to a quark. By this we mean all amplitude terms 
to LO(A) and NLO(A), as well as those at NNLO(A) that can interfere with LO(A). As we 
will argue in Sec. IIVDI in most cases of interest, there is no LO(A)/NLO(A) interference, 
and so we focus on the most important power suppressed terms which are NLO(A)xNLO(A) 
and LO(A)xNNLO(A). Just as in the strongly-ordered case, it is convenient to integrate 
down to SCETj+i when describing the emission of z-gluons. We obtain these corrections by 
doing our matching computations at higher order. We will show that there are two distinct 
types of subleading matching, and they have a different physical interpretation: 

• One type originates in matching QCD — )■ SCETi at higher orders. This generates a set 
of subleading terms that remain suppressed as we move down to lower-scale SCETj's. 



31 



We call them hard-scattering power corrections as they involve the details of the hard- 
scale process that created our original partons. Also, they are most important for 
partons radiated closest to the hard vertex. 

• The other type comes from the subleading matching SCETj — )■ SCETj+i. They in- 
volve processes described by the SCETj Lagrangian, but ones that get integrated out 
into higher dimension operators at lower scales. These corrections are ubiquitous. 
They do not depend on the hard-scattering details, and we can determine them for 
arbitrary SCETj — )■ SCETj+i once we have found them in SCETi — )■ SCET2. Fur- 
thermore, they relate to known 0{as) corrections to the q -^ qg splitting function, 
which exponentiate to sum part of NLL. For this reason, we call them jet-structure 
corrections. 

Determining the above to NLO(A) in the cross section will only involve single and double 
gluon emission. Thus, we will never need to compute in a lower-scale theory than SCET3. We 
perform all the necessary QCD— )-SCETi— )-SCET2— )-SCET3 matchings for these amplitudes 
in Appendices OlEl Below, we discuss the final results for the corrections, with Sec. IIVAI 
focusing on hard-scattering and Sec. lIVBI on jet-structure. For these portions of the paper, 
the matching is only done at tree level, though formulas in the Appendices include one-loop 
RG kernels. We give the effects of LL running on correction terms in Sec. IIVCI along with 
a discussion of how to include NLL resummation for the LO (in A) Wilson coefficients. In 
Sec. lIVDt we will study the amplitude squared and will see there is a great simplification of 
the interference structure in SCETjv, and hence for NLO(A) power corrections in a shower. 
Lastly, we give in Sec. lIVEI the NLO counterpart to our LO table in Sec. IIIICI We describe 
how our corrections from subleading operators relate to improvement of the parton shower 
with higher order resummation of logs, corrections at higher order in as, as well as corrections 
to spin correlations and interference. These effects are summarized in a shower reweighting 
formula, Eq. f ll04p . 

A. Hard-Scattering Corrections 

Just as in Sec. IIIIAt we begin by examining the matching QCD— j-SCETi for single gluon 
emission coUinear to the quark. For this case, all corrections are of the hard-scattering 
type. Beyond LO, we can have dependence on the process that creates the qq pair. For 
concreteness, we will consider the coupling of QCD quarks to the vector current, Jqqy) = 
qj'^q. The matching is performed in the center of mass frame with the initial virtual photon 
having momentum, p^ = {Q, 0,0,0). The full details of this matching calculation for QCD 
to SCETi are in Appendix O To reproduce the full QCD current, Jqcj), we need an infinite 
tower of SCETi operators increasingly higher order in A. However, to get the required 
amplitude to NNLO(A), we only need four: 



ATnNLO = C'i°Lo(^o) I dx{0\T{CsCETAx)Of^}\ qnognoqn) 

+ Ci^NLo("'0,'^o)(0|C!^^| qno9noqn) + C[]l{no,no){0\Ti^^ \ qno9noqn) 

+ C«(ni,n;)(0|O«|g„,(7„.g^), (75) 
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FIG. 9: Matching QCD to SCETi to SCET2 for one gluon emission which is either colhnear to 
the quark or is in its own direction (SCET graphs for emission colhnear to the antiquark are not 
shown). The figures represent operator structures that describe this process in each of the three 
theories. The QCD contribution is standard. In SCETi, we either emit a colhnear gluon through 
the time-ordered product of the Lagrangian with an two-parton operator, or from three-parton 
operators. In SCET2, the emission relevant for us only arises from higher-dimension three-parton 
operators. 



where 



)(0) 









XnoXn ; 
Xno 9^no±Xn 



Ti '[no, no) 

0^^\m,n[] 



X' 



no 



7/3 



' nn-L y «( 



no J 



Xf. 



Xm 9l3n[±Xn 



(76) 



-(1)/ 



Here we introduced a short-hand for the notation estabhshed in Eq. ( fT4l) . 7^ '[no, no) 
0[ ' ' {no,no )■ We give the expression for C| lo i^ Eq. ( I42l) . The amphtude from the 
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operator 0\ (no) is shown in the first diagram in the SCETi column of Fig. [9l those from 



i(i), 
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l(l). 



0\ {no, no) and 7^ {no, no) in the second, and that for 0\ '{ni,n[) in the third. 



)(i). 
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We call 0{ {no, no) and Ti'{no,no) "two-jet" operators as they are labeled with two 
distinct colhnear directions (no and n) (we do not denote the antiquark direction explicitly, 
following the convention in Eq. H3|) . They describe a gluon colhnear to the quark. We 
obtain the coefficients C{ {no, no) and C[ ■^{no, no) by expanding the QCD amplitude in the 
limit of small gluon momentum transverse to the quark's direction with the usual SCET 
proportionality: (no ■ ki,ki, kino±] 



{X^,1,\)Q. C['>{no,no) and C| 



(1), 



ni,n^, 



are derived 
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above Eq. (1C13P and given here: 

^i!r(^o, no) = -^7„^,x7fox7„"oX - ^ /^T^ox • (77) 

We use the same kinematic variables as in Fig. [51 For C| (rio,no) and C|^(no,rio), the 
initial current is not a spectator, so neither term is simply proportional to the 7^^ with 
which we started. This dependence on the details of the rest of process is a characteristic 
feature of hard-scattering corrections. There are an additional set of two-jet configurations 
corresponding to the gluon collinear to the antiquark. These are trivial to obtain by charge 
conjugation. 

The operator C| {ni,n[) is a three-jet configuration, as it describes three distinct di- 
rections. Whenever we have an operator where each field has its own index label, we can 
choose the n^ such that they are exactly aligned with the external particle momenta. We 
give the coefficient CJ {ni,n\) in Eq. (IC16p . 

Going to SCET2 for single gluon emission is straightforward. The basis of operators 
needed to reproduce the amplitude fl75|) is equal to f l76|) . but with SCET2 fields: O2 {no), 
O2 {no, no), T2 {no, no), and O2 {ni,n[). As the computations get more comphcated with 
subsequent emissions, we wish to minimize our effort by only including those terms necessary 
to give the corrections to a shower Monte Carlo. This means we are only interested in the 
following: 

1. We will need to keep those NNLO(A) contributions that can interfere with LO(A). 
These give terms at the same order as an NLO(A) operator squared. We do not 
compute NNLO(A) amplitude terms which have zero interference with the LO(A) am- 
plitude. A list of the necessary computations is found in App. |Dl 

2. Our ultimate goal is not a complete SCETj theory from which one can do computa- 
tions, but an improved shower algorithm. In Table [U we give a list of those ingredients 
needed to construct a map between SCETj and a LL parton shower. We will augment 
the map with items needed for corrections (Eq. I104[ Table HTj) . but will not calculate 
contributions which only contain redundant information for the shower amplitude. 

The latter point has important implications for the sorts of operator structures we need to 
consider. If we wanted to do computations in SCET2, then we would need all operators 
and Wilson coefficients to the order we are working. However, single gluon contributions 
in SCET2 where the gluon and the quark are collinear (inside a cone of angle ~ A^, i.e. 
O2 {no, no) or 7^ {ni, ni)) correspond to a quark which does not split until after the scale 
of matching SCETi — ^ SCET2. The corresponding no-branching probability, however, 
is already determined in SCETi from the one-loop RG kernel. Thus, the coefficients of 
these operators in SCET2 are not required. We only need to calculate those single gluon 
contributions where each field has its own index label in SCET2, which means C2 {ni,n[) 
ioT 0^2\ni,n[). 
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The matching equation for C2 {ni,n[) in SCET2 is: 

Ci'\n,,n\){0\O^^^\q^,g^^Jn) (78) 

= C[%o{no) j d^x{0\T{CscETAx)Of^}\qn,9n'Jn) 

+ Ci^NLoK, no){0\0'i^\qn^ g^'^ qn) + C[^\ni,n[){0\O^^^\qn, gn[ qn) 
+ Ci^|(no,no)(0|7^^^^|g„, g^'^ qn) ■ 

It is convenient to decompose Cg {ni,n[) as 

where the four terms on the RHS of Eq. (!79l) correspond to each of the contributions on the 
RHS of Eq. flTHj) . We calculated in Cj lo i^ ^^- ^|52|) using RPIi to rotate objects in the 
SCETi amplitude such that they can come from SCET2 operators that annihilate the given 
external state. The second through fourth terms can be calculated in a similar manner. 
Their values are derived in Eqs. ( 1D17I) - ( JD19I) : 






(80) 






+ 



C^NLoK'^^l) = [^{^n[±V' 



1 

Q 



ni-n[)qiki . 



7"0fcK-^i], 



W^lA + n'^-^i^l ■ n'l))Tn'^. 



^,{^^^^A-n' {n.-n\)^^y^,^)QsAni-n,] . 






Here ni and n'^ are aligned with the direction of the quark and the gluon, and Vi is defined 
in Eq. (]B12|) . In Eq. (!80|) . we have left off the running factors from evolution of the SCETi 
operators. The terms in Eqs. (1761) run differently. In particular, the two-jet and three-jet 
operators have different LL evolution. Therefore, it is important to decompose C2 as in 
Eq. ( 179|) . so that we can keep track of which SCETi evolution factor to include for each. 
The running of these operators is discussed further in Section IIVCI 

We also note the different G dependence of the terms, where G and G we introduced in 
Eq. (I5^ and the surrounding discussion. We can read off from Cg ^lo ^^^ origin as a three-jet 
term in SCETi, while the others come from two-jet operators. The G functions are necessary 
because without them SCET2 operators, {e.g. O2 (^1, n[)) can only tell that the quark and 
gluon are not collinear according to the SCET2 definition. By including these phase space 
cutoffs, we can keep the distinct origins of different contributions manifest. By adopting a 
smoothed step function, as suggested in Sec. IIII Al and given in Eq. (]D2ip . the amplitude 
squared for C2 C2 "^ill be continuous despite having different supports in different parts 
of phase space. An example of this is shown in Fig. [TOl The full expression for the plot 
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FIG. 10: Merging of the two-jet and and ttiree-jet squared amplitudes using a smooth theta 
function for the 7* — t- qqg process. 



Plots of the amplitude squared components from C2 O2 ■ 



|^999|iQ + |^W9|^^Q 2_jet (short dashed green), \Ai^3\ 



red) versus \k 



NLO, 3-jet 



(long dashed blue) , and sum (solid 
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The amplitudes are evaluated without running coefficients, and taking ki/qo 



0.4. The 62 parameter in the 0-function is 1.2, which for the above p fraction corresponds to rj 
and A = 0.08. 



0.5, 



is given in Eqs. (1D24P and (lD25p . To illustrate the effects of including hard-scattering 
corrections, in Fig. Owe plot the ratios i?LO = I^'^^'^Ilo/I^'^^^Iqcd and i?NLO = (I^^^^Ilo + 
l^^^^lNLO,2-jet)/l^^'^''lQCD versus the gluon perp momentum. Here, |^^^^|qcd is the QCD 
amplitude squared for one-gluon emission, |A''^^|lo is the SCET2 amplitude squared for 
one-gluon emission from the LO coefficient Cg lo^2 (from Eq. [52]), and |^'''^|NLO,2-jet is 
the NLO (A) amplitude squared for one-gluon emission in the two-jet region that comes from 
the coefficients C2 nlq and C2 nnlo (given in Eq. [80]). As we expect, including corrections 
up to NNLO (A) in the amplitudes squared extends the region where tree- level SCET2 and 
QCD agree. The advantage of using the one-gluon SCET2 amplitude over QCD comes from 
factorization properties that effect interference as well as renormalization group evolution. 
For example the one-loop running in SCET2 performs the LL Sudakov resummation. 

With two-gluon emission, the SCETi graphs will include jet-structure corrections in ad- 
dition to hard-scattering ones. It is straightforward to distinguish the types as the former 
result from taking time-ordered products of the SCETi Lagrangian with operators generated 
by the LO replacement rule, Eq. fISBl) . while the latter will come only from terms involving 
a power suppressed SCETi operator. To fully identify the subleading contributions to two- 
gluon emission, we must match down to SCET2 where the LO contribution is ffist uniquely 
identified. We already know that it comes from two applications of Eq. (1561) . 

In Fig. [121 we show the contributions to two-gluon emission in QCD, SCETi, SCET2, and 
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FIG. 11: Plot of the ratios of the amphtudes squared for 7* 



qqg, namely Rlo 
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FIG. 12: Matching QCD to SCETi to SCET2 to SCET3 for two gluons emitted cohinear to the 
quark direction (SCET graphs for other gluon kinematic configurations not shown). Once again, 
we depict the operator structures that lead to this process in each of the theories. Gluons drawn 
away from the central vertex are emitted by the leading order Lagrangian in that theory, while 
those coming from the vertex are due to higher dimension operators. 

SCET3. The first column in the SCETi category corresponds to the jet-structure corrections 
to be considered in the next section. In the second column we have a set of hard-scattering 
corrections from taking the T-product of the SCETi Lagrangian with the suppressed single 
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gluon operators we calculated above in Eqs. fl77|) . C| 0| and CJ^Tf . 

In considering the basis of operators in SCET2 we do not need operators such as 
T2 (ni, n[), since Vn'xBn'x = 0, with n'l lying along the gluon momentum. We can use 
RPI2 in SCET2 to make a coordinate choice where they are not necessary. As mentioned 
above in the single gluon matching section, our interest is only in calculating those terms 
needed to improve a shower algorithms, which precludes us from considering operators such 
as T2 {no, uq) or O2 {no, no), corresponding to an unbranched quark passing from SCETi 
into SCET2. Therefore, for double gluon emission we only need to calculate the coefficients 
of the following operators: 

0?{n,,n[) = Xrn9l3:'^±Xn (81) 

Of\n2, n2, n[) = Xna^^naX dK'^^Xn , 
Cf V2, n\, n\) = Xn29K{l. 9^n{±Xn , 
Cf V2, n[, n'2) = Xn29K{l. 9K'^±Xn ■ 

Thus in SCET2, we are interested in two-glnon operators where two fields can have the same 
label. When we pass to SCET3, we can restrict our interest to only C3 {n2, n[, n'2). 

We already gave the coefficients of O2 (^1, n\) needed to compute the leading power 
corrections in Eqs. ( IHOj) and f lD17P - (]D19p . We get an NLO(A) contribution to the two 
gluon amplitude by computing the matrix element, C2 nlo(0|'^{^scet2C'2 }\qQ99) (first 



SCET2 column in Fig. |T2|) . The contribution receives no further suppression as the 
gluon from £scet2 gives a tree-level vertex x propagator factor of A~^, just as with LO. 
There are also coefficients we need from two-gluon matching calculations for the oper- 
ator O2 (second SCET2 column in Fig. [12]). Putting in the index structures, these 
include Q nlo('^2, n[, n'2) for O2 {n2,n[,n'2), Q nlo('^2, n[, n[) for O2 {n2,n[,n[) and 

^2^NLo(^2, ^2, n'2) + C'2^NNLo('^2, ^2, n'2) for Cg (^2, ^2, "-'1) • We iuclude NNLO(A) for the 
last one as only it interferes with LO(A). In the next subsection, we give the jet-structure 
corrections. All hard-scattering contributions to these structures just listed are beyond the 
order we need except for C2 nnlo('^2, "^2, ''^i), given by the matching equation: 

Ci!NLo(^o,'^o) dx{0\T{CscETi{x)0^^\no,no)}\qn29n29n[qh) (82) 

- C'2^NLo(^2,rii) dx{0\T{CsCET2{x)O^^\n2,n'j^)}\qn^gn29n[qn) 

= C2^lmho{'^2,n2,n'^){0\OP{n2,n2,n[)\qn29n29n[qn) , 

where we subtract the Lagrangian emission graph in SCET2 from that in SCETi (C| {^lq 
is given in Eq. [77] and C2 nlo ^^ ^^- ISH])- The result for Cg nnlo(^2 5 n,2, n'^) is given in 
Eq. (]D61|) . It is straightforward to see why O2 only gets hard-scattering at NNLO(A) and 
higher. By definition, hard-scattering has to involve a suppressed operator from the QCD 
— )■ SCETi matching, and so we begin at NLO(A) at the lowest order. Including a second 
gluon, but demanding that we cannot write it as coming from a SCET2 Lagrangian emission 
takes us to one order higher, namely NNLO(A). 
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All the contributions we have discussed so far have come from the hard-scattering, single- 
gluon, suppressed operators in SCETi. There are also those with two gluons. That is to say a 
process where neither gluon comes from the SCETi Lagrangian, represented by the diagram 
in the third SCETi column in Fig. [121 One example is double _L-gluon emission from the 
antiquark, as shown in the third QCD graph of Fig. [121 We know from applying Eq. f l56|) 
twice, that LO for this process is at 0{X~^), counting only the tree-level vertex x propagator 
factors, as these are all we need to compare different qqgg processes. We readily see that 
double antiquark emission is ~ A*^ as there are no small virtualities or emission angles for 
this term. Thus, they are N^LO, and beyond this analysis. Besides antiquark vertices, 
we also have subleading emissions from the quark in QCD that arise from the suppressed 
SCET-spinor portion of the QCD quark propagator (c/. Appendix |X|) . If both emissions 
come from the suppressed propagator, once again, this is ~ A'^ at lowest order, and so we 
can neglect it. Mixed antiquark/suppressed spinor contributions are also N'^LO. 

Thus, we do not need corrections to double emission coUinear to the quark if they do 
not involve at least one SCETi Lagrangian insertion. We can extend this argument further. 
If there are no SCETi Lagrangian insertions, then the contribution goes like (9(A°), while 
LO goes like 0{X 2 ). Thus, to the order we are working, we only need the single gluon 
hard-scattering corrections given by Eq. (jHOl) . plus Lagrangian insertions. 

B. Jet-Structure Corrections 

The jet-structure corrections only involve contributions from the SCETi Lagrangian. 
These arise from the graphs in the first SCETi column in Fig. [121 We specifically designed 
our leading order replacement rule in Eq. (l56l) . so when used twice it only contains that 
part of double emission corresponding to the leading strongly-ordered limit. This occurs 
for the gluons having collinearities ~ A, A^, respectively. However, SCETi describes other 
kinematic situations and in this section we compute the corrections from them. 

The prescription for obtaining two-gluon jet-structure corrections is to compute the dou- 
ble gluon emission amplitude in SCETi coming from two Lagrangian insertions and take 
different limits on the relative collinearities of ^2, n2, and n[, where these labels refer to the 
null vectors exactly proportional the corresponding particle momenta. We can define: 

^NLO = C[%o{no) dXidx2{0\T {CsCET^Xl) CsCETi{x2)0'^^ }\qn2gn[gn'^qn) , (83) 

and then calculate, 

lip ,^NLO = C'?NLo('^2, n[, n'^){0\O^2^\qn2gn[gn'^qn) , (84) 

, lip ,, ^NLO = C'2^NLo("'2, %> ^l) (010^ l^n^ t/< t/n^^n) , (85) 



n'j-nj^A* 






lip ,^NLO = Q Lo(^2,^l)(0|r{/:sCET2C2 }|gn2^ni^n^?n) 
n2-n2~A* 



+ C'2,NLo(^2, n'l, ^2)(0|C2 \(ln29n'^9n'^qn) ■ (86) 

We note a few things about the above equations. Firstly, there is a correction to the LO 
Wilson coefficient obtained from the replacement rule (Eq. [56]) . We cannot get it purely 
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as a limit of ^nlo' ^'^ ^^ need to subtract off the LO contribution. Secondly, the limit in 
Eq. ( 184|) does not lead to an expansion of any part of A^q, as the scaling of the n-indices' dot 
products is exactly that from SCETi. Even though it just gives back the same expression as 
the SCETi amplitude, A'^Iq, the SCET2 result for C2 nloI'^z, "^i, ^2)02 tells us something 
more. This Wilson coefficient is proportional to Q52[n'i ■ n'2]Qs2[''^2 ■ n'2], where the G's only 
have support outside the phase space region of Eq. (185!) . as well as the strongly-ordered 
limit, Eq. (186|) . (see Eqs. (15^ and ( ]D21|) for the definition of 6, 6). The full results for the 



Wilson coefficients shown in Eq. (jHSD can be found in Eqs. (ID 3 YD . (lD52p . and (lD59p . At the 
amplitude level, given a particular phase space configuration for an external state, we will 
only ever need one of these terms for double gluon emission in SCET2. Squaring the result 
is straightforward as there will be no interference between them. 

We will now examine how to improve the matching of SCETj to SCETj+i, and show that 
the jet-structure corrections computed here generalize to that case. We first notice that the 
first two operators above do not interfere with the one giving LO, as they have different index 
structures. The subleading term in Eq. (1861) does inhabit the strongly-ordered region of phase 
space, but as we will argue in Sec. lIVDt LO(A)/NLO(A) interference cancels out of most 
observables of interest. Before proceeding, we note that our description of corrections to two- 
gluon emission gets even simpler when we match to SCET3. In SCET3, the only operator 
we need has distinct collinear directions for all fields. Thus, we can write all hard-scattering 
and jet-structure corrections to two-gluon emission we have found in the coefficient, C^ , 
for the operator C3 '{n2,n[,n2) = Xn29^n' ±9^n' ±^n, as we do in Eqs. (IE2p . The same will 
hold for i-gluon emission in SCETj+i. Our NLO(A) jet-structure operators therefore have 
the following form: 

Ci%ioin2, ni, n'JOf = hf Xn^gB^'^^gB^^^J^Xn , (87) 

where hj is given by Eq. flEllI) . Here / = {1,2,3}, and we distinguish the coefficients 
Cg j^Lo depending on which SCET2 operators they come from in order to properly account 
for their RG evolution in SCET2. 

When doing the LO matching for SCETj to SCET^+i, we found that the replacement rule 
to go from SCETi to SCET2 generalized to the case of z-gluon strongly-ordered emission. 
Similarly, we can take the above operator, Eq. (187|) . and recast it as a replacement rule for 
our original current insertion, C{ lq^i • ^^ takes the form of a 1 — )■ 3 replacement rule: 

Xno -^ hf Xn^gBn'^^gB^^,^^, (88) 

with contributions from 1 = 1,2, 3. 

If we want to consider the NLO(A) radiation of i+1 gluons, we can perform a very similar 
matching between SCETj and SCETj+2 to the one above for SCETi — )■ SCET3 to obtain 
an operator Ci^2 nlo^1+2 ■ Since the first {i — 1) emissions are strongly ordered, they 
completely factor out. Thus, the amplitude for the emission of the final two gluons will be 
identical to that for simple two-gluon emission. We can therefore take the {i — 1) gluon LO 
operator, C'j ^q Cf" , and use the replacement rule in Eq. (155]) . to obtain Q^2 nlo/^1+2 • 
Our NLO(A) replacement rule corresponds to violating strong ordering at any location in the 
shower, either by taking the j*^ and (j + 1)*'^ gluons to have the same parametric coUinearity 
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with respect to their parents, kj+i± ~ kj± (Eqs. HI] and |85]) ; or by including the region of 
phase space where the propagator between them is hard even in SCETi, and so we get no 
colhnear divergence as the quark and second gluon become colhnear (186 p . ^ 

It is not difficult to see that this gives an NLO(A) contribution for any j. If we have 
i-gluon strongly-ordered emission, the tree- level factors, CLo(^fc-i), (c/. Eq. [63]) will go as 
^-t{i+i)/2^ where the j^^ gluon contributes A"-'. If we violate strong ordering as we mention 
above for any two gluons, the product of their vertices times propagators goes like X~^^ 
instead of A"*^^-'^^). Thus, we can insert Xno — ^ hi Xn29^n'^±9^n'^i. instead of two successive 
Xno — ^ CLoXni-Sn'x's in Operator matching as a "defect" in strong ordering at any stage 
and obtain an NLO(A) jet-structure correction. The G-functions contained in the Wilson 
coefficients, Ci'-^-^q , allow us to read off at which step in the shower we violated strong- 
ordering. 

In App. [F] we show that an integrated version of h^ is related to the splitting function 
at NLO in a^, which serves as a cross-check on our computations. 



C. Operator Running 

Up until now, our discussion of matching has taken place mostly at tree-level. Connect- 
ing to the no-branching probabilities and log resummation in the parton shower however, 
requires that we include the anomalous dimensions needed for running. For this reason, 
our final expressions for Wilson coefficients in Apps. ]C]lE] include the necessary notation for 
evolution kernels. Identifying the power suppressed amplitudes as corresponding to pertur- 
bative corrections to more inclusive observables, it is natural to take only LLexp evolution for 
power suppressed or a.^ suppressed corrections, and include NLLgxp evolution only for the 
leading shower terms. For the former, we assume (without carrying out the proof in SCETj) 
that we must make the k^j^ choice for the scales /i| as in Eq. fl64l) . and that this accounts for 
the difference between LL and LLexp- NLLexp would require full one-loop, two-loop cusp, 
and NLL a^ running, plus any modifications to the evolution induced by subleading soft 



effects. If subleading soft effects are neglected then in the terminology of 60|, |6l|] , this gives 
the full collinear NLLexp resummation. The subleading logarithms coming from pure soft 
effects involve the exponentiation of nonabelian matrices. As mentioned earlier, we do not 
compute the effects of subleading soft SCETj operators here. (In fact, for more than three 
hard, colored particles, the problem is quite non-trivial [60|.) 

In this section, we determine the LLexp running for our subleading operators and discuss 
what is missing in our setup for a NLLexp evolution kernel for emission anywhere in the 
shower. To set the stage, we consider SCETi matched to QCD at the scale Q for the ffist 



At this point, one may ask why we do not go farther and consider the case kj^i±^ 3> kj^. In fact, we do 
not have to. Since the amphtude for i-gluon emission has an underlying Bose symmetry, we are free to 
partition phase space into i\ regions, each of which gives an identical contribution to the cross section. 
Thus, to get the final answer, we only need to integrate over one of them. While we can choose this 
region such that fcj+ij_ ^ kj±_ never occurs, we are forced to include fcj+i± ^ kj_\_- If we do not wish to 
partition phase space in this manner, then the Bose symmetry implies that the result for fcj_|_i_L 3> kj± 
can be obtained from the configurations already discussed. 

41 



order power corrections. We then run down to /i in preparation for matching to SCET2. 
The zero and single gluon operators in SCETi acquire the following running factors, U, {cf. 
the tree- level version in Eq. (I77|) ): 



Q 

qiKi \ giv / 



C['\m,n[) = -f/(2-i'°)K,n;; Q,/.) ( ?--^7„Vl>.7^ 



' 7r^^ - Qi ^iTJ + / Z'^sl T < 



n-pg)ki\ ^T / (ni-n'Jgi/ci 



't^ox)' (89) 



+ 

where the superscripts follow the convention in Eq. flT^ . We inserted the symbol in the 
second and third line of Eq. fl89|) since an operator with multiple fields sharing the same 
collinear direction can convolve the momentum fraction of p between the corresponding RG 
kernel U and momenta in the tree- level coefficient. This is because collinear fields that are 
in the same direction in SCET can exchange momentum while running down from Q to /x. 
The anomalous dimension of an operator is independent of which SCETj it is defined, but 
does depend on the field content and in particular how many different collinear directions are 
in the operator. Thus, the RG-kernel for Xno9^no±Xfi ^^ different from that of XnigB", j_Xn- 



In Ref. [2^, |25|], the LL part of U^^'''~^''^\Q, /i) was related to the Sudakov form factor, 
Eq. (1201) (up to accounting for the soft effects of angular ordering [75|). The cusp term 
in the anomalous dimension resums the LL, and comes from soft and collinear one-loop 
diagrams. The result from the soft diagrams is constrained by that of the collinear diagrams 
in order to cancel out infrared sensitivity that cannot be absorbed in local counterterms at 
the hard scale. Here we will use this same argument, but in reverse, in order to determine 
the LLexp anomalous dimension of various subleading operators. 

Due to the soft-collinear factorization, the soft structure only depends on the number of 
collinear directions. After making the field redefinition, operators like XnoXn and Xnofl''^no±^'^ 
both have Y^^Yn, and so both have the same soft divergences. Hence they have the same 
one-loop cusp term and the same LL anomalous dimension from the sum of collinear and soft 
loops. Thus, the leading-log resummation only depends on the number of collinear index 
directions in the operator, and not on the number of active partons. (At leading power these 
concepts are the same, but it is not so for the power corrections.) We therefore have 

US"''^ (no) = f/if'^^ (no, no) = t/g^^^ (no, no) , (90) 

where we give [/ll ' i^ Eqs. ( IT8l) and ( IT9l) . Thus, at LL order we have the full set of 
evolution kernels for subleading collinear operators, and we account for these factors in the 
appendices. Since this is a LL effect, we expect soft radiation and angular ordering to be 
incorporated in a manner identical to the evolution factor in the LL shower. 

An important consequence of this result for the LL evolution is that it justifies treating our 
hard-scattering corrections as improvements to the fixed-order, matrix-element calculation 
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that goes into a shower algorithm. Correcting the two-jet amphtude with either C| j^lo °^ 
C| ^, we see that the LL resummation is the same as that in the standard shower except 
that there is an extra parton already inside the leading jet. We thus get a shower correction 
just by using a matrix element improved by including our hard-scattering terms. This is 
unlike simply running a LL shower on higher order matrix elements, as different anomalous 
dimensions control different operators' evolution. Some, like those just mentioned with only 
no and n collinear directions, run like two-jet configurations, that is with a quark-antiquark 
Sudakov. Others, {e.g. C[ 0[ {n[, rii)) have three-parton running since they have three 
distinct collinear directions. This latter set corresponds to the usual implementation of fixed 
order corrections in parton showers, but the former is a novel type of shower improvement. 
On the other hand, the effect of jet-structure corrections is not to modify the initial 
scattering process, but to go hand in hand with the NLL change to the leading operators' 
running. Similarly to Eq. ( l63l) . we might anticipate the following Wilson coefficient for 
OJy ~ with evolution: 



.{N-l) 



^N,NLL,li^) 



N-1 



xU^"l-'\fXm-l,f^m)cl^{nm-l)] T^ (91) 



with a sum over all locations where the NLL evolution can be inserted: 

N-l 
^Ar,NLL, 1 = ,^ ^AT.NLL, l("^)- (92) 



m=l 

.(7V-1) .MN-1) 



One would expect to use C}^ j^ll i^n along with our real emission corrections (Eq. IE13P 



to correct a shower to resum at NLL the ratios of all emission scales (c/. Eq. 11041) . The 
complication we face for the calculation of f/^^L is that this correction to the evolution 
kernel must, in principle, be carried out in the same scheme used to distinguish the phase 
space regions for the jet-structure corrections, and hence can depend on the choice for the 
Q functions. In particular, we could have non-trivial operator mixing on the edge where 
the cutoff makes a smooth transition between operators with different numbers of jets, 
and we have not yet performed the analysis that would determine whether this affects the 
resummation at NLLgxp order. Furthermore, it is possible that power suppressed soft effects 
will also have implications for the subleading evolution kernel, and may make the nonabelian 
generalization of Eq. ( I92l) tricky. Our lack of an appropriate NLLexp evolution factor for the 
shower is due to these two issues. 

To setup the distinction between kinematic regions, we used Wilsonian type 6 functions, 
but from the point of view of evolution MS would be simpler. Although it is only indirectly 
relevant to our setup, it is nevertheless still interesting to consider how the NLL evolution 
kernel would arise in MS. As we discuss below in App. (Fj when integrated over phase space 
in dimensional regularization the jet-structure corrections give the real emission portion 
of Pqq , which is the 0{as) correction to the Altarelli-Parisi splitting kernel. Combined 
with known SCET results for single-emission at one-loop, we can recover all of the abelian 
portion of Pqq . Obtaining this expression is important conceptually. It validates our formal 
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expansion in A, showing that corrections to (9(A^), along with a set of known one-loop 
diagrams, capture contributions needed for coUinear NLLexp resummation. On the practical 
side, it provides a cross check on our computations. 



>(i) 



With Pqq in hand, we can extend the argument of [2J, |25| that the Sudakov factor 



gives the LL part of the the RG kernel U^'^'''^^{Q, ft) (Eq. [20]) to the NIL level, looking at 
[/(^'°'°)((5,/i) for running of the operator C| IqO[ . Using the Sudakov factor of \ZI] for 
quarks, we have: 



Aq{Q,fj,) = exp 



-^ I ^"^(^') 



dfjf 
7 



« J 1 + 2;^ 

az 

V7 1- z 



(93) 



where we recognize Pqq , Eq. ([3]). Performing the z integral and expanding in the limit of 
large Q gives: 



Aq{Q,fi) ^ exp 









'og|^ 



+ 



(94) 



which is identical to U^'^'^''^\Q, fj.) at one-loop. The term in the exponent proportional to 
log(/i'^/(5^) sums the leading logs in the parton shower. We also see that upon /i' integration, 
we get the double logarithm characteristic of the soft-coUinear divergence of collinear split- 
ting. Interpreting Eq. (p^ as an RG kernel, this log piece is coming from the one-loop cusp 
anomalous dimension, Cp. The factor of 3/2 is the remaining part of the one-loop anoma- 
lous dimension, and it sums part of the collinear NLL.^ In order to get the full collinear 
NLLexp summation, one also needs corrections corresponding to the two-loop cusp anoma- 
lous dimension. This is a known result in SCET for the operator XnXn, which we can relate 
to Pqq , by adding the subleading splitting function to the exponent of Aq{Q, /i). We wish to 
stress, however, that the ultimate goal of improving parton showers through resummation 
is to include all next-to- leading- logs. ^"^ In this paper, as mentioned previously we have not 
considered the effects of soft NLL, nor those related to the two- loop running of a^, which 
will affect collinear NLL. Our formulas in Apps. ICllEl include LL running for all subleading 
operators. In App. [F]we discuss the relation of our 1 — )■ 3 splitting amplitude with Pqq' in 
MS. The collinear-NLL- improved Sudakov corresponding to this is 



NLL 



(Q,/i) 



exp 



dfi' 



1-- 



dz [Pi^\z,as{f^')) + Pil\z,as{fi'))] 



(95) 



where Pqq given in Eq. (|3| 



and Pqq in 



62|. Once again, we integrate in z, expanding in 



^ Since Eq. ([94)) resums the NLL contributions expanded in the cross section (c/. Eq. [5]), Ref. [27| calls it 

the NLL Sudakov factor. 
^° At a practical level, while we see full collinear NLLoxp as coming from a straightforward extension of this 

work, pure soft NLL may only be possible at the leading orders in 1/Nc 60|, [61|. 
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large Q to get: 

Afi^(Q,/i) = exp 






-'(^'Xc, L, iO H- ' 



TX " V VQV 2 



+ -^C, [C, (y - yj - yC7,T.n.j log ^^ 



(96) 



where the term oc a1 reproduces the known result for the two-loop cusp anomalous di- 
mension. While including this MS NLL effect for "no-branching" was already possible, our 
result in Eq. ( 1E13I) allows one to modify the differential cross section for real emission to 
include the effects of Pqq\ as well. Without including both, one does not have a systematic 



improvement beyond LL. In [56|, the authors were able to get NLLcxp soft resummation 
by treating the subleading real and virtual effects in semi- inclusive observables for DIS and 
Drell-Yan. A full implementation in our framework with more exclusive observables must 
wait for computations that address the missing NLL ingredients mentioned above. 



D. Squared Amplitudes and Interference Structures 

As discussed previously, our series of matchings terminates with SCETjy, where each field 
has its own index direction. Further Lagrangian emission from these operators is physically 
meaningless, as the resolution scale is set ~ O(GeV), below which we stop computing in 
perturbation theory and pass to a hadronization routine. Thus, we match everything to the 
single operator O]^ {nN,n[, . . . ,n'j^_^) and all the information about the shower at LO 
and NLO is encoded in the Wilson coefficients. In this SCET^r, we square amplitudes and 
compute corrections to observables, as we detail in Sec. IIVEI As we saw in Sees. IIVAI and 
IIVBI for arbitrary A^, we only needed one and two-gluon computations to obtain leading 
corrections in A to the differential cross section. Using the LO replacement rule (Eq. 1561) will 
account for the rest of the multiplicity. Since the strongly-ordered emissions it describes have 
trivial interference, we should expect that squaring our results retains the simple picture we 
have for corrections at the amplitude level. 



1. Interference for LO^ and for Jet-Structure Corrections 

It is a general statement about SCET fields with different n index labels that they have 
no overlap in Hilbert space. As an example, we can take two different operators, O^ and 
On2 where all the fields in 0^ and (9„2 are identical, except those labeled by rii and n2 
{e.g. Xni versus Xn2)- For generality the field labeled by n2 may or may not be in the same 
equivalence class as rii. We thus have:^^ 



"^[mlJna] (gi, g2, •••, qm\0l^\0){0\0n2\qu ^2, • • • , qm) ■ (97) 



^^ By RPI, ni and n2 do not have to be exactly equal, but must concur up to an angle of O(A') in SCETi. 
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Amplitude squared for the LO SCET shower operator Cj^ lo^n 



FIG. 13 

drawing the less intuitive squared amplitude in SCETtv 

SCETi Feynman diagram in order to emphasize the simple ladder structure. 



. Rather than 
we illustrate the process here with a cut 



This relation between ni and n2 is simple when the difference is encoded in the collinear 
fields in operators. However, as discussed in Sec. lIVAt we also have to deal with situations 
where this information ended up in Wilson coefficients when matching SCETj to SCETj+i. 
It is to guarantee a relation like Eq. ( 1971) that our Wilson coefficients contain G-functions 
(c/. Eqs. [5^ and [D21|) . which will cutoff the overlap regions in phase space once we begin 
integrating. The amplitude squared is particularly simple in SCET^r, where we have only the 
operator CJ^ (^at, n'^, . . . , n'j^_^), and where each particle is defined in a different collinear 
direction. 

SCETjv (or SCETj, in general) easily distinguishes which configurations are strongly- 
ordered by the structure of their Wilson coefficients. This means that we have no interference 

and Cj^ NLo ^n where Cj^ lo i^ ^^^ ^O SCET at coefficient given 



between Cf "o^^"'^ 



and C[t Miri is in Eq. (IE12p . Even though the C's are the same, the 9-functions 



-'TV.NLO 



in Eq. (|63 

in the C"s enforce different conditions, where the former is strongly ordered, while the latter 

is not. Thus, in the analog of Eq. ( 1971) . the Kronecker delta will give zero. 

We get a further simplification when we square the NLO(A) contributions. Looking at 



Cj^ NLO i^ detail, we have: 



C 



(iV-l)J 
iV,NLO 



N-2 

EC. 

1=1 



{N-1)J 
Af,NLO 



(0 



(98) 
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where 



i-i 



CSTnLoW =Y1 [(n^LL ^\f^k-l, /ifc)cLoK-l))^LL^^(^fe-b /ife) ® /^/'''(nz+i, nj, nj+i) 



7=1 k=l 

N-1 



X ( n Ut'\f^k-u /i,)c^M^^-i))]r"- (99) 

k=l+l 

In C]y j^Lo (O5 'we have made exphcit that the /, (/ + 1)**^ gluons violate strong-ordering and 
come with the factor h"'^ of the subleading sphtting rule, Eq. (IHHj) . The sum in the last 
term over / counts the different types of NLO jet-structure terms given in Eq. (187|) . The c'^q 
are defined in Eqs. fl63|) and (IE7p . and the f/'s are running factors given in Eqs. (IT8|) - (!20|) . 
The complete explanation of the symbols in Eq. f l99p can be found in the discussion around 
Eq. flElSp . The convolution factor is explained below Eq. f l89|) . Since different / correspond 
to a violation of strong-ordering at different points in the shower, each of the C^j^h^q (^ ^ + 1) 
encodes a different G structure. Therefore, there is no interference for different values of I, 
and we have that the amplitude squared to NLO (A) for jet-structure corrections (we call 
corrections of C(A^) at the amphtude squared level NLO(A)) is just the sum of squares of 
the individual operators: 

I^^^^-^^^^"^ILnlo = l^^^^^-^^^iL + l^^^^-^^^^"'lLo , (100) 

where 

i^.(A^-iM|2^=ic7f-)n(o|or^)|g(iv-i)M)r, (loi) 

N~2 

1=1 
and \q{N — l)gq) indicates the state with A^ — 1 gluon emission. The simplification even 
extends inside each of the terms, since the j*^ gluon only gets contracted with itself. Dia- 
grammatically, this means there are zero nearest-neighbor crossings in the |LOp diagram, as 
we see in Fig. [131 ^-nd a maximum of one in the |NLOp, Fig. [TH We thus only slightly modify 
the factorized emission formula, Eq. (|2]). Even for an arbitrary number of gluon emissions, 
we at most have to take into account a single defect that involves a full two-particle phase 
space. 

We can see why terms that have non-trivial interference with more than two gluons are 
suppressed by looking at the propagators in the amplitude. The amplitude for i + 1 emissions 
has a factor l/qf x l/g| x ■ ■ ■ x 1/qf. The LO term comes form the strong-ordered region where 
I1 ^ I2 ^ ' ' ' ^ lh Eq. (II]). The jet-structure NLO(A) is given when gj ~ '?i+i' "which 
allows the two gluons kj+i and kj+2 to share the same region of the phase space and therefore 
interfere. To have an overlap of three or more gluons, we would need g| ~ g|_,_]^ ~ . . . ~ Q'j+fc, 
which is clearly suppressed beyond NLO (A). 

2. Interference for Hard- Scattering Corrections 

The corrections to the differential cross section to O(A^) involve squaring the subleading 
hard-scattering amplitudes as well. Unlike the jet-structure case, these involve amplitude 
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FIG. 14: Contribution to the amplitude squared of the jet-structure piece at NLO. We show a cut 
SCETi Feynman diagram to emphasize that the square of the SCETtv operator, Cj^ j^^q 0]y , 
contains only a single deviation from the simple ladder structure appearing at LO in Fig. [131 



terms up to NNLO(A). As we argued above, they only modify the gluons closest to the hard 
interaction. Thus, we will not need to sum over many terms as we do in Eq. (1101 p . In fact, for 
hard-scattering corrections, we only need to worry about interfering SCETj operators that 
arise from acting with the LO replacement rule Eq. ( 1561) on either C| Iq, C[ ^lo' andCJ ^, 
given in Eqs. fj42l) and (1771) . or 6*2 nnlo given in Eq. (ID65p . Since the 3'''^ through i^^ gluons 
arise from the LO rule for all three coefficients, they proceed as in the |LOp case. The 
interference to look at in detail is that of the first two gluons. In SCETjv, we have: 



,Aq{N-l)gqH,2 

1^ Inlo 



(I 

+ 



MN-1)-\MN-1)H 
^N, LO ^N, NLO 



' "-"Af.NLO "-"Af.LO I 



+ 



,^{N-l)j^(N-l)H 
l'-"Af.LO "-"Af.NNLO 



+ c 



N, NNLO *-" 



1) 
Af,LO 



["-^Af.NLO I 
l(JV-l) 



\)mO\^-'^\qiN 



(102) 



The Wilson coefficients are found in Eqs. (IE6p . (lESp . and (IE9p . respectively. Nontrivial 
interference in Eq. ( I102p occurs between the ffist two-gluon emissions. 

The interference between LO and NLO(A) simplifies in many cases of interest. For ex- 
ample for one-gluon emission. 



l^S 



lo/nloIa^j 



MlPq 



ki±u{nf, -rif,). 



(103) 



If we can cleanly separate the initial and final states {e.g. e+e^ — )• jets), then by a classic 
proof involving the Ward identity (reproduced, for example, in [76|), once we have integrated 
over final state vector quantities (we can keep scalars such as Zi unintegrated) , the resulting 
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differential observable depends on 5''^'^|Alo/nlo|^i/; which for Eq. fllOSp is zero. This is 
quite straightforward for leptonic initial states, and one may be able to extend it to certain 
hadronic ones as well. 

One can account for these corrections by modifying the hard-scale matrix element and 
then running a parton shower modified to include the different no-branching probabilities 
for different phase space configurations of the same particle content. In the next section we 
discuss using a reweighting to implement these corrections. 

E. Correction Summary at Subleading Order 

In general, our corrections avoid double counting issues, because all contributions, 
whether LO, hard-scattering, or jet-structure corrections are kept separately with distinct 
G structures. Given the SCETat amplitude for A^ + 1 final state particles with corrections 
implemented both for the branching and for the no-branching, one can consider reweighting 
a LL shower in order to implement our results. For correcting the abelian emissions off a 
single quark line, this weight factor would take the following form: 

[J(iV-l,0) + g(iV-l,0)] 

A{N- 1,0) ' ^ ' 

where A{N — 1,0) is the LL amplitude squared for A^ — 1 emissions from the quark line, 
J{N — 1, 0) includes the LL result along with power corrections and subleading resummation 
associated with jet-structure corrections, and H contains hard scattering corrections. With 
our LO(A) result, 

A(iV-l,0) = |<Lo)n(0|O{f-^)|g(iV-l)M)r, (105) 

but in general A could be whatever amplitude squared a particular shower algorithm has for a 
given configuration. Eq. f ll04p then reweights that particular shower to our NLO(A) corrected 
result. An example of shower Monte Carlo with an analytic expression for A{N — 1, 0) is 



GenEvA [54|. For a leading log shower without an explicit formula for A(N — 1,0) one can 
use Eq. f llOSp with the understanding that it is likely a good approximation to the shower 
output. For the terms in the numerator of Eq. fll04p we have: 

j{N- 1,0) = \d^''^^;^^'j^mo^^^^^^^^ 

H{N -1,0) = \Ci,''-l^o''\'mO^N~'^\q{N-l)gq)\' 

(106) 

where we give formulas for C}^ lo ^^^ ^n ~ ^^ ^^- dSS]) , Cj^ ^-^{^ ^ is discussed near 
Eqs. dHD and ([92D, C^^nlo is given in Eqs. ([98]) and (199]), cJ^^nlc? is given in Eq. (|E8]) . 
and Cj^ NNLO is given in Eq. ( IE9p . Our operators, O]^ , describe a process with A^ — 1 
emissions off the quark line. The A{N — 1, 0) amplitude squared in Eq. (llOSp is contained 
within the first term in J{N — 1,0). As discussed in Sec. IIVC| while we have worked out 
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the real emission terms (C^ j^^lo ) completely, we have yet to determine the subleading RG 
kernels needed for C}^ -^jj^ ^. 

We introduce the A, J, H notation to describe more general abelian processes. y4(j, k) 
gives the amplitude squared necessary for the LL shower of j gluons collinear to the quark, 
and k collinear to the antiquark. The correction, J(j, fc), contains the virtual and real 
corrections necessary for NLLexp resummation of collinear logs. Since it contains an implicit 
sum over insertions of a single defect, which can occur anywhere in the shower, it depends 
on the total number of collinear emissions. Including the hard-scattering contributions to 
NLO(A) only requires modification of the first two emissions, after which one simply uses 
the LO replacement rule, Eq. ( I56l) . In the general case we denote it by H' , which differs 
from the above by including corrections to antiquark emissions as well. These are easily 
obtained by charge conjugation. Thus, an abelian two-jet process with j + k gluons gets the 
following reweighting factor: 

wU. k) = [J{j, k) + H'U, k)] /A{j, k). (107) 

These weight factors are positive definite. All contributing terms are squares of amplitudes, 
except for LO(A) x NNLO(A) in H{N — 1, 0). This contributes in the same region of phase 
space as the LO(A) amplitude squared, and the sum of these terms is positive. In the full 
nonabelian case, with the presence of gluon splittings, one must sum over possible shower 
histories in writing down the analog of Eq. fll07p . Algorithms for handling this complication 



can be found in hj, 54 



In Table [TTl we list concepts that are addressed by our shower framework at subleading 
order, and associate these concepts with corresponding calculations in SCETj. This table 
provides a summary of our results which appear in the weights given in Eq. f ll07p . as well 
as pointers for future calculations. Since it is easier, in the table we use the language of 
SCETi and SCET2 to discuss the corrections, rather than referring to terms in the final 
SCETat. In SCETat, the features of the SCETi operators that avoid double counting and 
allow the various contributions to be distinguished are encoded by B functions in the Wilson 
coefficients, and the operator language makes the discussion easier. For the total differential 
cross section, we found at NLO(A) two kinds of power corrections. This includes a set 
of matrix-element corrections called hard-scattering corrections (Sec. IIVAI) . and a set of 



contributions that improve double real emissions that we called jet-structure corrections 
rSec HVB]) . 



In the the hard-scattering category, we have overall three different kinds of corrections. 
The first is due the the SCETi operator Xni^n'Xn that gives the SCET2 coefficient Cg ^lo 



in Eq. flHOj) . This is an improvement of the hard matrix element that takes into account the 
emission of an extra parton at the hard scale. The second is due to the SCETi operators 
Xno^noXn and Xno['P±'B„o]xn that give the SCET2 coefficients C^2,mlo and C^^-^^i^o in Eq. (EQ]). 
This correction also accounts for more partons, but it describes a situation where they are 
initially emitted close to the collinear quark. Therefore, they are corrections which improve 
the description of the first branching within the leading jet. It is important to note that 
because these two types of hard corrections occur in different regions of phase space they 
have different renormalization group evolution, and thus different Sudakov no-branching 
factors. The required LL Sudakov factors were determined in our analysis. For a full 
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Category 


Shower Concepts 


Quantity in SCET^ 


Found In: 


Hard Scattering 


Hard matrix elements 
with more partons 


Wilson coeff. of Xm-S^^^Xn 
in SCETi 


Eq. pil) 


Power correction to initial branching 
within the leading jet 


Wilson coeff. of Xno-B^Xn 
and Xno[P±S^o^Xn 


Eq. m 


0{as) hard virtual 
correction 


One-loop matching for 

XniXn 


See [24, 25j 


Jet Structure 


1 — > 3 Splitting functions 


Double gluon real 
emission in SCETi 


Eq. (|871) 


Combining 1—7-2 splittings with 
the various 1 — )• 3 splittings 


Compute weights from SCET 
squared amplitudes 


Eq. (I107P 


0{as) virtual correction 
for LO 1 ^> 2 splitting 


One-loop correction 
to 1 — 7> 2 replacement rule 


Left for 
future work 


No Branching 
Probabilities 


NLL Sudakov factor for 
leading branching 


NLL anomalous dimension 
for leading operators 


SeeSecUVCn 

Left for 

future work 


LL Sudakovs for 
subleading branching 


LL anomalous dimensions for 
subleading operators 


Eqs. (I20I90P 


Soft Emission 


Subleading corrections 
from soft gluons 


Include effects of soft 

emission from subleading 

SCET soft Lagrangians 


Left for 
future work 



TABLE II: Mapping between concepts in an NLO parton shower algorithm and computations in 
SCETj. For exclusive cross sections these ingredients would together yield results accurate to NLO 
in the power expansion (A), and with corresponding NLL resummation. 

NLO(q;s) treatment we also need a third type of hard scattering correction, the one-loop 
virtual corrections to the leading shower operator. For the required operator, XnoXn these 
types of corrections were discussed in Refs. 2J, |25 . 

For the jet-structure corrections, there are several ingredients to consider. We derived 
a replacement rule for two emissions 1 — )■ 3, Eq. ( 188|) . that involved three different types 
of terms. This correction takes into account emissions in a region of the phase space that 
is not strongly-ordered and automatically avoids double counting from multiple 1 — )■ 2 
emissions. ^^ In addition at NLO(as)/NLL we require the 0{as) virtual correction to the 
LO splitting rule. This would be derived from a one loop matching computation that should 



^^ The method by which we avoid double counting for two gluon emission should be obvious, coming directly 
from our implementation of the Q functions. Since + = 1, the double 1 — >■ 2 and 1 ^> 3 together cover 
all of phase space without double counting. For three emissions we have either i) three 1^-2 emissions, 
ii) a 1 — > 3 followed by a 1 —> 2 emission, or iii) a 1 — ;■ 2 followed by a 1 ^' 3 emission. Here there is an 
apparent combinatoric issue, as ii) and iii) both provide corrections for the middle gluon in i). However 
they do so in nonoverlapping regions of phase space. The same is true for more than three emissions. We 
thank J. Thaler for asking this question. 
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be straightforward, but was not considered here. 

We also discussed how no-branching Sudakov factors are associated with the operator RG 
kernels, and by extension their anomalous dimensions (Section lIVCp . To NLL, we need the 
NLL Sudakov factor for leading branching and the LL Sudakovs for subleading branching. 
These are associated to the full one-loop and two-loop cusp anomalous dimensions for the 
leading operators, and one-loop cusp for the subleading ones. At LL, we have determined 
all the Sudakovs for subleading branching (Eqs. [20] and [20]) • We have not yet calculated the 
NLL Sudakov for leading branching in the scheme with B-functions that is needed for our 
setup, as described in Sec. ilVCi 

The last item in the table is the treatment of soft radiation at NLO. This can be achieved 
by considering time-ordered products for the matching of QCD to SCETi and SCETj 
to SCETi+jthat involve soft gluons and subleading soft Lagrangians that are known in 



SCET [6J-[66| up to 0{X ). For the terms involving coUinear quarks they read 



4^ = i^nW) ^Ipts^{W^^Iftlin) + {L^JIiW)l^Ilt.{W^lin) , (108) 

4? = ^"« -^ ^gKWqus + h.c. , Cf^ = ^tPl-J—- tg ^« W qus + h.c. , 

4r^ = ^nTT •- n [^^ ■ ^ri, in ■ D^ + Qu ■ Aus]W q^s + h.c. , 

while the analogous pure glue Lagrangians can be found in Ref. |6J]. Here the expressions 
are prior to the soft field redefinition, and ig^]^ = [in-Dn, ip±\- One must then work out the 
effect that these NLO soft amplitudes have on interference. The associated soft calculations 
and investigations have also been left for future work. 

We also briefly comment on how the corrections in Table [TT] relate to those already im- 
plemented in parton shower codes in the literature. In most cases, the goal of these codes 
differed from the power suppressed corrections considered here. This makes a strict associa- 
tion impossible, but there is still a general correspondence that can be made. CKKW [27| is 
a LO(as)/LL procedure whose goal is to merge matrix elements involving multiple partons 
with a parton shower in a manner that avoids double counting. In our language, this corre- 
sponds to the real emission hard- scattering corrections in the first row of Table [IT] The XnoXn 
and Xni-B^^/Xn operators describe processes with different numbers of initial well-separated 
jets. In CKKW, a parameter ?/cut is used to separate the extra emission in the matrix ele- 
ment from emissions in the shower. In our analysis, the contributions from showering XnoXn 
does not interfere with the direct contribution from Xni-B^J^/Xn, and this is encoded by 9 
functions in the Wilson coefficient of SCET^r. CKKW carries out this procedure for several 
matrix element emissions, while we have only considered one. 



In MC@NLO [37[ and POWHEG [38], virtual and real matrix element corrections at 
NLO (as) are incorporated into the shower, with the goal of ensuring that it reproduces 
an associated cross section completely at NLO(as). The implementation includes careful 
handling of the cancellation of real and virtual IR divergences. Our goal was to implement 
corrections at NLO (A) and we discussed NLL, but for all emissions from the shower rather 
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than just the first jet needed for the NLO(as) cross section. At NLL, we would have 
only terms up to 0{as log) in the total cross section, and hence this does not encode the 
entire NLO(as) result. In our language, the corrections that contribute to the NLO(q;s) 
cross section correspond to the hard scattering corrections in the first through third rows of 
Table HTl In order to compute the NLO(q;s) cross section it is not necessary to distinguish 
between the terms in the first and second rows of the table, and these terms are indeed 
considered simultaneously in MC@NLO and POWHEG. The full NLO(as) virtual result 
are obtained in our language by including the items mentioned in the 3'''^ and 9*^ rows of 
Table HI 



The work of KRKMC 57H59| , on the other side, aims to improve the shower algorithm 



taking into account an exclusive version of the Altarelli-Parisi splitting function at NLO(q;s), 
Pqq ■ In our language, this corresponds to jet-structure corrections and we show in Appendix 
|F]how our replacement rule in Eq. (IHHj) is also related to Pqq ■ Hence our 1 — )• 3 emission 
corresponds to an exclusive version of Pqq , though in a different scheme. Part of the 
corrections in Pqq involve order a^ corrections to the 1 — )■ 2 splitting function, which are 
taken into account by 0{as) virtual 1 — )■ 2 matching corrections in our framework (6*"^ row 
of Table HTl) . In fact, in Sec. IIVDI we saw that SCETj also leads one to view corrections to 
the shower as a "defect" insertion just as KRKMC. In addition to these splitting corrections, 
in our framework the amplitude also involves no-branching probabilities given by evolution 
kernels that appear in the weight factors, which do not appear in the KRKMC weights. 
Keeping track of the evolution also determines the appropriate scale for evaluating as- 

V. CONCLUSION 

In this paper we developed a framework based on a tower of independent but related 



EFTs, the SCETj, to study corrections to the parton shower. The work of \24\ . |25| showed 
how to formulate the LL parton shower in terms of SCET, and how virtual corrections are 
straightforward to incorporate by one-loop matching. Our SCETj framework extends these 
ideas in a manner that makes it easy to deal with: double counting, the issue of disentangling 
coordinate choices from kinematic power corrections, and the construction of a complete set 
of operators for corrections at a desired order. The interference structures, and hence the 
leading corrections that give spin and color correlations, also appear in a straightforward 
manner in the SCETj setup. 

The SCETj are iteratively used to integrate out the characteristic scale, QA* for increasing 
i. This approach allows us to perform a systematic expansion which can correct both the 
hard-scale process that produces partons to setup initial conditions for the shower algorithm 
and the iterative shower itself. We described the parton shower through a set of operators 
Ol in SCETj, and used standard matching procedures to make the transition from SCETj 
to SCETj+i, where more partons become apparent. Performing the matching relied crucially 
on the RPI symmetry of SCET, and we extended the usual infinitesimal version to carry out 
the finite rotations that we needed. At LO, a simple operator replacement rule generates the 
LL shower, Xno ~^ ^toXniQ^n' l^ where clo is related to the standard LO sphtting-function. 
Also, angular ordering and coherent branching for LO soft emissions emerge naturally in the 
SCETj framework. A summary of ingredients required for the shower with power corrections 
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at NLO(A) are given in Table [TTl including both calculations carried out here, as well as those 
left for future work. The main results of our paper are: 

1. At NLO(A) we found two kinds of branching corrections: hard-scattering and jet- 
structure. The hard-scattering corrections depend on the hard process and appear 
near the top of the shower tree. They came from matching QCD to SCETi at higher 
order. Since they only occur at the top of the shower, one can treat these as a modified 
form of matrix-element corrections. A subset of these corrections correspond to the 
usual implementation of fixed-order matrix elements, while the remaining ones give 
power corrections to the initial branching in the LL shower. These two types require 
different Sudakov factors. This effect is apparent for the kinematic power corrections, 
but is beyond NLO(q;s) for the fixed order counting. 

2. The jet-structure corrections are independent from what happens at the hard scale, 
hence they are universal for any process we want to study. They come from matching 
SCETj to SCETj+i at higher order for any i. They can appear anywhere in the 
shower tree and they take into account emissions in regions of the phase space that 
are not strongly-ordered. For these corrections we found that the NLO(A) operators 
are related to the LO operator via a replacement rule for two emissions: Xm ~^ 
^iXn29^n' ±9^n' ±- '^^^^ NLO(A) rule automatically avoids double counting with the 
iteration of two LO operator replacements. 

3. The SCETj picture allowed us to easily take into account interference for the NLO(A) 
power corrections. Once we reach the final SCETat theory, all the fields are labeled in 
a different collinear directions. Because in SCET we can only contract collinear fields 
that share the same collinear direction, in SCETjv calculating the amplitude squared 
becomes very easy. Kinematic information that is encoded by the shower history from 
passing through earlier SCETj's is encoded by 6 functions in the final SCETtv Wilson 
coefficients. We demonstrated that when emitting an arbitrary number of partons, 
the non-trivial part of the amplitude squared involves at most four fields. 

A comparison of how these SCETj results relate to earlier parton shower literature that goes 
beyond LL is given in Sec. IIVEI 

The framework developed here allows for systematic improvement to arbitrary orders in 
the kinematic expansion. There are still several important steps to take, though, before this 
picture can lead to a practical implementation, including additional computations that we 
outlined in Sec. IIVEI We list here three topics which are natural next steps, and which we 
believe should be straightforward to approach: 

1. This work has only considered q ~^ qg splittings and an abelian theory. One should 
include the full nonabelian results and compute the coefficients required for gluon 
splitting as well. This is required to properly treat color correlation corrections in a 
manner determined by the NLO(A) interference pattern. For collinear particles we 
expect that one can include the dominant part of these effects by considering nearest- 
neighbor interference since this arises from the kinematic expansion, and thus leaves 
the rest of the shower as before. 
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2. Only a subset of the terms required for a full NLLexp resununation were considered 
here. We determined the LLgxp evolution for subleading operators, but did not carry 
out the computation of the NLLgxp evolution of the leading operator in a scheme that 
is consistent with our power corrections (we only considered it in MS). In order for 
a consistent treatment as a probabilistic process, the real emission probabilities and 
Sudakov no-branching corrections must go hand in hand. Furthermore, once these 
evolution factors are determined, the reweighting discussed in Sec. IIV El must be tested 
in an actual shower Monte Carlo. 

3. Since soft modes in SCET can communicate between different collinear jets, they carry 
the ability to spoil their factorization. Fortunately, this does not happen for their LO 
interactions, which yield angular ordering and coherent branching of soft gluons in 
SCETj. It is open question as to what extent NLO soft couplings can be factorized in 
the shower tree and the necessary SCET computations were discussed but not carried 
out here. The treatment of soft NLO interactions in SCET in other contexts has 
always led to factorized structures, so we remain optimistic that such effects will be 
tractable for the shower. 

Future investigation of these items is well warranted. 
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Appendix A: More SCET basics 

Soft-CoUinear Effective Theory describes the interactions of collinear and soft quarks and 



gluons 20|-|23| . As we mentioned in Sec. Ill At to define the collinearity of a particle, the 
momentum is decomposed along two light-cone vectors, n and fi, with ra^ = 0, ra^ = and 
n ■ n = 2 

V^ = n-p—+p—+p^^, (Al) 

where p = n ■ p. A particle is collinear to the direction n if its momentum scales as: 

(n-p,p,px)~(A2, 1, A)g, (A2) 

where Q is the hard scale of the process, and A ^ 1. A particle is soft if: 

(n■p,p,p^)~(A^A^A2)g. (A3) 
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We obtain SCET from QCD by expanding in powers of A and integrating out modes harder 
than ~ Q'^X'^. Both Eqs. (jMl) and (jMl) imply that p"^ = p{n ■ p) + pi < Q^A^ 

In addition to the expansion, we also want to divide the quark and gluon fields into 
separate soft and coUinear modes. For the collinear case, the fields are indexed by n, and 
two collinear sectors are distinct iini-rij ^ A^. In addition, we introduce a momentum-space 
lattice for the (9(A°) and 0{X) momenta in order to facilitate carrying out the multipole 
expansion with respect to the O(A^) momenta. To divide the QCD fields in this way, we 
split the momentum of a collinear particle into a "large" part p'^ and a residual one A;^ ~ A^ 

p^ = p^ + k^ , where p^ = n-p h p^ . (A4) 

We can pull out the large momenta p from the fermion field by the phase redefinition 

^{x) = Y,e-'^'^^n,p. (A5) 

For a collinear particle along n, d'^ipn,p{x) ~ A^. The four component field, ipn,p, has two 
large components, ^n,p, and two small components ^n,p, that can be separated using the 
following projectors: 

These satisfy the relations, 

~j~ Sn,p ^«,i5 5 r' <,n,p U . \-^' ) 



Similarly, we can define a collinear gluon field, A'^Jx). Pictorially, we can think of C,n 



X 



and A'^ Jx) as fields that create a particle whose three-momentum lies inside a cone with 
opening angle ~ A about the three-direction n. V!^ is the momentum operator that picks 
up the large components of the momentum, 'Pf^C,n,p{x) = p^ ^n,p{.x). Collinear fields always 
appear with a sum over p, and both label and residual momenta are separately conserved. 
Therefore it is often useful to abbreviate the notation as 



Sn / ^ Sn,i5 y ^n / ^ ^n,q ■ yi\C 



The SCET collinear Lagrangian, £„, describes the interaction between the collinear fields 
$,n and A'^{x). It is derived from the QCD Lagrangian by integrating out the field, ^n- At 



LO, for the kinetic and purely collinear interaction terms we have [2l|, |22 



C^^^=^^n■^ + gn■A^ + ^Ip^^W^^Wl^Ip^^)^^r., (A9) 

' n ^ 
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where we intrinsically sum over the large, label momenta, p. The in ■ d derivative picks out 
the O(X^) momenta. The collinear derivative, D!^, and collinear Wilson line, Wn, are defined 
as |23| : 

^D^ = Vii + gA^^ , 
Wn{x) = \ J2 exp( - J-n-A„(x))j . (AlO) 

perms. 

The leading order coupling of collinear quarks to soft gluons is eikonal, 

4°)=a^n-^|en, (All) 

while the Lagrangian for purely soft quarks and gluons has the same form as full QCD. 



The LO collinear Lagrangian for gluons has similar properties and is given in Ref. 23 



The interactions between soft and collinear particles, such as the one in Eq. (lAlip . can be 



removed from the Lagrangian by the field redefinitions 23 



Sn -^i^n^n-, ^n ~^ ^n ^n ^n ! (^12) 

where the soft Wilson line Yn is defined in Eq. (!67|) . This causes soft interactions to be 
represented by Wilson lines in operators, as in Eq. ( 169|) . 

Now that we have split up gluons according to a momentum-space lattice, the gauge 
structure of the theory has become more complex and involves global, collinear, and soft 
gauge transformations. Fortunately, with the collinear Wilson line, it is possible to construct 
fermion and gluon fields that are manifestly invariant under collinear gauge transformations. 
The definitions are: 

Xn{x) = Wl{x)U^) , S^(x) = - [Wlix) tD^Xx)WM] , (A13) 

9 

where the derivative in i3^ does not act outside of the brackets in its definition, and we always 
have n-Bn = 0. In the n-An = light-cone gauge, Wn = 1 and Bl^ = A'^. One can construct 
collinear operators out of just three objects: the fermion field, Xn-, the perpendicular gluon 
field, -B!^! , and the perpendicular momentum operator, "P^i . All the other operators, like 



n- Bn, or n-d can be written in terms of these three using the equation of motions [67 . 



Appendix B: Finite RPI 

Even though SCET explicitly breaks Lorentz invariance, the symmetry returns at each 
order in A by reparametrization invariance (RPI). RPIj is the version appropriate for SCETj. 
As usual, we define p as collinear to the direction n in SCETj if its components scale as 
{n-p, p, p±) ~ (A^*, 1, A*) Q , where Q is the hard scale and A ^ 1 (c/. Eq. IA1|) . The 
vector n has physical meaning as its 3- vector subset, n , is the direction where most of the 
momentum is allocated. The direction p is therefore inside a cone of opening angle A* 
around n , (cf. Fig fTSl) . By contrast, n is an auxiliary vector only needed to decompose the 
momentum. The parameter A gives the amount of coUinearity to n. The decomposition is 
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FIG. 15: In SCET, a particle is collinear to the direction n if it is inside a cone centered in n and 
of opening angle A. 



not unique since we can shift n by an amount A and the particle will still be collinear to it. 
This means that if we move n inside the cone in Fig. [151 P is still collinear to it. This is 
called a reparametrization invariance (RPI) transformation of type-I. Thus, if a particle is 
collinear to n, it is also collinear to any direction n' related by a type-I transformation. To 
be more formal, we can divide the space of light-cone vectors, {ui}, into equivalence classes, 



{[rij]}, where 



rii 



{n e [rij 




< A^'}. The meaningful objects in SCETj are the [uj 



By extension from SCET [63|, two collinear sectors in SCETj, rii and n2, are distinct if 



ni-n2 > A 



2i 



(Bl) 



Just as in regular SCET, we can write the external state with the n-label to which each 
particle is collinear. For working in SCETj, we give a subscript to indicate the appropriate 
definition of collinearity. For example, \qni)i is a state with one quark, collinear to ni that 
can be annihilated by any Xn such that ra ■ ni <C A^*, or n E [rii]. 



For each {n,n}, the type-I RPI infinitesimal transformations are 



13 



(I) 



n'^ 



n^ + Ai: , 



n^ — )■ n^ 



(B2) 



where A^j^ ~ A and n ■ A„x = ^ ■ ^n± = 0. These transformations preserve the relations 
n^ = 0, n^ = and n- n = 2. -^^ 

The general problem of matching SCETj — )■ SCETj+i is our need to rotate the direction 
n of objects in the amplitude (such as spinors and vectors) to n' that is close enough to 
the particle momentum such that p is collinear to n' in SCETj+i. Thus, RPIj is crucial for 
matching as it determines how formerly identical SCETj configurations wind up in different 
SCETj+i terms. Any transformation in RPIj/RPIj+i is therefore of consequence. By con- 
trast, the choice within SCETj+i is purely a convention we may use to our convenience (c/. 
Fig. H]). For example, we can pick n' as that direction np such that p has zero perpendicular 
momentum in the Up — n frame: 



nt. 



n'^ 



p 



p—^ + n^-p — 
^2 ^2 



(B3) 



^•^ Infinitesimal does not refer to the expansion in A. 
^"^ It is also possible to rotate n ^)- n + e± where £± ^ 
type-Ill transformation takes n — >■ e"n and n —> e~ 



A", which is a type-II RPI transformation. Finally, a 
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This is satisfied for: 

< = n^ + 2 ^ - n^^^ , (B4) 

^ p p'' 

witfi p± defined in tlie n-frame. Unlike Eq. (1B2I) . this RPIj transformation is finite. It is 
easy to check that rip = 0, Up-n = 2 and that Pn±=P^~ np-pn^^/2 — pn!^/2 = 0. 

We can derive similar relations for other quantities. To see how the quark field transforms, 



we use the RPI invariant fermion field 67 



^n = (i + ^^)e.. (B5) 

V n ■ Un 2 / 
Since (IBSP is invariant under RPI, ipn = ipup and we can write. 

Multiplying ( JB6I) by the projector ^^/4 we get the finite RPIj relation 

a = ff„.. (B7) 



The relation flB7P is in agreement with the spinor equation (AT) in 2J] upon setting ni = n2. 
Objects with a full Lorentz index, like p^ or 7^^, are RPI invariant as there is no reference to 
the light-cone vectors n and n. Those in the perpendicular direction though, such as p^ or 7^, 
are not, as _L is defined with respect to n and n. Using the relation 7^ = 'y'^—n^ Tfi/2—n^ fl/2, 
we derive the expression 

7^^ = l1-n^^-^-pl^-+ n^^^i ■ (B8) 

^ p p p^ 

We now focus on those transformations needed for one-gluon emission. As in Sec. IIIIAl 
we consider the case of a virtual quark with momentum go emitting an external gluon and 
quark with momentum ki and gi, respectively. In Fig. [T6r A). we portray this kinematics for 
one-gluon emission where the initial quark go comes from a QCD current g7'^g. We call no, 
n'^ and ni the directions where go, fci and gi zero have perpendicular component, that is: 

no n^ 

In, 

— n 
ki = h^, 

gi = gi^, (B9) 

Using Eq. ( 1B4I) . we can relate n'^ and rii to no, 

n^ = n[; + 2^^^i^ - n'^-^^i^ , (BIO) 
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go,^o, 



qi,ni 




9i,m 



A . 

92,^2 /C2,^ ki,n[ qi,m k^,n'^ ■■■ ki,n[ 
[B] (C) 



FIG. 16: Kinematic variables for one-gluon emission (A), two-gluon emission (B) and i-gluon 
emission (C). The n's are defined such that the corresponding particle's momentum has no per- 
pendicular component along the directions n — n. Note that qi is a different vector for single and 
double emissions. 



where we have used the equahty {ki 



where 



'no- 



-{q 



lJno± 



. Some useful relations are: 



<& 



kf q( qi kf 



?1 



Un 



kiTi^ + qin'i 



— n^ (ni-rii) 



i^Qih 



nr 



^2\ 



(9l)no- 



qih , jv^ _ , ,.qiki{kf-qf 

qo ^ 



^3 



4gg 



7 



"0- 



7 



n 



M^-^i 



I no- 



ql 



ki 






u //I 

?2'[ — n-i 



\/ni ■ n[ 



and \vf\ = 2. Another useful relation is 



go' = (g, + k,r "' 



n^-rii- 



qiki 



(Bll) 



(B12) 



(B13) 



We can express all quantities of interest in terms of the vectors n[, rii and the momenta gi 
and ki. 

In two-gluon emissions, the kinematic variables are assigned in Fig. [T6](B). We define uq, 
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n']^, n'2 and rii as follows (note that gi and ni are different from above): 



rin n^^ 



% = (io^ + nQ-qo—, 



7, "o 
I 

— 71 

Qi = qi^ + ni-qi — 

la 

V. T "-2 
k2 = k2— , 

A* 

- «2 



2 
Eq. (jBlOp is still valid, and we can similarly define n'2 and n2 as: 



(B14) 



A;2 kl 



77-2 =n'l — 2 — =-^^ n^ 



712 = ni+ 2 — ^-^ ""- — =2"^ ' (B15) 

?2 q'2 



where {k2)n^_i. = — (Q'2)ni±- Also, Eq. (IBlip is still valid, and we get a new set by sending 
^ 1 and 1 ^ 2: 

ifi fi o (92)na9i .„^„. 

77,2-^2 = rii- 77-2 r^ = rii- 77-2— = -2 ^^^^ , (B16) 

"^2 ^2 ?2 "^2 

k2n'^ + q2n'^ _. . ,.^2^2 
< = T n^" {n2-n2)—2 , 



( ^/. 92^2 / T^2 -ill I^ (l2k2{kl-ql) 

9i ^ ^1\ 



<x = 7„V - n"^ - (*):a| + "" ("-'4)|ji<. 



where 



no — TiX^ 



< = ^= , (BIT) 

a/?72 ■ n'2 

and \v\\ = 2. We can write rii ■ 7i[ and vi in terms of 77,2, n'^ and n'2 so that once again we 
only need to work with external quantities: 

, fc2(n'2-n;) + g2(n2-n;) , qih 

7li-7l.= (?72-n2)— 5-, (B18) 

?i _ % 

^ 2 k2qi n^ + 2 gig2n^ - A;2g2 (^2 ■ n'2) n^ - 2 (gO^n'/" 



Eq. (JBISP for two emissions is modified to: 



go = (?2 + fci + K2) = ni-n2^- + n2-?72-— + ni-?72^— . (B19) 
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Other useful relations are 



a _ a _ ({£)no^ {i^\ /Mill + Mll^^ 



fci k2 y \ ki 



n-#f^ 



2 ' J,2 



'no-L I \^^/?ii_ 



2 



gi = (A:2 + 92)^ = ^■^2^ , 

(A:i + g2)'=n;-n2^. (B20) 

For i gluon emissions, Fig. [161 0) . n'^ is parallel to the fc-gluon, rii parallel to the external 
quark, and rik is the light cone vector such that the k^^ virtual quark has zero perpendicular 
momentum with respect to (n^, n). To calculate nj, n[ we can iterate the formulas above up 
to i emissions. That is we can calculate rit, n[ from ni_i using Eq. (IBlOp with — )■ (i — 1), 
1 ^i. 



Appendix C: Matching QCD to SCETi 

To study the process oi q ^ qg emission, we match the QCD current, 

'^QCD=?V?, (CI) 

to SCETi operators for a final state with a quark, antiquark, and gluon. The particle 
momenta are gi for the quark, pq for the antiquark, and ki for the gluon, (c/. Fig. [T71) . We 
do the matching in the center of mass frame with 

p^ = qi+Pq + ki = {Q, 0,0,0). (C2) 

SCETi, being equivalent to the usual SCET, is formulated as an expansion in the parameter 
A. The current in Eq. (ICip matches onto an infinite series of SCETi operators. We will 
perform the matching up to NNLO(A) for one gluon emission, and focus only on the cases 
when the gluon is either coUinear to the quark or has its own direction. Obtaining the result 
for gluon-antiquark collinearity from our work is a simple exercise in charge conjugation. 
We can construct the SCETi operators out of a few building blocks: the quark field Xn, the 
gluon field i3"j^ and the perpendicular momentum operator V^j_, plus Dirac structures. Xn, 
B^j_ and V^j_ all scale ~ A. The basis of SCETi operators for one emission up to NNLO(A) 
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is 67 : 15 



Of\no) = XnoXn^ 



O'l\no,no) = Xn^gKoLXn 



' 71(1 ± iJ'^nn 



"0- 



Xn ; 

0\'\n,,n\) = Xn.gB':,,^^Xn, (C3) 

Following the convention of Eq. (fl^ . we do not write the antiquark direction as it is always n. 
Ol is the LO operator and scales as A^, OJ (no, uq) and OJ (tt-i, tt,'^) are NLO(A) operators, 
scaling like A^, and 7^ ~ A^. 

In SCETi, two particles are collinear if they are inside a cone with opening angle ~ A, 
equivalently pi ■ P2 ^ {Q^Y/v'^- The factor rj ~ 1/2 represents the average p fraction 
taken by the daughter from the mother as discussed in Sec. II C. Usually, we formulate this 
condition with dimensionless quantities, Up-^ ■ Up^ < A^/r/^, where n^. is exactly proportional 
to the particle momentum. To distinguish a "two-jet" from a "three-jet" state, we label 
the external states with the direction to which the particles are collinear. A state \qna)i 
indicates a state where a quark with momentum gi is collinear to the direction no, that is 
(gi,no ■ gi, (i?i)no±) ~ (1; -^^^ ^)Qj and the subscript, 1, tells us the state can be annihilated 
by any operator, Xn, where n and no are in the same SCETi equivalence class, {[?^]}. As we 
will see, when we match to lower-scale SCETj, we will change this number appropriately. 
A two-jet state with a collinear quark and gluon, and an antiquark is given by |g„Q gno qn)i- 
The fact that the quark and gluon share an index label implies that qi ■ ki < {Q\Y /tj^. 
A three-jet state is indicated by |g„i Qn'^ qn)i, where each particle is collinear to a different 
direction. The operators Ol (no), 0[ (no,no) and 7^ (no,no) can only create a two-jet 
state, whereas 0{ {ni,n[) is for three-jets. Multiplying the terms in (1C3J) by the Wilson 
coefficients, we have: 

+ C«(ni,<)0« + ..., (C4) 

where the ellipses indicate higher order terms in A. When it is unambiguous, we will only 
write the n-labels in the Wilson coefficients, as above. We begin by looking at two-jet 
operators in detail. Here, because we are in the center of mass frame, the two jets are 
back to back. We define the kinematics as follows, the antiquark is exactly parallel to 
n = (1,0,0,— 1), while the quark and the gluon are collinear to no = (1, 0, 0, 1), such that 



^^ Ti {rLi,n'i) encodes redundant information that can be obtained with RPI. For example we can choose 
the directions ni and n'l to ahgn perfectly with particle momenta such that e.g.. 'Pni±Bni± = 0. This is 
not possible for 7^ {no,no). 
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FIG. 17: Matching QCD to SCETi for the two-jet configuration: In the first column there are the 
two Feynnian graphs for one-gluon emission in QCD, labeled by the 4-momenta. In the second 
column there are the two Feynman graphs in SCETi that reproduce the same amplitude in the 
case the quark and gluon are collinear along the direction no- The first graph comes from the 
operator 0[ (no) with the insertion of the SCETi Lagrangian, the second graph comes from the 
operators 0[ (no, no) and 7^ {no, no). 

9o = 9i + kl has no component perpendicular to no and n, and: 



fi 



2 



n^ 



?i^ +nQ-qi'-^ + {qi^^ 



"o-L ' 



^ ^0 



n'^ 



where (no ■ qi,qi,qi±) and (no ■ ki,ki,ki_[_) scale as (A^,1,A), and (gi) 



-{kl 



(C5) 
, I by 



momentum conservation. The Wilson coefficients are defined through the equation 

(Ol^QCDkno^no^n)l =Ci!LoK,^o) /c?a;^(0|T{£sCETi (a;)Ci°^} l^no ^no ^n)! (C6) 

+ C'~l\no,nQ){Q\0'l\qn^gno qn)i + C'J^|(no,no)(0|T;^^^|g„o 5f„o g^)i . 

Calculating the C's for this two-jet process goes as follows. We decompose the QCD ampli- 
tude along no and n, using Eq. fICSp . and we write the QCD spinor in terms of the SCETi 
spinor, Eq. f ICSp . Expanding in A up to NNLO, on the RHS we compute the amplitudes for 
the three SCETi terms. The coefficient C| /o '^^is already determined from matching QCD 
to SCETi for zero gluon emission, it is: 



C 



(0) 
l,LO 



r 



(C7) 



The coefficients C| and C{j- come from solving Eq. (IC6p at NLO(A) and NNLO(A), re- 
spectively. Since Cj and 7^ are at different orders in A, there are no ambiguities. 
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In order to do the matching, we need the relation between the QCD and SCET spinors. 
Using Eq. (1B5I) . we can write: 



u{p) 



i>i 



2p 



Un{p) 



(CJ 



where u{p) is the QCD spinor and Un{p) is the SCETi one. It easy to see that the SCET 
spinor satisfies: 



-u„, = 



""^n '^n ) 



Yl ^nK = P^ 



_lfi 

2 



(C9) 



The QCD amplitude for 7* — )■ qqg (shown in Fig. [T7]) is: 



% Po 



(CIO) 



Using Eqs. (JCS]) & (JCl]) in ( ICTOJ) and expanding to NNLO in A we get: 



AQia _ A<iig j_ Alia _, Ama 



(Cll) 
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AZlO = 9[- + YjQU^ol'noA^l)no^{<o± - ^^n- )V. 
-9-Q-no{hrno.[lno.-^^n)v,. 



(C12) 



We already know C| £05 ^'^d it is easy to determine the other two Wilson coefficients to 
reproduce Aq^^, they are: 



C^NLo(^0,no) = g« - ^'')7°o± ' 



.(1) 



Cif(no,rao) 



-*- /i ,9 o _ 

- 7, 'noX 'no± '"o-L 



gi<5 



/"T^ox > 



(C13) 



where we have used the relation qi + ki = Q. 

For the three-jet operator 0| {ni,n[), the matching was already done in 67|, but we 



will translate it to the notation used here. In this case, we need three distinct directions in 
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SCETi to describe the three external particles, and there is no small parameter to expand 
in. This means that the amplitude for this operator is exactly equal to the tree-level QCD 
amplitude for a qqg process. One may wonder then, why we simply do not apply this 
everywhere instead of just the three-jet region. The answer has to do with factorization and 
running effects. The RG kernels of our two-jet operators, 0{\ 0{\ and 7^ , will resum 
the large coUinear logarithms of those configurations (c/. Sec. lIVCp . It is for this reason 
that we gain by keeping track of them as separate contributions. 

Even though they are all in independent directions, we need only four independent vectors 
to decompose the particles. In the center of mass frame, Qq = qi + ki is back to back with 
the antiquark, pq oc n. We decompose go along (rio, n) such that it has no component 
perpendicular to them: qQ = n ■ go^o/2 + ^o ■ go^/2. Using Eq. (IBlOp we can define rii 
and n[ such that they are parallel to gi and fci, respectively and such that the quark is 
decomposed along (ni, n) and the gluon along {n[, n), both without ± components. Unlike 
the two-jet case, where {qi)no± ^ -^j since the quark was collinear to no, here (gi)no± > A in 
Eq. flBTOJ) . We have: 

""' (C14) 



n- 



M 



— //,n • tr^- 



P^ = no-Pq—, 



1 i 2 ' 

where n- qi, Uq- pq and n ■ ki are 0{Q), and rii ■ n[ > y^ jrf . With this setup 7^ ^(ni, n'^) = 
The matching is therefore given by: 

(Ol-^QCol^ni 9n[qn)i = c[^\ni,n[){0\0'^'' {q^, gn[ qn)i , (C15) 

and the Wilson coefficient is: 



{ni-n[}qiki '^ 



+ 



[irf^ - Qi ^iTJ + 7 ^/f' T K 



n-pq)ki\ ^T / (ni-n'i)gi/c 



7^;^, (C16) 



where the subscript T applied to a generic four vector /^ means: fj. = f^— p^ (/ ■ p^)/^^, 
and p^ is defined in Eq. I \C2\} . 

Before moving on to lower scale SCETj, we note that all the Wilson coefficients in SCETi 
are of order A°. This will change with SCET2 as these factors will determine the relative 
importance of different contributions. As we discussed at the very end of Sec. IIVAI we do 



not need to compute any suppressed two-gluon operators in SCETi to the order at which 
we are working. Their Wilson coefficient will be 0{X^). Matching this contribution to a 
two-gluon SCET2 operator will leave this factor unchanged as there are no further emissions 
from it. The field content in SCET2 will scale ~ A^. As shown in Eq. (ID7p though, LO in 
SCETo is at A^ 
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Lastly, we described the effects of adding running effects in Sec. IIVCI In tlie next 
Appendix we will match SCETi to SCET2. Before doing it we have to run the SCETi 
operators from Q down to /ii, where we have the first emission: 



«o-L 

V {ni-n\)qiki ^ 






7"^ • (C17) 



V{n-pq)ki\ ^T / {ni-n\)qiki 

For the definition of the running factors U^'-'^'''\Q, fj,) see Eqs. ([l8])-([20]), and ([90]). The 
convolution symbol, ®, is only relevant beyond LL; that is beyond the level required here. 

Appendix D: Matching SCETi to SCET2 

1. One-Gluon Emission 

We now match SCETi to SCET2 for one and two-gluon emissions, starting with the 
former. The basis of SCET2 operators necessary for the matching up to NNLO(A) is equal 
to Eq. del]), but defined in SCET2: Cs^^K), C^^^(no,no), T^^^Vo, ^o),0?^(ni, n;).!^ In the 
previous section, we matched QCD to SCETi for one emission and found either a two-jet 
(kno 5'no ^n)i) or thrcc-jct configuration (l^'m 5'n'^ ^n)i), depending on the collinearity of the 
external particles. When we go to SCET2, our definition of collinearity becomes stricter. 
Particles with momenta pi and p2 are collinear only if pi ■ p2 < Q'^X'^/rj'^, where 77 ~ ^ is the 
average energy loss factor between mother and daughters discussed in Sec. Ill CI As a result 
of this change, a two-jet configuration in SCETi can be matched both onto |g„o gno qn)2 and 
llni Qn' (ln)2 in SCET2. The three-jet configuration in SCETi can, of course, only go to the 
three-jet state |g„^ (?„' qn)2 in SCET2. The matching is given by 

^^CD =Ci°[oK)0r + C;;UK,no)Or^ + C|^)(ni,n;)0« (Dl) 

= Cf (no)Of + Ci'\no,no)0^^'> + ci'\n,,n[)0^^^ (D2) 



^^ As before, we do not consider operators like O2 {no,n) that describe a gluon collinear to the antiquark. 
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where we give the decomposition into both SCETi and SCET2 operators. The eUipses 
indicate higher order terms. If we close Eq. (JD2I) with the state |g„(, gno <in)2, we get 

Ci^Lol'^o) / rfa;''(0|T{£scETi(a;)Cf^}|g„ofi'„ogn)2 

= C^°\no) / dx'^{0\T{CsCET2{x)0^2^}\qnognoQn)2 

+ C?Vo, ^o)(0|C2^^|g„o gno Qn)2 + C*^ r(^0, '^o)(0|7^^^^|g„o 9no Qn)2 ■ (D3) 



Since the structure of the operators in Eq. (ID3p is the same on the LHS and RHS, we simply 

get: 

C^'\no) = C[%o{no) , 
C^^\no,no) = Cj^J^Lo(no,no) , 
Ci]^{no,no) = C[%{no,no). (D4) 

Acting on Eq. (ID2p with the state jg^^ (?„' qn)2, we have: 

Ci°LoK) / dx'^{0\T{CsCETi{x)Of^}\qn^ gn[ qn)2 + Cj^i,LoK,^o)(0|Cf^|g„l gn'^ qn)2 

+ C[^\ni,n[){0\O\^^\qni gn[ qn)2 + Cj^r("'0>^o)(0|T/^^|g„i 5'n; 9^)2 

= C^')(ni,n;)(0|o5')|g„, gn'^ qn)2 • (D5) 

We decompose C2 {ni,n[) as 

C5^^(ni, n;) = C^%o{ni,n[) + C^]l^^o{ni,n[) + C^^J^^o(^i, n[) + C^%i^o{ni,n[) , (D6) 

where Cj lo i^ ^^^ coefficient that reproduces the first term on the LHS of Eq. fIDSp . etc. 
All the SCET2 coefficients in Eq. f lD4p scale as A°, like in SCETi, but we will see that those 
in Eq. f lD6P scale differently, giving the hierarchy indicated in the subscript. We will show 
that: 

Cf (no)Of ~ A^ , Ci]Un„n[)Oi'^ ~ A^ , (D7) 

C^'\no, no)Oi'^ ~ A^ , C^l'toi^un'M'^ ~ A^ , 

Ci%{no, n,)Ti'^ ~ A« , &^;^{^^{n,,n',)0^^^ ~ A^ , 

C^2'NNLoK,^'l)^2'^~A^ 

In the second column we have only one operator O2 {ni,n[) and we have decomposed its 
coefficient according to Eq. f lD6p . The matching does not conserve the power counting, as 
coUinear SCETi fields scale as A, but in SCET2 they go as A^. For example, we have that 
the LO operator in SCETi is CJ 0[ ~ A^, but for the LO operator in SCET2 we have 

Cf ) of ~ A^. 
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If we want to calculate a cross section for a fixed number of external particles, then we 
need all the SCET2 operators in Eq. ( 1D2I1 . Our interest, though, is in improving shower 
Monte Carlo, and so we only calculate operators needed for that (c/. discussion in Sec. lIIIA|) . 
To reproduce the LL emission of two gluons, the only higher dimension operator we need 
is 6*2 Lo(^i''^i)^2 • The operators O2 (no) and O2 {no,no) only tell us about the no- 
branching probabilities already determined by the one-loop cusp anomalous dimension. For 
example, O2 (no) describes a quark which does not emit until after the scale of matching ki±. 
For this reason, we call C2 lo(^15 ^^1)^2 ^^^ LO operator. Naively, two-gluon contributions 
from O2 (no) and O2 {no, uq) are lower order at tree-level, but this does not take into 
account the exponential suppression from running. The dominant contribution to showers 
comes from strong-ordering, not "every emission as coUinear as possible." Thus, we build 
our shower around C2 l^Q{ni,n[)02 ■ The coefficients C2 nlo('^1'''^i) ^^^ ^2 NNLo("'i'^i) 
give corrections for one emission. We therefore obtain a correction if we run a LL shower 
based on a matrix element computed with one of these suppressed terms. 

We now turn to calculate the terms in Eq. ( JD6I1 . in three steps: first we calculate 
the amplitudes in SCETi on the LHS of fIDSP : second we rotate it using the finite RPIi 
transformations defined in App. [B], so that the necessary operators overlap with SCET2 
states; and third we calculate the Wilson coefficients necessary to match the two sides. We 
do it order by order and we start by calculating the coefficient C2 lo- The first term of the 
LHS of dni]) is 



Al^S = U^''"''^ {no; Q,f^i)9 U [< + "^ ^; ° J <^^n , (D8) 

where f/*^^'°'°)(no) is the running factor {cf. Eqs. [181 [13 and [901) . and fii ~ XQ is at the 
scale of the emission. In fIDSI) . we have omitted the terms proportional to n" as they are 
unnecessary for matching. Gauge invariance constrains all appearances of n ■ An to come 
from the Wilson lines in x a^nd B. The amplitude is written in terms of objects projected in 
the hq and n directions. As discussed in Appendix [Bj these directions are not suitable for a 
SCET2 states, but we can use the formulas f lB7p and write (JDSP in terms of the directions 
ni and n[ where the quark and gluon have zero perpendicular component, this gives 

Y f^i Qi J % 

In (ID9P we have rotated the spinor in the ni direction, 7no± in the n'^ direction and we 
have written no in terms of ni, n'l and (gi)no±- We have dropped all the terms proportional 
to n" and we made use of relations yl^ = and ^m^i = 0. Since the gluon momentum 
is parallel to n'/*, only the polarizations in the perpendicular direction with respect to n^ 
are physical, thus we can neglect the term proportional to n'l in Eq. (1D9I) . The SCET2 
amplitude (0|xni5''B"/^±Xn|gni5'n;gn)2 is 

(0|Xni^i3"/±Xnkni^n;gn)2 = Q Un^e^.^^Vn , (DIO) 

where in Eq. flDlOp we have explicitly written the polarization vector for the gluon. From 
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Eq. (1D9|1 and Eq. f IDlOp . we can see that the LO Wilson coefficient is 



C^%o = U^'''''\no; Q. l^^YloinoY-^ll^ , (Dll) 

% 



where 



The difference with Eq. (1571) is that we replaced no in terms of external vectors. G^j [ni ■ n'^] 
is the phase space cutoff that guarantees {rii-n'^ < A^/r^"^/^ we will say more about it below. 
Since this comes from matching to a SCETi operator, (gi)„o_L ~ A and q^ ~ A^, thus C^ lo 
scales as A^^. Using formulas (IBlip . we can write fIDlip only in terms of rii and n'^, this 
gives 



CxLo = t/^''°'°H^o; Q, /ii) (|v/^^^^< + ^V^^^^ii^i<',^ 







{ni-n'^)qik 



^b'n'^-n'lfV^^^^A)^ (D13) 

where v^ is defined in Eq. f lB12|) . qi + ki = Q. For /ii, as explained in Sec. IIII At we take 
it at the scale of (A;i)„o_l as in Eq. ( 1M|) for fc = 1. Since |ff| = 2, the power counting of 
(ID13P is given by the scalar product rii -n'^, that is C(A^). In a similar way, we can calculate 

*-'2,NLO '^^^'^ *-'2,NNLO- 

We have done the matching starting from the vector current Jqcd — Ql^Q- If we had 
started from a general structure, qT^q, the results (IDlip for Cg lo would have been the 
same upon the substitution 

7^o± ^ r'^ . (D14) 

We can obtain Cg lo(^1' ^1)^2 from the SCETi operator Xno^^Xn by running down from Q 
to /ii, multiplying by the factor U^'^''^'^\nQ; Q, /ii), and subsequently using the replacement 
rule 

(Xno). ^ (c^oK)),.(XnJ,^S"'^^ ■ (D15) 

The coefficients C2 nloj ^2 nl'o' ^2 nnlo' however, are sensitive to the particular QCD cur- 
rent. This is why we refer to them as hard-scattering corrections, denoted by the superscript, 
H. 



^^ The factor of 77 ~ ^ tracks the average energy loss between mother and daughter. In choosing appropriate 
values for the parameters Sk in the numerical implementation of Q it is important to track these rj factors 
in the scaling of ni ■ n'l. 



70 



For the NLO(A) and NNLO(A) amplitudes in the second and third hne of the LHS of 
Eq. (1D5I1 we have 



ATr. = U^'^'''\no,no;Q,fii'"' 



\qqa 

^NLO 



rin — n'^ 



Q 



-9Unolno±^n 



Uq — n^ 



^?j^Lo = f^^'''''n^o,^o;Q,/ii) 



f/(2'i'i)(no,no;g, /ii) 
1 1 \ 1 



Q 9Umln[±Vn 



9Uno{klZ,±lno±Vf. 



The SCET2 coefficients needed to reproduce the amphtudes in Eq. flD16P are: 
where 



(D16) 



(D17) 
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X ©52 [rii ■ n[] 



Q 



^1^1 +gi< f^ , qih 



(DI8) 
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gi /ci / Q 



7^o± (l^l)no±7°i±- -7:: (fcl)^0±^<- 
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The coefficients scale 6*2 j^lq ~ A° and Cg nnlo ~ -^^ ^^ discussed below Eq. flS^ . we have 
a convolution because SCET fields coUinear to the same direction can exchange longitudinal 
momentum during the running. However, this convolution is only needed beyond the LL 
level considered here. 

For the coefficient Cg ^lc *^^ matching comes from the SCETi three-jet operator where 
Hi ■ n\ ~ \^/r]^}^ Since the n-labels in C| (ni,n[) are already parallel to the external 
particles, we can simply write: 



a 



2, NLO 



^(1) 



{ni,n[) = C[ '{ni,n[)Qs2[ni ■ n[] 



(D19) 



^^ With our conventions where rii-n = 2, two well seperated directions ni and n2 really do give ni -712 ~ 16. 
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where 0^2 ['^i ■''^i] only has support for (ni -n'^) > \/rj^, where it is equal to 1. Knowing that 
for this term, rii ■ n\ ~ A^/r/^, C2 isho{ni,n\) scales ~ A°, and 

Ci%Zini,n[)0^^^r^\'. (D20) 

In keeping with our conventions, we keep track of dependence on r^ ~ | for our G func- 
tions and their dot product arguments, where the various 2" factors affect where the step 
function turns over. We do not include them in the power counting for operators, where 
A parametrizes strong-ordering and the deviations from it. Accounting for rj here is cer- 
tainly possible, but in the end we always will compare amplitudes with the same number of 
external particles, so t] factors from operators will not play any role. 

The operator O2 {ni,n[) only knows that rii ■ n[ > A^/r/^, it is not able to distinguish 
its two-jet contributions, Eq. ( 1D13I) and Eq. ( ]D17p . from its three-jet one, Eq. ( 1D19I) . This 



information must then be in the Wilson coefficients, and we have put it in the functions 
9 and 6, first described in Sec. IIIIAi We can think of 0^2 [^] as usual theta function: 
QfeN = 0[62 — x] and 052k] — ^ ~ ^S2[^]: but for integrating phase space, this can lead to 
numerical problems. Instead, we can use a smoother theta function, such as the following, 
plotted in Fig. [7] 

if a; < A - a 

eA,,(x) = <| .^MpJlle^+ ^JA^^XntU + Sisn(.-A)+i ifA_a<x<A + a , (D21) 

1 a X > A + a 

The parameter A determines where the function switches from to 1, and a governs how 
fast it does it. For the SCET2 coefficients, we have A ~ 62- In order to have rii ■ n'l < A^/ry^, 
we need A^/r/"^ ^ 62 < l/?7^, so we choose 62 = A/r^"^. When we go down to lower SCETj, 
in general the Wilson coefficient has to encode the that either rii ■ ^j < X^^^"^^ /r]'^ or 
Ui ■ rij > X^^"^^^^ /r]^, in order to do so, we will use G5. where 5i = X^^~^/r]^. To see how this 
6 works, we look at the amplitude squared up to NLO(A).^^ The LO amplitude squared is 

\A'''\lo = \Ci]Unr,n[)\'G{q^,h,k2,p,), (D22) 

where 

G{qi,ki,k2,Pq) = 2{qmgn[qn\O\^^\ni,n[)\0){0\O'y\ni,n[)\qn^gn'^qn)2- (D23) 

The NLO(A) amplitude squared is 

1^ InLO ~ 1^ lNLO,2-jet T^ 1^ lNLO,3-jet; l-L^^^J 



^^ We perform some trivial azimuthal integrals in order to eliminate some terms that will drop out of typical 
observables. Also, by NLO(A) corrections for amplitudes squared, we mean suppressed by two powers of 
A. Since there are no odd powers of A in the expansion, this means NLO(A) for the cross section. 
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where 



1^ ltoNLO,2-jot 



\^ lNLO,3-jet 



— (C'2,Lo(^1''"'i)^2,n'nLo(^1'''^i) 

+ ^2,NNLo('^l' ""-1)^2,10(^1' '^l) 

+ |C2Sl'S(^l>^l)nG'(9l>^l'^2,Pg) 



(D25) 



In Fig. Owe plot the ratios \A^^^\Io/\A^'%cd and i\A^'~^\lo + l^^^"^lLo,2-jet)/l^''"'lQCD 
versus {ki)no±, We note that including NLO(A) corrections extends the region where tree- 
level SCET2 and QCD agree. In Fig. [TUl we plot the the merging of the two-jet and and 
three-jet amplitude squared using the theta function. Although we have not undertaken 
any systematic study of how our phase space cutoff enters observables, we take Fig. [10] as 
visual evidence of minimal sensitivity. Lastly, in Fig. [18] we plot |A''^^|lo 



+ I^'^^Inlo with 



Amplitude ^ 
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without running 
with running 




Gluon Perp 
oQ Momentum 

FIG. 18: Plot of the SCET2 amplitude squared up to NLO, I^'^^^Ilo + I^^^'^InlO' ^ith (green) 
and without (red) running factors versus {ki)no± ^^^ ki/qo = 0.4. 

and without running factors. As expected, the latter is suppressed relative to the former. 

2. Two- Gluon Emission 

We show the Feynman diagrams corresponding to the operators needed for matching 
two-gluon emission in Fig. [19] As discussed at the very end of of Sec. IIVAI and in App. [C] 
we do not need the two-gluon, SCETi operator, 0[ {no,nQ,no) at this order. Thus, the 
ones in (1C3P are sufficient. 
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FIG. 19: Matching SCETi to SCET2 to SCET3 for two emissions to the two-jet configuration in 
SCETi. We organize here by column number: (1) QCD Feynman diagrams; (2) SCETi diagrams 



from the operator O} (no); (3) SCETi diagrams from the operators O} (no,no) 



andT/ (no, no) 



(4) SCETi diagram from the operator 0\ (no, no, no), this operator contributes only at N^LO to 

the SCET2 matching; (5) SCET2 diagram from the operator O} (n2,n'^); (6) SCET2 diagrams 

(2) (2) (2) 

from the operators Og (n2,n'^,n2), O2 ('^2,ni,n2) and O2 {n2,n'i^n\); (7) SCET3 diagram from 

(2\ 

the operator 0\ {ni^n'i^n^)- 



The SCET2 basis has the following two gluon operators: 



Of\n,^n\^n,) 

OP{n2,n[,n2) 
OP {no, no, no) 



Xn29K{L9K2^Xn , 
Xno9B^O±9Ko±Xn ■ 



(D26) 



The last operator in (ID26I) is not necessary for the matching at NNLO(A). It can only be 
closed with states \qno9no9no'ln)2 having both gluons collinear in SCET2. Its coefficient can 
only come from the SCETi operator 0\ {no,no,no). Any contribution involving SCETi 
Lagrangian emission that matches to a higher-dimension operator in SCET2 will necessarily 
have some partons in different SCET2 directions, e.g. {no,no,n[). Since CJ (no, no, no) ~ 
A°, and the matching does not change this, O2 (no, uq, no) contributes at N^LO. The Wilson 
coefficients of the operators (ID26I) are defined such that 



•J r\ 



QCD = <[oK)Of ^ + Cf Hno, no)0« 



+ Cf)(ni,n;)0« + .-- 

Cf\no)Of + Cf Vo,no)Or + C^'|(no,no)r2^'^ 

+ C5'^(ni, n;)Oj') + Cf (n2, n;, n2)Of 

+ Cf\n2, n\,n\)Of^ + cf\n2, n\,n'2)Of^ + ■■■ 



(D27) 
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where we have written the QCD current in terms of SCETi and SCET2 operators. The 
eUipses indicate higher order terms. 

We divide the subleading Wilson coefficients in two categories: jet-structure and hard- 
scattering, labehng their contributions with the superscripts J and H. As mentioned pre- 
viously, the latter come from suppressed operators in the QCD — )■ SCETi matching and 
depend on the details of the hard partons' creation. The former are subleading terms from 
the SCETi Lagrangian that correct Eq. fISB]) as we match to lower-scale theories. They are 
completely independent of the initial hard process. 

We have seen in the previous section that the LO single gluon coefficient x operator is 
C2 {ni,n'-^)02 ~ A^, Eq. (1D7I1 . We are interested in calculating the amplitude squared 
to NLO(A). We therefore only need to calculate those NNLO(A) contributions that can 
interfere with the LO amplitude. These operators are of the form O2 {n2,n2,n'i), as the 
others in Eq. f lD26|) are not strongly-ordered. ^° 



We now calculate the coefficients in flD26|) . starting with Cg (^2,^1, rig), which we de- 
compose as: 



where 



Cf (n2, n[, n'2) = C^%ioi^2, n[, n'2) + Cf ""(na, <, n'2) , (D28) 

<^2^NLo(^2, %, n2){Q\Of\qn2gn'^gn'^qn)2, = (D29) 

CiI'loI^o) dXidX2{0\T{CsCETi{Xl)CsCETi{x2)Of'\nQ,no)}\qn2gn[9n'^qn)2 



and 

C^^^^{n2,n[,n2){0\O^2^\qn2gn[gn'^qn)2 = (D30) 

+ C[^\no,no) dx{0\T{CscETAx)O'l^}\qn^gn'^gn'2Qn)2 

+ C[]l{no,no) dx{0\T{£sCETi{x)Ti^^}\qn2gn[gn'^qn)2 

+ C[^\ni,n[) dx{0\T{CscETi{x)O^^^}\qn2gn[gn'^qn)2 ■ 

We decompose Cg (722,^,^2) as 
Cf ^^(^2, n[, n'2) = c5^)Slo(^2, n[, n'2) + C'f {^^^0(^2, ^i, ^2) + C'2^nnlo(^2, n'^, n'2) , (D31) 



where C2 nnlo(^2, ""-i, n'2) is the coefficient that reproduces the the second line in Eq. ( 1D30p . 

etc. 

(2) 
Since O2 (^2, '^i, "^2) does not interfere with the LO operator, we only need the coefficient, 

C2 NLo(^2,''^i,'"'2)- We also calculate Cg nnlo(^2,'"'i,'^2) though, because it will be useful 



^° In principle, we also have T2 {n2,n2,n'i). However, the field content alone makes this A^°, but all our 
correction operators have Wilson coefficients at 0{X^^), so its contribution is beyond NNLO(A). 
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FIG. 20: Amplitudes for two emissions in SCETi from the operator O^ 



(0) 



later. We prove below that these coefficients and their corresponding operators are of order 
A^ and A^, respectively (in Eq. \D7\ we show that LO is at A^). C*2 nloI'^s, 't-i,'"'2) and 
Cg NNLo(^25 n'^in'2) come from two-jet operators in SCETi. Thus, they both contain factors 
of 052 [n2-n'^Qs2 [^2 ■''^2]0(52 [^2 '"'^i]- We ffist described these phase space cutoffs in Sec. IIIIAl 
and made use of them in previous section on single-gluon matching. The subscript, 82-, 
constrains the argument to be < A/r/^. 

To calculate the coefficients, we proceed as with one-gluon emission: on the LHS of 
Eqs. flD29P and flD30l) we calculate the SCETi amplitude and rotate it along the directions 
77-2, n'l, n'2 where the quark and the two gluons are aligned using the finite RPIi described 
in App.dH]); on the RHS we write the SCET2 amplitude and calculate the Wilson coefficient 
necessary for the matching. We decompose the SCETi amplitude: 



.(0), 



(0)- 



^NLO = C*! ('^o) /(ia;iCia;2(0|T{/:sCETi(xi)£sCETi (3:2)01 }\qn29n'^gn'^qn)2 



(D32) 



m 



Am99 _ Amaa 1 Aiigg 1 Aiiaa 



^NLO 



^NLO.A 



NLO,B 



^NLO,C ' 



(D33) 



where A, B, C correspond to the three graphs in Fig. [201 Using the SCETi Feynman rules, 
we have: 



Amgg 

^NLO, A 



l/"o 



no 



13 , 13 



I^)no± 1,^2''"0-L f3 



qi 



Q2 



nn + 



qi 



-1; 



no- 



01% M 

„2 „2 ^noX"^" ' 



Aqqgg 

^NLO.B 



f/(2'°'°)(^o;Q,/ii)/^. 



q2 + k 



,4 + ^'^ 



^iJno-L 



no- 



AQigg 

^NLO,C 



q2 + ki 
f/(2'°'°)(no;Q,/ii)/w 



q2 + h go ^ 
{q2 + hyqf-'^^"'' 



no 



1 



■ „ „ (^2 + ^l)"0^ , (^2)»0^ .,« 

'^O + 7no± ~ TT ' = lno± 



q2 



q2 + fci^"°^^"«^ + g2 + ^2^""^^"°^. 



(D34) 



^2 7no±'^"' 

vo 



where gi = ^2 + ^2 and go = g2 + ^i + ^2- As before, we do not write terms with n" and n'^, as 
they are not necessary for the matching because the operator n-An is constrained by gauge 
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invariance to be only in Wilson lines. Now we rotate the amplitude (1D34P to the directions 
77-2 and n'l and n'2 parallel to the quark and the two gluons, as described in Eq. ( 1B14I) 



Amgg 

^NLO, A 



f/(^'°'°)K;Q,/ii)/n„,[|V;^^^2" + |v/^^^y74_ 



(D35) 



X 



?i 



h 



:7?^i « 



il / T a , "'l / Tfl a 



2 ^2 'yno±'""' ' 



Aiigg 

^NLO, B 



kn 






h 



92^2 



go 2 1 gi(g2 + Ki, 



V'^2-n'2 7°u^ 



i^2 



/Clfc2 



+ 



go(g2 + h 

q2h 



V«l-<7ni±Y 



?2 



y2 /^ , ,„ 

— V"-2-^2'^2 



~».^ ^^ /„ r>' „^ 



qo 



k2ki 



— In ■n'—i^ 



-V^i-^yT^j 



92 + /^l go ;, 



^Sc = f^^''°'°n^o;Q,/ii)/^. 



1 



g2 + fc 



ir^"i^^"2^ 



+ 



q2 + k 



rM±M± 



(g2 + fci)2g2'-o^^^' 
go 



go 



and 111 ■ n[ is defined in terms of n2, n[ and n'2 in Eqs. (IB 181) . 

The vectors Vi and ^2 are defined in Eqs. flB12p . (IB17P and (IBlSp . The values of q^, q\ 
and (gi + ^2)^ are given in Eqs. flBlOp and flB20p . As with the one-gluon emission, we can 
neglect the terms with n'" and rig as they are orthogonal to the B"^, _,_ and B^, j_ fields. The 

SCET2 amplitude for {{]\Of\n2,n\,n'2)\qn29n{9n'^qn) is: 



{^\Xn29K'^±9^i'^^Xn\qn29n'^9n'^qn) = 9 



'^n2^n{±^n'^±'^n ■ 



(D36) 



In Eq. (1D36P we have explicitly written the polarization vectors of the external gluons. For 
the jet-structure corrections, we get: 



CfNLo(^2,ni,n2) = [/(^■°'°)(no; Q, fn) d(°''^{n2,n[,n'2)<ds2[n2-n[]Qs2[n2-n'2]Q52[n'2-n[] 



(D37) 



where 



,Ja/3 



J a/3/ 



jJaPt 



d^ "^(^2, n\, rig) = d^j{n2, n[, ^2) + d^^^{n2, n[, ^2) + d^^\n2, n[, n'2) 



(D38) 
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with 



,Ja/3 



jj an / / / \ 



X 



-Qi 



^Jn2-r)!^vl + -^\fn^-'n:^ '-^In' ± 






'2 /3 

2 '"2- 



(D39) 



T?^l, 



yi / T a , '^i / r f 1 Q 



?0 



% 






,Ja/3/ / /\ yi / r a , '^2 / ^ f 2 a 

"IB (^^2,^1,^2) = — V'^l-'^l^l + — V'^2-n2— 7n', 



+ |v/^^rx|7„V + 



gi • - 2 
q2k2 



qo 



^J^^2il-Ji 



kik2 



qo{q2 + h 
q2 



+ 



X 



qi{q2 + ki. 

q2 / T p ki I J 13 

— \/n2-n2V^ - —^Jnl■n\vl 

qi _ qo 

q2k2 I —J2_ /3 

fe + ^i)g/'''"''2^«2^ (g2 + ^i)go 
^2 + fci go 



^2/^1 I — jix fi 

^- — A/rii -rii — 7 , 



7 

'71 



,Ja/9 



jj aa I I I \ 

"i,c (^2,ni,n2) 



(g2 + fci)2g2''^°^' 

1 a /3 , 1 B a 

- ^ T 7n;±7„M + -— Tr7n'±7n;± 
g2 + Kl ' ^ g2 + «^2 ^ ' 



go /. 

2 7no± • 

go 



(D40) 



The 9 functions in Eq. ( 1D37P show that C2 nlo(^2,'^i, ^1-2) comes from the two-jet SCETi 
operators. To examine the power counting of C2 ^-^Q{n2,n[,n'2), we have to consider that 
this coefficient comes from matching SCETi to SCET2 in the region where n2-n[ ~ r22 ■""-2 ~ 
n[ ■ n'2 ~ y? jif^ thus we have 



.(2)J 



C2.NLo(^2,ri'i,n2) ~ A 



-2 



(D41) 



,(2) 



and since this muhiphes O2 ~ A , by comparison with Eq. (ID 71) we see that we get an 
NLO(A) contribution. 

We proceed similarly to calculate the coefficient C\ nnlo(^2,''^i,'"'2) and show that it is 
C(A~^). We decompose the SCETi amplitude: 



.(1) 



(i)~ 



m 



^NNLO = C'l,NLo('^0,^o) / rfx(0|T{£sCETi (a;)0} }|gn2^n;^n^gn)2 



A1<199 



Aiigg I Am99 

•■NNLO ~ ^NNLO, A "•" ^NNLO, B i 

where A, B correspond to the two graphs in Fig. |211 We have: 



(D42) 
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FIG. 21: Amplitudes for two emissions in SCETi from the operator 0\ 



(1) 



AQQ99 
^NNLO, A 



-f/(^'^'°^(no,no; Q,fii)g'^Un2 



Q2 I T B 

— \/n2-n'^v'2 



, ^2 / Tf2 P 



<ll 



2 "2- 






'"'^^? Q 






gi 



(D43) 



go 



gi 



+ -\/^l-^'l^7ni± + -7-— TrTV^2-n'2 7ni±17 



A;iA;2 



qo{q2 + A;i 



/ f^ ' 

V^r<7ni±Y 



^ gi n(;-?2^^ 



7„;± 2 



gf Q 



where in Eq. (ID43P we have already rotated the amphtude to the directions n2, n'^ and n'2. 



From Eqs. f lD36P and (ID43P we can see that the Wilson coefficient Cg nnl( 



NNL0l'^2,n;,n'2) is 



Co mmt''o('^2,%,^2) =f/^^'^'°^(no,no; Q, iJi)d^"''^{n2,n[,n'2 



-'2,NNLO^ 



X 653 [n2 ■ n[] Qs2 [^2 ■ n'^] Q52 K2 " ^1] 



(D44) 



where 



,Haf: 



]HaP 



di "'^(^2, n'l, n'2) = d^ A {^2, n\, n'^) + d^ b i^2, n\, n'^) 



Ha/B, 



(D45) 



with 



^i,A (^2,^1,^2) 



92 



+ 



92 + k2 



\fn^2 



2 "^2 



q2 + k2^ '2 ^"2- 



7 



(ji rzo — n^ 



k. 



jHa/B/ / /^ yi / T a , '^2 / T ^2 a 

Lgo qi ^ 



2 



+ ^V/^^|7„V + 



go 



g2^2 



gi(g2 + ^i 



=-v/;^^7:,_ 



fciA;2 



go(g2 + ki) 



V^r^'i7n;± 



2 



7 



^ gi ra(^ - ra^ 



(D46) 
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To get the power counting of Cg nnlo("'2 5 ""-i, ^^2)) as in the previous case, we have to consider 
that the matching is done in a region where n2 ■ n[ ~ n2 ■ n'^ ~ n[ ■ n'^ ~ A^/r/^. This imphes: 

C^']Siloin2,n[,n',) r. X-\ (D47) 

which justifies its labehng as NNLO(A). 

We now turn to calculate the coefficient C2 {ni,n[,n[). We will proceed as above. We 
decompose C2 {ni,n[,n[) as: 



where 



Cf H^2, n[, n[) = C^^J/lo(^2, n\, n[) + Cf ^^(na, n[, n\) , (D48) 

C2,NLo(^2,<,O(0|O2 ("'2,<,0|gn2fl'n',fi'n',gn)2 (D49) 

= C[%o{no) dxidx2{0\T{CscETi{xi)CscETi{x2)Of\no)}\qn^gn[gn[qn)2 , 



and 



We further set: 



Ci^^^{n2,n[,n[){0\O^^^\qn2gn[9n[qn)2 
= Ci^NLoK,'^o) dx{0\T{CscETi{x)O\^^}\qn^gn'^gn'Jn)2 

+ C[]l{no,no) /rfa;(0|T{£scETi(a;)7^^^^}|gn2fi'n;fi'n;9n)2- (D50) 

Cf^"{n2, n[, n\) = C^^n^lo(^2, n[, n[) + Cf J^fLo(^2, n\, n[) , (D51) 



where Cg nlo('^2, ""-i, ""-i) is the coefficient of the contribution that reproduces the second line 

in Eq. (ID50I) . etc. We will only calculate C2 nlo('^2,'"'i, 't-i) and show that it scales as A~^. 
This is the only operator of this form that we need to calculate the amplitude squared at 
NLO(A). 

To calculate the amplitude on the RHS of Eq. flD49l) . we can use Eqs. flD35p . which are 



written in terms of n2, n[ and 112 that are parallel to the external particles, and take the 
limit n'a ■ n'l — )■ X'^/f]'^. In this case the two gluons are collinear in SCET2. Thus, we can 
define Cg ^lo('^2,'"'i,'"'i) as 

C^%ioi^2,n[,n[) = f/(2A0)(^^. Q^^^)di-P(n2,n[,n[)QsAn2 ■ n[] , (D52) 
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where 



jj cxB / I I \ 






(D53) 
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gi 



2go 



J gi g2 fci('^2-^) + ^2 hin2-n'2) 



n'2 f 7:?i. 



, ^1 / tA a , ^2^2 



A;i/i;2 



+ 



X 



+ 



go(g2 + h 
q2k2 



nvn',iZ.J-^ 



X 



\/^2-n'2 7n;. 

g2 / 

— V' 

gi 



/3 



712-712^2 



^1 / T /3 

go 



/ rf2 /3 

g2 + ki)qi " 2 "2- 



^2/i;i / rl^i /3 

^- — \/n^ ■n\ — 7 , 

(g2 + A;i)go^ ' 2 ^"2 



g2 + /i;i 



2go 



(g2 + kiY q2 ki{n2-n[) + q2 k2{n2- 



Tin 



q2 + k 



-7n'±7„u + 7 



g2 + h 



/3 a 

7n'±7„i± 



go 



q2 ki{n2-n\) + §2 /c2(n2-n2) 



X 7. 



noJ 



In Eqs. (1D53P there is a difference in the notation between the LHS and RHS. On the LHS, 
we have labeled the quark with n2 and the two gluons with n'l because the coefficient ( 1D52I1 
is for the operator O2 {n2,n'i,n'i), where the gluons are collinear. On the RHS, n2, n[ and 
n'g are the directions parallel to the quarks and the two gluons as defined in Eqs. (IBlOp and 
(jBlSp . We encode that the two gluons are collinear using the function on the RHS of 
Eq. flD53l) with ^3 = X^/t]*. It restricts that n[ ■ n'^ < X^/rj^. On the RHS of Eq. flD53l) 
we could decompose n2 in terms of n[ and avoid inserting the 9, but it is convenient to 
leave n'^ explicit because it will make the matching easier to SCET3. We notice that the 
RHS of Eq. (]D53P is just equal to the coefficient Cg nlo(^2,''^i,''^2) defined in Eq. (]D37p 



with the substitution q^ — )■ q2ki{n2-n\)/A + g2/i^2(^2'''^2)/4. Knowing that n'^ 



Uo 



xVv\ 



rti ■ n'2 ~ y? jvj^ and n\ ■ n'^ ~ A^/r^^, it is easy to check that Eq. ( 1D53I) scales as A ^. The 
information that Cg lo (^^2 , ""-i , ''^i ) comes from a two-jet SCETi operator, is encoded in the 
e-functions of Eq. ( itell . 

we decompose it as: 



(2') 
For the coefficient Cj {n2,n[,n2) 



Cf (n2,n;,n2) = cSlUn2,n[,n2) + Ci'^'' in2,n[,n2) 



(D54) 
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where 

<^2^NLo(^2,%,n2)(0|Cj^^|gn2^„'^^„2gji)2 (D55) 

= C'i°Lo(^o) dXidx2{0\T{CsCETi{Xl)CsCETi{x2)Of'^\qn2gn[9n2Qn)2 
- C2!lo(^1'^i) dx{0\T{CsCET2{x)O^2^}\qn2gn[gn2qn)2 , 

and 

Ci^^^{n2,n[,n2){0\O^2^\qn2gn[9n2qn)2 (D56) 

= C[%^o{no,no) /rfa;(0|T{£scETi(x)Of ^}|g„2fi'nifl'n2gn)2 

- C'2!NLo(^2,n'i) /(ix(0|T{£scET2(a;)C2^^}|gn2^n',^n2g^«)2 
+ Cj^|(no,r2o) dx{0\T{CsCETi{x)Ti^^^}\qn2gn{gn2Qn)2 

- <^2!nNLo(^2, %) dx{0\T{CsCET2{x)O^^^}\qn2gn[gn2qn)2 , 

We write C2 {n2^n2,n[) as 

Cf ^(na, ra'i, ria) = C'Slo(^2, ?^i, ^2) + C'2^n^lo(^2, n[, ria) , (D57) 

where C2 nnlo('^2,'"'i,''^2) is the coefficient of the contribution that reproduces the the sec- 
ond and third hne in the Eq. f lD56p . and C^ ^slq ^^^ fourth and fifth hne. As for the pre- 
vious cases, the coefficient Cg ^lo(^2 5'"'15'"'2) scales as A~^, C2 nnlo(^2 5 '^i5'"'2) as A^^ and 
C^ N3Lo("'2' '^'i-i^'i) as A°. Since O2 (^2, "^i, "^2) interferes with the LO operator, O2 {ni,n[), 
to have the amphtude squared up to NLO(A) we need both C2 ^lo ^^^ ^2 nnlo- ^^ start 
with C2 NLO- '^o calculate the amplitude in the second line in Eq. (JD55I) . we use Eq. ( 1D35P 
and take the limit n2 ■ n'2 — ?■ A'^/r^^ with n2 ■ n'l ~ n'l ■ 112 ~ }? Irf. (We could alternatively 
take the limit n2 ■ n'^ — )■ X'^ jrf with n2 ■ n2 ~ rig ■ n'^ ~ }? j-rf^j It is easy to check that 

1™,, ,f^lA^('^2,^l,^2)(0|Crkn2^ni^n2^n)2 

= C2\o{n2,n^ dx{0\T{CscET2{x)O'i^}\qn2gn[gn2qn)2- (D58) 

With Eq. (]D58|) . we can write C2 ^lo('^2,'"'i,''^2) as 



cSloK,n;,n2) = t/(2'°'°)K; Q,fii)di"^in2,n[,n2)Qs2[n2-n[] (D59) 
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where 
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The scahng of the dot products of n's in this configuration make the coefficient 



^2 NLo('^2,'"'2,'"-i) ~ A ^. As previously for Cg ^lo(^1'''^1'''^i)' ^^ prefer leaving ( ID59P 



in terms of ^2, n']^ and n'g. To calculate C2 nnlo(^2, ^2, n[) we proceed in the same way. We 
have 
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(2)H 
2,NNLO 



{n2,n2,n[) = U^^'^'°\no,no] Q, fii)d^°''^{n2,n[,n2)Q52[n2 ■ n[] 
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A;iA;2 



go(g2 + /i^i! 

In Eq. flD6ip we use the fact that, 



V^l-<7ni±lT 



7: 



^ gi n[^ -n^ 



"i^gi' Q 



n2=n2 



(D62) 



jHap, 



l(2)| 



hm <r(n2,n;,n'2)(0|Or^|g„2^n;^n2g.)2 

n2-n2— s>A*/r;'' 



c. 



(1) 

2,NLO 



(1)1 



0?x(0|T{£sCET2(a;)C2 ]\(ln29n'^9n2qn)2 ■ 



(D63) 



In Eqs. ( JD601 ID62I) there is again a difference in the notation between the LHS and RHS 
similar to Eq. (ID53I) . Since Cg lo(^2,'"'i, ^^2) and C'g NNLo("'2,'n.i,n2) come from SCETi 
two-jet operators, we include the appropriate 0-functions in Eqs. f lD59t ID6ip . 

We have that all the NLO(A) terms for two gluon matching come from the SCETi opera- 
tor, 0\ (no), and are jet-structure corrections. At NNLO(A) we have only hard corrections. 
Before matching SCET2 to SCET3, we have to insert in the coefficients the SCET2 running 
factors. Below we list all the needed SCET2 coefficients to NNLO(A) that we have calculated 
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with the appropriate RG kernels. From the matching of one-gluon emission, we have: 

% 
Ci]]S^o{ni,n[) =f/(^'^'°)(ni,n'i; /ii,/i)f/(^'^'°)(r2o, ^o; Q, fii) ® cJ'nlo('^o,^o) , 

'^2^NNLo(^i'^i) =U'''^'^'^\ni,n[; fii, fi)U^^'^'^\no,no; Q, /ii) ® cjNNLoK,no) , 



where the coefficient in ( JD64I) without the SCET2 RG- kernel is defined in Eq. ( IDllI) . the 



second and third in Eqs. flD17p . and the last in (ID19p . From the matching of two-gluon 



emission we have the coefficients: 

X d( {n2 , n[ , n^ ) Qs^ [n2 ■ n[] O^^ [^2 • n'^] Qs^ ^2 " ^1] > 
Ci%ioin2,n[,n[) = U^'^'''\n2,n[,n[; /ii,/i) ® f/(2A0)(^^. g^^^) 

X (^2(^2, n[,n[)Q52 K ■ n[]Qs2 [^2 ■ ^2] , 
cSNioi^2,n[,n2) = U^^^'^'\n2,n[,n2; /ii,/i) ® f/^'-^'^^K; Q,/ii) 
X d'^{n2,n[,n2)Q 321^2 ■ naJG^Jn^ ' ^2] > 



df{n2,n[,n2)Qs2[n2 ■ Ti'aJO^aK " '^s 



2J 5 



where the coefficients without SCET2 running are defined in Eqs. flD37t ID52t ID59t ID6ip . 



The RG kernels are given in Eqs. f lT8| [T9| and|90]). As discussed below Eq. fl89|l . we have a 
convolution because SCET fields collinear to the same direction can exchange longitudinal 
momentum during the running. 

Appendix E: Matching SCET2 to SCET3, SCETat 

We match SCET2 to SCET3 before proceeding to the general case and listing a set 
of master operators for SCETjv The SCET3 operators necessary for matching up to two- 
gluon emission are: O3 (rio), O3 (^0,^0), 0^'{ni,n[), 0^\n2,n[,n2), 0-^\n2,n2,n[), 
O3 {n2,n'i,n[). We have seen that to describe the parton shower for one emission, we 
only need the coefficient of the SCET2 operator, O2 {ni,n[). Similarly, in SCET3 we need 
the coefficient of the operator Og (n2, n[, n^). We can follow the same steps from App. iDJto 
calculate the Wilson coefficients C3 {n2.,n'^^n'2)- In this way, it is not difficult to show that 

Cf\n2.n',,n'2) = Cflo 

^{2)H,a M2)H,b ^(2) J a ^{2)Jb ^{2)J c 

"T '-^3,NLO "T "-^S.NLO "T '-^3, NLO "•" "-"S.NLO "•" "-"a.NLO 

^{2)H,a ^{2)H,fe .p.-, X 

"T '-'3,NNLO "T "-^S.NNLO ' K^^l 



where 



Ci%o{n2, n[, n'2) = C%{n2, n'2)C%{n,, n\) , (E2) 
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c, 
c. 



^{2)H,b 

'-'3,NLO 
^(2) J, 2 
'-'3,NLO 
^(2) J, 3 
'-^3,NLO 
{2)H,a 
3, NNLO 
{2)H,b 
3, NNLO 



("-2,^1 

{n2,n[ 
{n2,n[ 
{n2,n[ 
in2,n[ 
{n2,n[ 
{n2,n[ 



n 

^2 

n: 
n: 
n: 



^3,Lo('^2,'^2)C'2,NLb(^l'^l) ' 
C3,Lo('^2,'"'2)C'2,NLb('"'l'''^l) ' 

C2^NLO (^2, n[ , n'^)&s, [n2 ■ n'^] Q53 [^2 ■ ?^'i] ©53 [^2 ' ^1] 

C2, NLO (^2, < , ^1) 0<53 [^2 ■ ^2] ©53 [^2 ■ <]0<53 [^2 " ^'l] 

C'2!^NLO (^2 , n'l , ^2 ) 6^3 [^2 • n'^l Qsa [^2 ■ n[] Qs, K ■ n[] 



^(1) 



(1) 



^3,Lo("'2, ^2)C'2,NNLo(^l' ^l) 



C^NNLO (^2, n[ , n2) 653 [^2 ■ n'^l Q53 [«2 ■ ^1] ©53 W2 ■ n[] 



and 



.(1) 



^3.1.0(^2,^2) 



il2 



'ni_ 



-7„' 



92 
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9! 



2 4 



-0<53K-^2" 



(E3) 



On the LHS of the equations in the first, second and third hne of flE2l) we can write rii 
in terms of n2, n'g and n[ using the formulas in (JB16I) . The SCET2 coefficients C. 



(2)J 
2, NLO 



and C^^NNLo ^^^ defined in Eqs. (|D65p. C^^J^lol'^i,^), C'2^iSo('^i'^i)' ^"^^ C^^J^NLo("'i''^i) 
are given in Eq. flD64l) . and C^^J^o(^i'^i) in Eq. ([52]). As with any SCET^ -^ SCETi+i 
matching, we encode the definition of colhnearity from the higher scale theory in the lower 
one by G functions (c/. discussion in Sec. IIII Al) . Some of the SCET2 coefficients above 
already contained such factors as a result of matching to SCETi. In Eq. ( ]E2p . we write 

out the new ones that appear with 9^3, with ^3 = X^/r]^ according to our usual convention. 

(2) 
Since all the coefficients above multiply O^ , the scaling of contributions comes from them 



(2) 
alone, with C3 lo ~ 



A ^, the NLO terms ~ A ^, and NNLO going as A" 



At LO, the contribution in SCET3 is given by the replacement procedure on the 
LO contribution in SCET2, C2 lo(^1;'^i)^2 ■ ^^ multiply it by the running function 
U^^\ni,n[] /ii, ju) and apply the replacement: 



(x. 



-)■ 



where Clo('^i) is 



Clo(^i) 



,(92 



k. 



(cLo(^l))ii(XnJirf^- 



(E4) 



+ 



Q2 



- Qsa [^2 ■ 



Un 



(E5) 



Eq. (lE4p has the same structure as Eq. (IDlSp . If we go on with the matching down to 
SCETat, we find that the LO result would be given by applying the above replacement A^ — 1 
times. At SCETtv we could match everything to the operator O}^ ~ '{nN-i,n[, 
and the LO coefficient is 



,n' 



N-ly 



N-l 



c 



[N-l) 
N,LO 



U"" 



'-"'(nn.n;, 



fc=l 



/ 



i;/ifc_i, /ifc)cLo(^fc-i)r'' 



(E6) 
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with Hk ^ {h)nk-i± given in Eq. ([6lD and 

'k+l)no± ^ 

fcfc+1 qk+1 I 4 






where Sk = \^^~^ /r]^. 

At NLO(A), we have two kinds of corrections: hard-scattering and jet-structure. We 
notice that the NLO(A) hard-scattering terms in SCET3 are just given by those in SCET2 
with the apphcation of the replacement rule ( IE4p . If we go on with the matching down to 
SCETjv, we find that we get NLO(A) hard-scattering by applying the above replacement 
rule N — 2 times to the SCET2 hard-scattering operators. Thus, we can consider this as a 
correction to the matrix elements that we pass to a LL shower: 

^NLi"" = {Ci'f^oinu n\) + C(^}fLo(^i, <)) (E8) 

N~l 
X 

k=2 






This approach also works for hard-scattering at NNLO(A). Since we did not get C^ nisJlo 
from a replacement rule, it contains one less factor of c^q. 

(N-2 
fe=2 

where the coefficients Cg nnlo('^15'"'i) ^^^ ^2 nnlo(^2,''^i,'^2) are defined in Eqs. (]D65p . 

The NLO(A) jet-structure corrections in SCET3 are given by C^ nlo(^25 "^i, '"'2)C'3 , where 
/ = {1, 2, 3}, are given by the LO SCETi operator Xnol^Xn in three steps: First, we multiply 
it by the running factor U^^^uq] Q, fii), second, we apply the replacements 



where 



(XnJ. ^ {hf),,in2,n[,n',)ixn,),gB<^gBf^ , (ElO) 



h1^{n2,n\,n2) = d"'^ {n2,n[,n2)Qs3[^-2-^'2]^S'An2-n[]Qs3W2-^'i] ^ (EH) 

hf{n2,n[,n'^) = df{n2,n[,n[)QsA^2-n'2\QsA^2-n[]QsA^2-^i] > 
hf{n2,n[,n'2) = d'^'^{n2,n[,n2)Q52[n2-n2]Qs3[n2-n[]Qs3W2-'iT''i\ ■ 



The dj coefficients are defined in Eqs. (ID38| ID53[ ID60p . Third, we multiply the op 



erators that come from applying Eqs. (lElip by the second running factor. This de- 
pends on the SCET2 operator so each replacement rule (jElOp is followed by a differ- 
ent factor: h" by f/*^^'^'°^(n2,n'j^,n2; /Ui,/i2), /i2 by U'^'^''^'^\n2,n[,n[; ^i, ^2) and h^hy 
U^'^''^'^\n2,n'i,n2; /xi, ^2)- Since these corrections are independent of the initial hard pro- 
cess, we would encounter the same calculations we have done just now for SCETi to SCET3, 
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at any matching SCETj to SCETj+2- Thus, the NLO(A) jet-structure coefficients for the 
SCETat operator are: 



C 



{N~1).J 
Af.NLO 



N-2 

Ho. 



{N)J .,^ 
Ar,NLOl'J 



1=1 



where 



^(N-l)J 
•-"Ar.NLO 



3 l-l 

(o=5:[(n 



X U 



1=1 k=l 



TT f|(2,fc-l,0) 
k=l 



{nk-i,n^, 



, 'n-fc_i;/ifc-i, /^fc)cLo('"'fe-i, 



N-l 



n f/(^'^-^'°)K-i,n;,, 



,'^fc_i;/ifc-i, /^fc)cLo('"'fe-i. 



■p/i 



k=l+l 



with 



Ui^^\f^i, /i/+i) 



rC+i) 



f/2 (/Ui, /i/+l) 



rG+i) 



f^3 (^«> ^m) 









and 



n',,n', 



+u'H^'H+i) ^ ^1 (^m,^/,^/+i)05,+ih+r^m]05,+ih+r^/]05,+iK+i-<] 
hf{ni+i,n'i,ni+i) = df{ni+i,n[,n[_^_-^)Qsi+Ani+rn[+i]esi+A'^i+rni]Qs,^,[n'i_^^-ni] 



(E12) 



(E13) 



(E14) 
(E15) 



(E16) 



The coefficients dj here are equal to the coefficients dj defined in Eqs. (ID38| ID53[ ID60p 
upon the substitution {n2,n[,n2) — > {ni+i,n'i,n'i_^_^) and ^3 — ?■ 61+1. 



Appendix F: O(a^) Correction to Splitting Function 

One of the cross-checks on our results is the rederivation of (the abelian part of) the 0(q;^) 
correction to the q ^ qg splitting function, Pqq . This follows from obtaining the NLO(A) 
correction to two-gluon emission. For comparison, we have chosen the classic result of Curci 



>(i) 



et al. [62[. The full expression for Pqq involves many real and virtual contributions. Here 



we will only explicitly calculate the ~ C'j^ component of Pqq and show it agrees. (Obtaining 
the full result requires additional non- abelian diagrams.) Ref. 62(| splits the abelian, two- 
gluon, real emission contributions to Pqq into two topologically inequivalent diagrams, the 
box and crossed graphs. Fig. [221 We calculated each of these individually. 

The SCETi amplitude contains three graphs for two-gluon emission. These are shown 
in Fig. [20I and we give the corresponding amplitudes in Eqs. f lD34p . In order to obtain 



Pqq\ we will need to square the amplitudes and partially integrate over phase space. Thus, 
we need to choose an explicit kinematics. We redraw, in Fig. ( 123|) . our vector labels for 
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FIG. 22: Two distinct real emission contributions to Pqq drawn as amplitudes squared. They are 
referred to as the box (L) and crossed (R) contributions. 




FIG. 23: Kinematics for double gluon emission. This particular diagram corresponds to the "A" 
graph of Fig. (pOj) . 



two-gluon emission. We choose a somewhat nonstandard assignment for our variables. This 
is to aid in the comparison with [62]. The final state parton shower occurs for timelike 
virtual particles, and momentum fractions decrease the farther we are from the initial hard 



scattering. By contrast, [62|| considered a DIS-type process where the shower is spacelike. 
Since the radiation in that case comes from initial states, the momentum fractions decrease 
toward the hard interaction. Only at LO in a ,^ ar e the spacelike and timelike splitting 

At higher orders, this gets violated. 
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functions equal, by the Gribov-Lipatov relation 
but there is a straightforward conversion procedure, detailed in [62|,|78|. We, however, choose 
our kinematics such that our variable relations are equivalent to those for a spacelike process. 
For example, Pqq is a function of x = go/^2- In a spacelike process, x G [0, 1]. Rather than 
convert our answer, we will also define x as above, even though this means for us x G [1, oo). 
Other integration variables will have their ranges shifted so that they have the same relation 
with X as in DIS, and thus they enter into our expression in the same way. Lastly, we do 
not do the phase space integration for g2- While this is necessary for the timelike splitting 
function, the analogous particle for a spacelike process is a fixed initial state. Thus, for 
comparison purposes, we can leave it undone. Our vectors are as follows (note that this is 
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a different frame from tlie one used previously for matching): 

ga = {p,0,0,p} 

ki = <-Zip- -^^, ki_i_cos{(j)i), ki_i_sm{(f)i),-zip+ -^^ 
I Apzi Apzi 

fc2 = < -Z2P - , fc2±, 0, -Z2P + 



ApZ2 4:pZ2 

J , gg + |ki± + k2±P ^ -^ ql+\ki± + k2±l\ ,„,. 

go = Sa;p+ ,kix + k2±,xp >. (Fl) 

4ps 4px I 

Before proceeding, we wish to note some things about our assignment. First of all, while it 
is redundant to include go = ^i + ^2 + ^2, we will integrate over d^qo and wanted to present 
our parametrization. We see that x = \ — z\ — Z2- This is consistent with the spacelike 
case, but here, ^1,-22 ^ (~oo,0], hence the minus signs in k\ and k2. Additionally, only the 
relative azimuthal angle between kx and k2 is physical. Thus, to simplify our formulas, we 
fix ^2 in the x — z plane. 

As a last step before squaring and integrating, we will introduce our measure and inte- 
gral parametrization. While one could integrate the full final state phase space including 
the antiquark, we instead exploit the factorization of the the cross section into a hard in- 
teraction "H, a radiation-function /C, and fragmentation functions qB^rix) which determine 
how the partons arrange themselves into hadrons. Schematically, a = l-L ® {lCi^o{x,q^) + 
/Cj,NLo(a^5 g^) + . . .) ®IiqB,F{,x) =1-L® {'R-lo + "T^nlo + •••)• For our computations we need 
only integrate the phase space for 7?., and it will remain independent of the details of H. 
Taking d = 4 + e: 

^LO = E- fU^-^4%dq'FF\\Ci:~S{0\Ot%i^-l)9q)\' 



i=i 



X 






7^.,NLO = E- |i{^^^4^^dq'Vv\\Cf-^'^i{Q\Ot%{^-l)9q)? 

i i 

x6{x - goM)<5(g^ -iq^ + J2 ^^f) i'^rcYS^'\qo - J2 ^3)), (F2) 

i=i i=i 

and the qi phase space and spin-sum are moved into Ti. We define Zj analogously to Eqs. (IFip . 
The setup we describe in the body of the paper uses Wilsonian cutoffs in phase space, both to 
keep the contributions of different operators distinct via O's and to cutoff soft and coUinear 
divergences via some shower resolution parameter which keeps configurations outside of 
nonperturbative regimes. In the shower language the 11 ^^^^(x) term in TZ signifies the 
hadronization model and may depend on more than just x variables, and the /C term signifies 
the infrared finite fully differential shower computations. In Eq. flF2l) we are integrating 
over _L-momenta to carry out the perturbative comparison with Curci et al. Here we are 
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implicitly in the MS scheme, and it is the perturbative IR divergences in 7^ that get absorbed 
by qB,F{x)- The TZ terms that we need consist of only the 1/e^ and 1/e portions of the 
corresponding operator expectation values. The non-pole contributions from C^ lq ^^^ 
W*NLO (Eqs. IE6l and [E13l) . along with higher-order corrections are in higher order terms in 
the TZ functions. The hard-scattering corrections are in "H. The reason we extract only the 
pole terms is that these are precisely what give the expression for Pqq and Pqq ■ In addition 
to selecting the pole part, we also define PP to remove those portions of the matrix element 
which enter into H such as the final quark spin-sum, current P, and antiquark quantities. 

In this MS factorization scheme, we need to define our correction operator differently 
than in Apps. [D] and [El Since Pqq requires the calculation of two-gluon emission, we find 
it simplest here to calculate in SCET3 where only C^ ^lo i^ ^^- ^E2[) contributes. This 
corresponds to taking limits such that only its G-function equals one, while the other jet- 
structure coefficients are zero. Since we integrate it over all of phase space, which includes 
the strongly-ordered limit, we need to subtract the LO contribution. This just comes from 
C3 l^Q{n2, n[, 112)0^ , but we take care to only remove the pole parts consistent with MS. 
We can thus write the subtraction as: 



'v 7 MT n / ^^-LJ-Am A:o rtn-t^-T 



J, NLO ~ / "'^^ki,k2,qo 



\Ci%fo{n2,n[,n',){0\Of^\qggq)\' 



MS 



(F3) 



where C^ j^lo ^^ evaluated such that = 1 over all of phase space. We will describe the 
subtraction portion in detail below, but first we concentrate on the correction term. 

By fixing the virtuality of ^q = g^, we can obtain an expression without having to know 
its exact limits, which will depend on the details of the hard scattering. For Pqq , one only 
needs to calculate one-loop corrections to single emission and tree-level double emission, and 
we now specialize to the latter case. We perform the d-dimensional integration over d'^qo 
and rewrite the integral in terms of kij_ and k2± dependent functions with zi 2-dependent 



coefficients. Using the same parametrization as Ref. 79|, we can write 

dzi dz2 d'^-^ki^ d'^-^ 
(IGtt^)^ J ^ Zi Z2 IT vr 

x6{q'^ - (ai ki^ + 02 kal - ki_L ■ ka^ )) 

1 Z'./ N , n/ Nki^-ka^ , ^, ,ki_L-k2 



^rNl^S = 7T7r^/^.^--^^^^^^(l---^i-^.) (F4) 



X—{A{Z^,Z2)+B{Z,,Z2) " ^ +Ciz^,Z2) ^ 2 

r V ki^ ka^ 

+ D{z,,Z2)^^^^^^+ E{z,,Z2)^ + F{zuZ2)^) -[LO], (F5) 

where oi = —(1 — Z2)/zi and 02 = —(1 — Zi)/z2. The functions A, B, C, D are defined in 



79( 1 ■ and their corresponding kij^ integrals are finite. We can check the intermediate step of 
their integration with 79(|. The terms in our q^ 5- function have a relative sign compared to 



theirs, as our q^ > 0. As a computational aside, we found it easiest to pass to a change of 
variables: {u = ki±k2±, w = ki±/k2±). Then the (5-function just enforces: 

U = Uq = rTT- {^^) 

aiw^ + a2 — 2w cos((pi) 
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Performing all but the dzi integrals in TZ, we get Table UTTl which corresponds to [79|'s Table 
5. 



Function of ki_|_ 
in integrand of equation (JF4|) 


Contribution to 7Z multiplying 
(16^2)2^ / dzi dZ2 6{l -zi- Z2-x) 


1 


A{zi,Z2) 


kl_Lk2X 


-^^B{Z^^Z2) 


kixk2j_ 

k2^ 


-^C(Z1,^2) 


(kix-k2x)^ 

klik2l 


(1+21^1^ 


{1-Z^){1-Z2) 


)D{zi,Z2) 



TABLE III: Purely finite contributions to IZ 



We thus reproduce the earlier result. 

The E, F functions multiply integrals that lead to single e poles after the dki±^ integrals 
(and double poles after integrating g^), and so we must be more careful in treating them. 
These double poles correspond to the LO contribution, which we are explicitly subtracting 



,(1) 



as it does not contribute to Pqq ■ We discuss the subtraction in detail below Eq. ( ]F12p . 



For now we concentrate on the divergent integrals multiplying E and F. When we did our 
computations for Table (1III|) . we were helped by the finiteness of the expressions under the 
dkix integration. We could thus take e — )■ for these terms, which greatly simplifies their 
integrals. By contrast, we will need to keep the e-dependence of the E, F terms, which 
results in an intractable computation. To get around this, one can introduce subtraction 
functions, which simply reproduce the e poles (these are merely a computational aid and are 
not related to the subtraction of LO). We will need to take care that they do not remove any 
finite pieces. Secondly, since their full contribution to 7?. is oc 1/e^, we will need to include 
for E and F any terms oc e that multiply -^ or -^. These arise from doing Dirac algebra 
in (i-dimensions. 

To do the integrals in TZ which multiply E and F, we will change variables to m, w, and 
perform the u integration as well as the trivial 02 azimuthal one. We get for this contribution 
to7^: 

TAe^f = TTTT^^- I d(f — - — -d(j)idw5{l- X - Zi- Z2) 



(167r 



2^2 



TT 



wu: 



2+e 




2g2 



-E{Z^,Z2) + 



u 



2+e 



2w q- 



:nzi,Z2)]-, 



(F7) 



where Uq is defined by equation flF6p . We only need the leading poles in e, and so rather 
than performing the w and 0i integrals for the functions multiplying E, F, we will define 



subtraction functions to reproduce the poles of 



,2+s 



2g2 



"0 
2wq^ 



respectively: 



Sf 



S, 



3L 

2al 
2al 



w 



;«; + !)' 
w + w^) 



(F8) 
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Function of kij^ 
in integrand of equation ()F4I 



Contribution to TZ multiplying 

J dzi dz2 5{1 — zi — Z2 — x) 



(167r^) 






2x zi 
(\-z-iYz2 e 



1-eln 



1-Z2 
21 



+ 



22(1-22)^ 



212:2 +x (in 



22(1-22) 



2:1a; 



E{zi,Z2) 



K2j 



23:2:2 



{1-ziyzi € 



1-eln 



1-21 



22 



+ 



21(1-21)^ 



21^2 + 3; (in 



2l(l-2l) 
22 a; 



Fizi,Z2] 



TABLE IV: Contributions to TZ\e f 



Integrating these in w gives us a pure 1/e term. Subtracting them from the functions in 
equation flF?^ : 



Ae 
Af 



g2y^3 



2q2 



2w q^ 



2(02 + aiw'^ — 2wcos(0i))^ \a2 + aiw^ — 2wcos(0i) 



w q 



2w{a2 + aiw'^ — 2w cos(0i))2 V'^2 + aiw'^ — 2w cos(0i) 



w q 



(F9) 



leads to finite integrals, allowing us to pass to the e — )■ limit prior to integration, making 
the calculation tractable. After integrating w and (pi, we want the e"^''' pieces as these turn 
into the single and double poles upon doing the g^ integral and contribute to TZj^^lo- The 
e° piece has one contribution besides that from {Ae,f — SE,F)\e=o {<Se,f contributes a pure 
1/e pole). Our w integration goes from to 00, and we obtained Se,f by expanding Ae,f 
in the appropriate w — ;■ 0, 00 limit to pick up the pole, while carefully regulating the other 
integration limit so as not to contribute its own spurious divergence or any subleading terms. 
However, we see that in equation f lFOp . taking these limits actually results in factors (aiw)~^ 
and {w/a2Y- Expanding the af" to LO in e does not affect Se,f- Nonetheless, since the 
subtraction functions have 1/e poles, including the NLO part of the e-expansion will yield 
an e° contribution. This 0{e^) term is not in AE,F\e=o since they send Mq -> 1- Thus, we 
have the following addition to the contributions from the integration oiTZ\E,F'- 



Be 
Be 



-eln(ai) 



-eln(a2) 



q^ 



w' 



2al{w + l)' 
,2 



Q 



w 



2al{w + w^)' 



(FIO) 



^'° contributions after w and 0i integration come from: Se,e + Be,f + 
Q. For integrating the first two terms, we leave the full e dependence as this 



In the end, our e 

{Ae,f — <Se,f 

was tractable. Collecting everything, we can obtain the counterpart to Table [XTTl for E, F, 

(Table lYD- 

Having set up this much of the integration, we can take the amplitude squared from the 
process of interest and decompose it in terms of the A{zi, Z2), B{zi, Z2), etc. basis. We then 
simply have to read off the results from Tables [XTTl and HVl and perform the zi^2 integrals. 
One of these is made trivial by the remaining x-dependent (5-function. As mentioned at the 
beginning of this Appendix, |62| recognizes two topologically distinct contributions, which 
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we shall refer to as box and crossed (c/. Fig. [22]) . because of their appearance as cut two- 
loop diagrams. We can identify them in our calculation by their color structures {Cp and 
Cp — \Cp Ca, respectively). In fact, we can already calculate the entire crossed contribution 
as it only involves terms from Table llllt having no double pole contribution to 7^ and thus 
requiring no subtraction of LO. Determining the box graph, however, involves treating the 
LO subtraction properly. 

As this subtraction is one of the more subtle points of the computation, we will present 
it in some detail. Its handling is tied up with what one means precisely by a "subleading 
splitting function." At LO in a^, the definition is clear. The same splitting function that 
gives us the probability for a 1 — > 2 radiation also determines the running of parton densities: 

Q'^/(^, Q') = ['^ P^^ (f' ^s{Q')) f{z, Q% (Fii) 

where the 0{as) part of Pqq, Pqq is given by Eq. (jS]). To determine Pqq\ we have had 
to calculate a 1 — ?■ 3 splitting, thus the probabilistic interpretation in terms of radiation is 
nontrivial as it involves a mix of 1 — )■ 2 and 1 — ?■ 3 processes. At the level of Eq. fIFlip though. 



we see that we are just correcting PDF evolution. In addition to the real-emission calculation 
that we are pursuing, one can alternatively determine Pqq from the anomalous dimension of 



certain twist-2 operators |80L l8l[|. Ref. [62| made a comparison to this approach and found 
agreement to 0{a1). Since Pqq is thus a two-loop object, it has the scheme dependence 
one would expect at this order, and so we need to make sure that we compute in the same 
scheme, which is why we do our LO subtraction in MS. In SCET, one could attempt the same 
cross-check from a straightforward two-loop calculation after fixing to one's renormalization 
scheme of choice. 

We will now show how to subtract the LO portion in the calculation of T^jnlo- We get a 
double coUinear pole associated with the strongly-ordered emission of two gluons. We want 

(2') (2") 

to write this as removing the emission coming from our LO operator, Cg lo^s ■ ^^ with 
any subtraction scheme, while the pole is unambiguous, we need to make sure to remove 
the appropriate finite pieces. We note that c^q defined by Eq. fl57j) contains NLO(A) pieces 
(in SCET3 power counting) which come from the offshellness of the intermediate quark. It 
is true that the LO replacement rule, Eq. (l56l) . gives only the splitting function times the 
logarithmic, collinear divergence. Nonetheless, the Wilson coefficients given by Eq. (163|1 for 
offshell quarks have additional terms. From the point of view of amplitude matching, this 
poses no problem. However, if we want to copy [62] 's scheme, then we can only subtract poles 
associated with the pure LO result after integration. As an operator subtraction in SCET3, 
this means we need to change C3 ^q. In order to recover the correct splitting function with 
no NLO contribution, we will need to project the offshell quark momentum to an onshell 
one with the same p-fraction. This alone, though, does not specify the spatial orientation 
of the vector and will not necessarily kill the subleading terms. To do that, we write the 
replacement rule, but in the limit that the offshell quark's daughters are exactly collinear 
with it. Equivalently, if we are in the frame determined by ra = {1, 0, 0, —1}, we can project 
the quark momentum along n = {1, 0, 0, 1}, i.e. qi — )■ ^n = q[. Since the replacement rule 
also makes reference to the quark's parent's momentum, we also need to project it to what 
it would be if it had emitted an onshell quark with q'^. Thus, gj_i — > ki + q[ = g^_p In the 
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end, this changes our replacement rule coefficient for the j^^ quark to: 



-LO 



"i+i 



nj,+ 



where g'^ = g^ n'^,/2. Thus Clq has the same form as clo but with a different orientation for 

its momenta. This changes the expression for Cg lo^s ^^ involve c'lq instead of clo (c/. 
Eq.[63D 

After the dq'^ integration, the 1/e term in 'R-'j'^lo will allow us to read off Pqq ■ As a re- 

minder, we need this subtraction operator because our NLO(A) term, C3 nlo('^2,'"'i,''^2)C'3 
is supported over all of phase space, and thus contains LO portions. We therefore have 



'^J, NLO 



dIiki,k2,qo^^ 



\Ci%il,{n„n[,n',){0\O^^^\qggq)\' 



-0 



LO 






Oi2 I 



LO 



,(2) 



O'.'T^lqggq)]' 



MS 



(F13) 



The MS indicates that we are only subtracting pole parts of the LO contribution with no 
finite pieces. However, there is still an ambiguity over which pole parts we subtract, since 
the LO contribution has a double pole from its two coUinear divergences, but we are at some 
liberty to decide which single pole parts we remove as well. As we expect, this subtraction 
operator squared takes the form of a convolution of two splitting functions: 

„2 



dU 



-^LO 



)(2) 



io'"'{0\Ofnqqgg)\' 



MS 



g2 ^ill) dq^ dyxp 



y)Hq 



-l+e/2 



a 

2^ 



X -P. 



y 



(0) ^ . Pi,qq(^/y) 



d,qq 



Tr 



g2fmt 



(F14) 



where the trace contains those terms that get passed to the hard function, "H, along with the 
g2 phase space by the projector PP. This includes the final quark spin-sum and phase space, 
the current F which is a spectator for both LO and jet-structure corrections, and quantities 
related to the antiquark (c/. Eq. |27]). What may seem surprising is that the two splitting 
functions live in different dimensions. The reason for this particular scheme for regulating 
phase space has to do with the alternate, two-loop method for calculating Pqq , which was 
the original approach. For that result, in MS we would subtract a pure pole counterterm, 
regulate the loop integral in li-dimensions, and leave external particles in 4d. Since the phase 
space integrals are related to loops by cuts, we see above that our y-mtegral is, in fact, in 
(i-dimensions, but the splitting involving two external particles is left simply in four. 

Looking at the SCETi diagrams for the process (Fig. |20]), the amplitude 



-LO 



-LO 



02 



mo. 



(2) 



T^lqggq) comes from a subset of diagrams A^ and B^. The expression 



for subtraction is thus: 



PP 



dn{\c'^o'c'^o'{0nT\qqgg)\ 



MS 



dq'^ dzixp 



Zl 



(g2)-i+^/2 «2 



1 



27r2 X + Zi 



\ x+zi I 



X + Zl 



i'^ 



X 



X + Zi 



X + Zi 

l + {x + Zxf 
X + Zi — 1 



+ Zi ^ Z2, 



(F15) 
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where we now act with PP, dropping the trace from Eq. (1F14P and keeping only those terms 
needed for the computation of 7^ and Pqq . We can note several things about this expression. 
For concreteness, we discuss the ^i-dependent term corresponding to graph A^, Fig. [201 The 
p fraction of go relative to gi is x/{x + Zi), and that of gi to g2 is a; + ^i, in terms of the 
variables in Eq. ( ]F14p . y' = x/{x + zi). Performing the integrals leads to double and single 
poles. For later use, we write down the result of doing the dq^, dzi integrals, where one of 
latter is trivial since we have 5{1 — x — Zi — Z2) sitting inside the phase space integral (c/. 
Eq.EH). 



2xp 



a 

2^ 



rfnpp 

1 2 



c'Lo^4o"^(0|OrTlgg-^^)| 



(2)n 
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e^ X 



-[[-2[2{x' + l)\og{\) + {x-lf) 
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+ (x" - 1) Li2 
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(x - If) log(A) - (3x2 + 5) (a; _ 1) log2(x) + 2(x - 1)^ log(x)) 
+6x (x^ + 1) (x - 1) log2(x - 1) - 2x(x - 1)^ 

-2x(x + l)(x - Vf log(x - 1) log(x)((x - ifx) 



(F16) 



where we have done the dzi integrals between 1 — x + A and —A to regulate soft divergences. 
All A-dependence cancels out of the final answer, which gives us a consistency check on the 
scheme. 

Before comparing Pqq , we can check our setup with Pqq\ by looking at the 0{as) con- 
tribution to IZi^o We see that [62| gets the following contribution: 



(0) ^ (as\ 2 l+x^ 
"« V27ry el-x' 

Calculating in SCETi, we get the following amplitude squared: 

2k - -2 



A, 



1^Q9 



With our definition of T^lo in Eq. (1F4I) . we get: 



go V ki qih qi J 2 



2 1 + x2 



(F17) 



(F18) 



p(0) 



m 



a. 



27r/ e X 



(F19) 



The overall minus sign between Eqs. (1F17I) and (1F19I) is due to the difference between the 
spacelike and timelike processes. It arises in the dzi integral. Even though the Zi dependence 
is the same in the two calculations, and the integration limits are the same, and 1 — x. 
For us, 1 — X < 0, but in [6^, it is positive. 
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We will compare the different contributions to double emission separately. In SCETi, 
the C graph in Fig. [20] will give box and crossed terms when interfered with itself and the 
A and B ones. We identify the crossed contribution by inserting the color structure and 
taking those terms proportional to Cp — \Cf Ca- As mentioned above, it only contains the 
integrals in Table IIIII In terms of its notation, we have: 



Function defined in Eq. (|F4p 


Value in crossed diagram 


A{zi,Z2) 


16x{x^+xzi + {zi-l)zi+l) 
zi(x+zi-l) 


B{zi,Z2) 


8{x'^{zi-2)-xzi+zi-l) 


Cizi,Z2) 


8{x{x^+{x-l)zi+2)+zi) 

2l 


D{Z,,Z2) 


16 (x2 + 1) 



TABLE V: Contributions to crossed amplitude squared diagram 



The box contribution additionally contains the functions in Table IIV[ though we are only 
interested in the finite parts. Their Zi dependence is: 



Function defined in Eq. ()F4 



A{zi,Z2) 



B{zi,Z2) 



C{zi,Z2) 



D{zi,Z2) 



E{zi,Z2) 



F{zi,Z2) 



Value in box diagram 



12x^ + 8xzi + 8{zi - l)zi + 12 



8(zi-l)(x^ + (zi-2)^i+2j" 



X + 21 —1 



8{x+zi)[2x^ +2xzi+zl+l) 



.^L. 







4 

+4e 



2x^+6x^ zi+x'^ [7zl+2)+2x(2zl+zi)+zf+z'l) 



(x^{x+zi-lf+z'f{{x+zi-lf+x+zi)+xzi{x+zi~l)'^) 



{zf-2zi+2){x'^ + {zi-l)'^) 
(x+zi-iy^ 



+4e 



{x'^{(zi~l)zi+l)+x(zi-l)({zi-2)zi+2)+(zi-l)'^{{zi-l)zi+l)) 
(x+zi-iy^ 



TABLE VI: Contributions to box amplitude squared diagram 



For the crossed contribution, we perform the multiplication in Table IIIII with the functions 
defined in Table IVJ and integrate dzi, having already done the trivial dz2 integral. We again 
use a cutoff to avoid soft divergences, thus its range is between 1 — x + A and —A. In the 
end, we obtain: 



P 



(1) 



qq crossed 



2^ 



1 + x' 



X 



(41n(x - 1) - In^(x) - ln(A)) - 2(x + 1) ln(x) 



(F20) 
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The A-dependent pieces will cancel against those from the box contribution. The other terms 
agree with [62|] up to the previously discussed minus sign, and wherever ln(l — x) appears 
in the spacelike calculation, we get ln(a: — 1). Since our integrand and integration region are 
real, the imaginary pieces generated by ln(l — x) when making a; > 1 all must cancel. 

The box calculation proceeds similarly using the functions defined in Tables IIIII and IVII 
We also include the terms proportional to E{zi,Z2) and F{zi,Z2) and we have subtracted 
the appropriate contribution, Eq. (1F16P from that given by C^ ^lo- Doing all this, we get: 



qqhox 



2nJ 



1 + x' 



X 



(ln(A) - ln(a: - 1)) + 2(2s - 1) ln(s) 



(F21) 



The soft divergent pieces cancel against the crossed contribution, and once again we agree 
with 62| up to an overall sign, and the continuation ln(l — x) — >■ \n{x — 1). 



In addition to these real emission contributions to the Cp portion of Pqq , there are also 
single-emission, one-loop diagrams, shown in Fig. [2H We can account for their contribu- 






FIG. 24: Single emission, one-loop contributions to Pqq ■ 



tions in SCET easily. We have already derived the tree-level expression for single emission 
(Eqs. IF17I and IF19|) . Furthermore, both the quark wavefunction renormalization and the 
vertex renormalization are the same in SCET as in QCD 21|. Thus, we recover the entire, 
gauge-invariant, oc a^ Cp contribution to the splitting function, in agreement with Ref. [62J, 
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Here we have written Pqq with its usual sign conventions for spacelike evolution. 
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